Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,772,467 C→T 6.3% L63L (CTG→TTG)  mtlA → mannitol‑specific PTS enzyme: IIA, IIB and IIC components

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,772,4670CT 6.3% 196.4 / 2.6 63L63L (CTG→TTG) mtlAmannitol‑specific PTS enzyme: IIA, IIB and IIC components
Reads supporting (aligned to +/- strand):  ref base C (26/33);  new base T (2/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

GTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGGGCGGCTGGTGCA  >  NC_000913/3772321‑3772611
                                                                                                                                                  |                                                                                                                                                
gTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCgctgct                                                                                                                                                 >  7:51674/1‑148 (MQ=255)
     cTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCgg                                                                                                                                           <  7:37710/149‑1 (MQ=255)
      tCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGt                                                                                                                                          >  5:78055/1‑149 (MQ=255)
        aGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGt                                                                                                                                          <  4:305086/147‑1 (MQ=255)
        aGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGt                                                                                                                                          <  6:307736/147‑1 (MQ=255)
        aGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTta                                                                                                                                        <  8:51674/149‑1 (MQ=255)
        aGCAACATGGTGATGCCAAATATCGGCGCGTTTAGCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGGTACCGAACGAGACGCTGGCGAAGCGGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGt                                                                                                                                          <  3:71677/147‑1 (MQ=255)
           aaCATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATAc                                                                                                                                     <  3:13479/149‑1 (MQ=255)
             cATGGTGATGCCAAATATCGGCGCGTTTATCGCGGGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCg                                                                                                                                   <  1:241356/149‑1 (MQ=255)
                 gTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCggtg                                                                                                                                <  4:35092/148‑1 (MQ=255)
                       ccAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGACAGCGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCt                                                                                                                         >  4:75200/1‑149 (MQ=255)
                        cAATTATCGGCGCGTTTATCGCGTGGGGTATCATCAACGCGTTATTTATTCCAACAGGGTTGTTACCTAACGAGACGCTGGCGAAGCTTGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATAc                                                                                                                                     <  8:168897/136‑1 (MQ=255)
                           atatCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTgg                                                                                                                       >  2:119368/1‑147 (MQ=255)
                              tCGGCGCGGATATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGAGAAGCTGGTCGGCCCGATGATAACTTAACACCTGCCGCTGCTGAACCGTTATACCGGTGGTAAGCTCGAAgga                                                                                                                  >  4:204426/1‑148 (MQ=255)
                                    cgTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGaa                                                                                                            <  6:17711/149‑1 (MQ=255)
                                         aTCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCggtg                                                                                                                                >  4:267516/1‑124 (MQ=255)
                                         aTCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCggtg                                                                                                                                <  3:267516/124‑1 (MQ=255)
                                          tCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTggc                                                                                                      <  1:119368/149‑1 (MQ=255)
                                                    gttcatcaCCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGGGTAGTCg                                                                                            <  3:204426/147‑1 (MQ=255)
                                                      tcatcaCCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGt                                                                                          >  7:278182/1‑149 (MQ=255)
                                                      tcatcaCCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGt                                                                                          <  8:290399/149‑1 (MQ=255)
                                                        atcaCCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGc                                                                                        >  8:124402/1‑149 (MQ=255)
                                                           aCCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGt                                                                                                 <  2:226792/137‑1 (MQ=255)
                                                           aCCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGAGt                                                                                                 >  1:226792/1‑137 (MQ=255)
                                                            ccGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCAt                                                                                     >  5:182905/1‑148 (MQ=255)
                                                             cgcgTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAggcggc                                                                                                               <  2:48845/121‑1 (MQ=255)
                                                             cgcgTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAggcggc                                                                                                               >  1:48845/1‑121 (MQ=255)
                                                              gcgTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTATGCGGCGAACGTGGCGGCGTAGTCGGTGCCATcac                                                                                  >  2:138107/1‑149 (MQ=255)
                                                               cgTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGACCGTGGCGGCGTAGTCGGTGCCATcacc                                                                                 >  5:56761/1‑149 (MQ=255)
                                                                 ttatttatTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGACGGCCCGAAGAACACTTATCTCCTGCCCCTGCTGATCGCTTATACCGGTGGTAAGCTGGGAGGAGGCGAACCTGGCGCCGTAGTCGGTGCCAAcaccac                                                                               >  3:65195/1‑149 (MQ=255)
                                                                  tatttattCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTa                                                                                                                             >  1:66859/1‑102 (MQ=255)
                                                                  tatttattCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTa                                                                                                                             <  2:66859/102‑1 (MQ=255)
                                                                    tttattCCAACAGGGGGGTTACCGAACGCGACGCTGGCGAAGCTTGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCAccacca                                                                             <  6:108222/148‑1 (MQ=255)
                                                                        ttCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGc                                                                        >  8:63422/1‑149 (MQ=255)
                                                                             aCAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTAt                                                                   <  6:78055/149‑1 (MQ=255)
                                                                                     gtTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGAAAACTTATCTCCTGCCGCTGCTGATCGGTTATACAGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGgcg                                                           >  8:425329/1‑149 (MQ=255)
                                                                                          cGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGa                                                       <  1:138107/148‑1 (MQ=255)
                                                                                            aaCGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACAt                                                    >  4:177050/1‑149 (MQ=255)
                                                                                             aCGAGATGCTGGCGATTCGGGTCGGGCCTATGATCACTTATCTCCTGCCGCTGCTGAGCGGTTATACCGGTGGTAAGCTGGGAGGCGGCGTACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGAct                                                     <  6:56761/147‑2 (MQ=255)
                                                                                             aCGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGTTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCAccacca                                                                             >  4:29829/1‑123 (MQ=255)
                                                                                             aCGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGTTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCAccacca                                                                             >  4:24416/1‑123 (MQ=255)
                                                                                             aCGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGTTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCAccacca                                                                             <  3:29829/123‑1 (MQ=255)
                                                                                             aCGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGTTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCAccacca                                                                             <  3:24416/123‑1 (MQ=255)
                                                                                               gagaCGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACa                                                     >  5:325508/1‑145 (MQ=255)
                                                                                               gagaCGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACa                                                     <  6:325508/145‑1 (MQ=255)
                                                                                                  acgcgtgCGAAGCTGGTCGGGCCGATGCTCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGAt                                              <  3:75200/143‑1 (MQ=255)
                                                                                                     ttGGCGCAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGGCGGCGCAGACATGCCGATGt                                            <  2:373053/147‑1 (MQ=255)
                                                                                                     cTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGtt                                           >  1:373053/1‑149 (MQ=255)
                                                                                                        gCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGgcgc                                                          >  2:367205/1‑131 (MQ=255)
                                                                                                        gCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGgcgc                                                          <  1:367205/131‑1 (MQ=255)
                                                                                                        gCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCt                                        >  2:356234/1‑149 (MQ=255)
                                                                                                         cGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTc                                       <  7:63422/149‑1 (MQ=255)
                                                                                                             gCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGt                                                                                          <  7:90022/94‑1 (MQ=255)
                                                                                                             gCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGt                                                                                          >  8:90022/1‑94 (MQ=255)
                                                                                                                             atCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTc                    <  3:177050/148‑1 (MQ=255)
                                                                                                                              tCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAg                       <  1:329927/144‑1 (MQ=255)
                                                                                                                              tCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAg                       >  2:329927/1‑144 (MQ=255)
                                                                                                                                  ttATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTgg              >  1:242668/1‑149 (MQ=255)
                                                                                                                                     tctcCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGggc            <  2:281843/148‑1 (MQ=255)
                                                                                                                                     tcgccTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGggc            <  2:281846/145‑1 (MQ=255)
                                                                                                                                          tGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGGGCGGCTgg      >  7:11036/1‑149 (MQ=255)
                                                                                                                                            ccGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGGGCGGCTGGt     <  8:11036/148‑1 (MQ=255)
                                                                                                                                             cGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAg                       >  6:140528/1‑129 (MQ=255)
                                                                                                                                             cGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAg                       <  5:140528/129‑1 (MQ=255)
                                                                                                                                              gctgctGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGGGCGGCTTGTGCa  <  7:425329/149‑1 (MQ=255)
                                                                                                                                              gctgctGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGGGCGGCTGGTGCa  >  5:219508/1‑149 (MQ=255)
                                                                                                                                              gctgctGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGGGCGGCTGGTGCa  >  5:192903/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTTTCCTCAGCAACATGGTGATGCCAAATATCGGCGCGTTTATCGCGTGGGGTATCATCACCGCGTTATTTATTCCAACAGGGTGGTTACCGAACGAGACGCTGGCGAAGCTGGTCGGGCCGATGATCACTTATCTCCTGCCGCTGCTGATCGGTTATACCGGTGGTAAGCTGGTAGGCGGCGAACGTGGCGGCGTAGTCGGTGCCATCACCACCATGGGCGTTATCGTCGGCGCAGACATGCCGATGTTCCTCGGTTCTATGATTGCAGGTCCGCTGGGCGGCTGGTGCA  >  NC_000913/3772321‑3772611

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: