Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,968 G→C 6.2% G404G (GGC→GGG lacY ← lactose permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,9680GC 6.2% 546.7 / 19.5 162G404G (GGC→GGGlacYlactose permease
Reads supporting (aligned to +/- strand):  ref base G (77/74);  new base C (3/7);  total (81/81)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 171 total aligned reads displayed.

TGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGC  >  NC_000913/361824‑362116
                                                                                                                                                |                                                                                                                                                    
tGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggg                                                                                                                                                 <  1:1796051/150‑1 (MQ=255)
 gCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAAGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggggc                                                                                                                                               <  1:1516197/150‑1 (MQ=255)
          tATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGc                                                                                                                                       >  2:261412/1‑150 (MQ=255)
          tATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGc                                                                                                                                       >  2:2412138/1‑150 (MQ=255)
          tATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGCCCCGGCCCGCtttcc                                                                                                                                       <  2:1372384/150‑5 (MQ=255)
           aTTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCg                                                                                                                                      <  2:477808/150‑1 (MQ=255)
           aTTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGCCCCGGCCCGCtttccc                                                                                                                                      <  1:713350/150‑6 (MQ=255)
           aTTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCGCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGCCCCGGCCCGCtttccc                                                                                                                                      <  1:529070/150‑6 (MQ=255)
            ttCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGt                                                                                                                                     <  1:190461/150‑1 (MQ=255)
            ttCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGt                                                                                                                                     >  1:1661150/1‑150 (MQ=255)
            ttCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGCCCCGGCCCGCtttccct                                                                                                                                     <  2:455636/150‑7 (MQ=255)
             tCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTg                                                                                                                                    <  2:2442663/150‑1 (MQ=255)
             tCTTTCCGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTg                                                                                                                                    >  2:2371067/1‑150 (MQ=255)
                ttCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAac                                                                                                                                 >  1:1666462/1‑150 (MQ=255)
                      cATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGaa                                                                                                                           >  2:120118/1‑150 (MQ=255)
                        ttGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAAt                                                                                                                         <  2:1559195/150‑1 (MQ=255)
                          ggCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGCAAATTa                                                                                                                       >  1:767185/1‑150 (MQ=255)
                           gCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTaa                                                                                                                      <  1:1281993/150‑1 (MQ=255)
                           gCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTAGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGAAGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTaa                                                                                                                      >  2:134373/1‑150 (MQ=255)
                            cATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAg                                                                                                                     >  2:68502/1‑150 (MQ=255)
                             aTGTTCAATGCGATGACTCCTTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGTGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAgg                                                                                                                    <  2:2435151/150‑1 (MQ=255)
                             aTGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAgg                                                                                                                    <  2:1093162/150‑1 (MQ=255)
                               gTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTg                                                                                                                  >  1:128497/1‑150 (MQ=255)
                               gTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTg                                                                                                                  <  1:192542/150‑1 (MQ=255)
                                  cAATGCGATCACTCCGTTATTATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAg                                                                                                               <  2:94118/150‑1 (MQ=255)
                                  cAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAg                                                                                                               >  2:1210034/1‑150 (MQ=255)
                                      gCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCa                                                                                                           <  2:1685680/150‑1 (MQ=255)
                                      gCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCa                                                                                                           <  2:378486/150‑1 (MQ=255)
                                       cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                                                                                          >  2:1308933/1‑150 (MQ=255)
                                       cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                                                                                          >  1:480233/1‑150 (MQ=255)
                                        gATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCGTCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCCAAGCGTGAACACGGAAATTAAGGTGAAGCCCAgc                                                                                                         >  1:2462655/1‑150 (MQ=255)
                                        gATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAgc                                                                                                         >  1:2256548/1‑150 (MQ=255)
                                             cTCCGTTATGATATGTTGGTCGGATACGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCAGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccac                                                                                                    >  1:733957/1‑150 (MQ=255)
                                              tCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccacc                                                                                                   <  2:2116644/150‑1 (MQ=255)
                                                cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGcc                                                                                                      >  1:2387314/1‑145 (MQ=255)
                                                     tGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCa                                                                                            <  1:1210034/150‑1 (MQ=255)
                                                     tGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCa                                                                                            >  2:1660217/1‑150 (MQ=255)
                                                     tGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCa                                                                                            <  1:2430456/150‑1 (MQ=255)
                                                       atatGTTGGTCGGATAATTCACTCGCGCCGCATCCGACATTGATTGCTTAATCGACTTCATTCACCTCACGACGAAACAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGc                                                                                          >  2:1002119/1‑150 (MQ=255)
                                                       atatGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGc                                                                                          <  2:1471014/150‑1 (MQ=255)
                                                         atcaatgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAc                                                                                        >  1:1130432/7‑150 (MQ=38)
                                                          tGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc                                                                                       >  1:1563043/1‑150 (MQ=255)
                                                           gTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCa                                                                                      >  1:2021446/1‑150 (MQ=255)
                                                           caatgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCa                                                                                      >  2:1265579/5‑150 (MQ=255)
                                                             tGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                    <  1:600493/150‑1 (MQ=255)
                                                             atgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                    >  1:1221517/3‑150 (MQ=255)
                                                             atgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                    >  2:1944824/3‑150 (MQ=255)
                                                             atgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCAACTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACCGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                    >  2:59630/3‑150 (MQ=255)
                                                              ggTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGCCCCAGCACCAGAt                                                                                   >  1:1711020/1‑150 (MQ=255)
                                                              ggTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGAt                                                                                   <  2:2322317/150‑1 (MQ=255)
                                                                tCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAATCGACTTCATTCACCTGACGACGCAGAATGGAAAGCGGTCCGCGGCCGCTAAGCATGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATaa                                                                                 >  2:1058265/1‑150 (MQ=255)
                                                                 cGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAg                                                                                <  2:80447/150‑1 (MQ=255)
                                                                   gATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACACAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAgcg                                                                              >  2:877581/1‑150 (MQ=255)
                                                                       aGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCt                                                                          >  1:774293/1‑150 (MQ=255)
                                                                        ggCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTg                                                                         >  1:2444983/1‑150 (MQ=255)
                                                                          cgcTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGa                                                                       >  1:1707068/1‑150 (MQ=255)
                                                                                 gccgcATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGAt                                                                <  2:2192249/150‑1 (MQ=255)
                                                                                  ccgcATCCTACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATg                                                               <  1:948189/150‑1 (MQ=255)
                                                                                   cgcATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                              <  1:699525/150‑1 (MQ=255)
                                                                                       tCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTc                                                          <  2:1517951/150‑1 (MQ=255)
                                                                                       tCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTc                                                          <  1:126765/150‑1 (MQ=255)
                                                                                        ccGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTca                                                         >  1:155762/1‑150 (MQ=255)
                                                                                        ccGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTca                                                         <  1:631348/150‑1 (MQ=255)
                                                                                        ccGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTca                                                         <  2:1130432/150‑1 (MQ=255)
                                                                                         cGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcat                                                        <  1:1343463/150‑1 (MQ=255)
                                                                                          gACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcata                                                       <  1:1660217/150‑1 (MQ=255)
                                                                                            cATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcataca                                                     <  1:120118/150‑1 (MQ=255)
                                                                                            cATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcataca                                                     >  2:2097810/1‑150 (MQ=255)
                                                                                              ttgattgCTTAAGCGACTTCATTCACCTGACGACGCAGAAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCATCACCAGATAAGCGACCTGGAAACCGATGCTTTcatacatt                                                   >  2:1328306/1‑149 (MQ=255)
                                                                                                   tgCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGcc                                              <  1:2268667/150‑1 (MQ=255)
                                                                                                    gCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcatacata                                                   <  2:1671567/144‑1 (MQ=255)
                                                                                                    gCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGccc                                             >  2:723272/1‑150 (MQ=255)
                                                                                                     cTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCg                                            >  1:1684241/1‑150 (MQ=255)
                                                                                                       tAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGGGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGcc                                          <  2:725760/150‑1 (MQ=255)
                                                                                                        aaGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCa                                         >  2:792760/1‑150 (MQ=255)
                                                                                                          gCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                              >  1:80270/1‑127 (MQ=255)
                                                                                                          gCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                              <  2:80270/127‑1 (MQ=255)
                                                                                                            gACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAATTTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATTCTTTAATAGATATTGCTTGCCAGTAc                                     <  2:1919255/150‑1 (MQ=255)
                                                                                                             aCTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACa                                    >  1:698499/1‑150 (MQ=255)
                                                                                                                  attcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACAt                               <  2:774293/150‑1 (MQ=255)
                                                                                                                   ttcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATa                              >  2:2075816/1‑150 (MQ=255)
                                                                                                                   ttcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATa                              >  2:2115446/1‑150 (MQ=255)
                                                                                                                    tcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATaa                             >  1:764327/1‑150 (MQ=255)
                                                                                                                     caCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGACCACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATaaa                            >  1:270898/1‑150 (MQ=255)
                                                                                                                      aCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATaaaa                           >  2:1392493/1‑150 (MQ=255)
                                                                                                                       ccTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATaaaaa                          <  1:808755/150‑1 (MQ=255)
                                                                                                                         tGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatc                        <  2:1855493/150‑1 (MQ=255)
                                                                                                                           acgacgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatcat                      <  1:1002119/150‑1 (MQ=255)
                                                                                                                                 cagcagGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAg                >  2:1058825/1‑150 (MQ=255)
                                                                                                                                    cagGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCg                    >  2:2370314/1‑143 (MQ=255)
                                                                                                                                     agGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGc            >  2:359231/1‑150 (MQ=255)
                                                                                                                                      gggAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCt           >  1:1425498/1‑150 (MQ=255)
                                                                                                                                      cgcttaGCGGCCCCGGCCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGACACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAACTCATCGCCAGTTGCt           >  1:359590/6‑150 (MQ=255)
                                                                                                                                        gAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTa         >  2:1990845/1‑150 (MQ=255)
                                                                                                                                          aaGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTaaa       >  1:355560/1‑150 (MQ=255)
                                                                                                                                             cGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAaagaa    >  2:1776474/1‑150 (MQ=255)
                                                                                                                                             cGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAaagaa    <  2:1563043/150‑1 (MQ=255)
                                                                                                                                              gggccggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAaagaag   >  1:1954320/1‑150 (MQ=255)
                                                                                                                                               ggccggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGc  <  2:557776/150‑1 (MQ=255)
                                                                                                                                               ggccggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGc  <  2:1666462/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGC  >  NC_000913/361824‑362116

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