Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,971 G→C 6.2% P403P (CCC→CCG lacY ← lactose permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,9710GC 6.2% 548.9 / 19.6 162P403P (CCC→CCGlacYlactose permease
Reads supporting (aligned to +/- strand):  ref base G (78/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 = 3.29e-01

Only 100 of 174 total aligned reads displayed.

TGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGCAGAA  >  NC_000913/361824‑362120
                                                                                                                                                   |                                                                                                                                                     
 gCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccgggg                                                                                                                                                    <  1:918865/150‑1 (MQ=255)
 gCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAAGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggggc                                                                                                                                                   <  1:1516197/150‑1 (MQ=255)
   cTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCgggccggggcc                                                                                                                                                  <  2:2120496/150‑1 (MQ=255)
     gCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGc                                                                                                                                                >  2:1343463/1‑150 (MQ=255)
          tATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGc                                                                                                                                           >  2:261412/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)
            ttCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGt                                                                                                                                         <  1:190461/150‑1 (MQ=255)
            ttCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGt                                                                                                                                         >  1:183087/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)
                        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:2094478/150‑1 (MQ=255)
                            cATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAg                                                                                                                         >  2:68502/1‑150 (MQ=255)
                             aTGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAgg                                                                                                                        <  1:1998713/150‑1 (MQ=255)
                             aTGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAgg                                                                                                                        <  1:175409/150‑1 (MQ=255)
                             aTGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAgg                                                                                                                        <  2:1093162/150‑1 (MQ=255)
                                  cAATGCGATCACTCCGTTATTATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAg                                                                                                                   <  2:94118/150‑1 (MQ=255)
                                  cAATGCGATCACTCCGTTATGATATGTTGGTCGGATAATGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTTAACACGGAAATTAAGGTGAAg                                                                                                                   <  2:2114496/150‑1 (MQ=255)
                                  cAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAg                                                                                                                   >  2:1210034/1‑150 (MQ=255)
                                   aaTGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGc                                                                                                                  <  2:1516301/150‑1 (MQ=255)
                                       cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                                                                                              >  1:480233/1‑150 (MQ=255)
                                       cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                                                                                              <  1:2427717/150‑1 (MQ=255)
                                       cGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAg                                                                                                              >  2:1308933/1‑150 (MQ=255)
                                        gATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAgc                                                                                                             >  1:2256548/1‑150 (MQ=255)
                                         aTCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAgcg                                                                                                            <  2:2326730/150‑1 (MQ=255)
                                             cTCCGTTATGATATGTTGGTCGGATACGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCAGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccac                                                                                                        >  1:733957/1‑150 (MQ=255)
                                              tCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGccacc                                                                                                       <  2:2116644/150‑1 (MQ=255)
                                               ccGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAATGTGAAGCCCAGCGccacca                                                                                                      <  2:1967351/150‑1 (MQ=255)
                                                cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGcc                                                                                                          <  2:2387314/145‑1 (MQ=255)
                                                cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAg                                                                                                     >  2:2430456/1‑150 (MQ=255)
                                                cGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAg                                                                                                     <  2:365288/150‑1 (MQ=255)
                                                     tGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCa                                                                                                <  1:2430456/150‑1 (MQ=255)
                                                     tGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCa                                                                                                <  1:1210034/150‑1 (MQ=255)
                                                       atatGTTGGTCGGATAATTCACTCGCGCCGCATCCGACATTGATTGCTTAATCGACTTCATTCACCTCACGACGAAACAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGc                                                                                              >  2:1002119/1‑150 (MQ=255)
                                                       atatGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGc                                                                                              <  1:521839/150‑1 (MQ=255)
                                                         atcaatgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAc                                                                                            >  1:1130432/7‑150 (MQ=38)
                                                         atGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAc                                                                                            >  2:188627/1‑150 (MQ=255)
                                                          tGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc                                                                                           <  2:1725140/150‑1 (MQ=255)
                                                          tGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc                                                                                           <  2:145509/150‑1 (MQ=255)
                                                          tGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCAcc                                                                                           >  1:1563043/1‑150 (MQ=255)
                                                           gTTGGTCTGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCa                                                                                          <  1:2093998/150‑1 (MQ=255)
                                                           gTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCa                                                                                          <  2:723316/150‑1 (MQ=255)
                                                           caatgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCa                                                                                          >  2:1265579/5‑150 (MQ=255)
                                                            aatgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTCCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAg                                                                                         >  2:290401/4‑150 (MQ=255)
                                                             tGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                        <  1:600493/150‑1 (MQ=255)
                                                             atgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                        >  1:1221517/3‑150 (MQ=255)
                                                             atgTCGGATGCGGCGCGAGCGCCTTATCCGACATTGATTGCTTAAGCGACTTCATTCAACTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACCGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGa                                                                                        >  2:59630/3‑150 (MQ=255)
                                                                tCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAATCGACTTCATTCACCTGACGACGCAGAATGGAAAGCGGTCCGCGGCCGCTAAGCATGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATaa                                                                                     >  2:1058265/1‑150 (MQ=255)
                                                                 cGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAg                                                                                    <  2:80447/150‑1 (MQ=255)
                                                                 cGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAg                                                                                    <  1:1311484/150‑1 (MQ=255)
                                                                   gATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACACAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAgcg                                                                                  >  2:877581/1‑150 (MQ=255)
                                                                       aGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCt                                                                              >  1:774293/1‑150 (MQ=255)
                                                                          cgcTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGa                                                                           >  1:1707068/1‑150 (MQ=255)
                                                                              cgcgCCGCATCCGACGTTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAAcc                                                                       <  2:2237822/150‑1 (MQ=255)
                                                                              cgcgCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTGAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAAcc                                                                       <  1:258904/150‑1 (MQ=255)
                                                                                  ccgcATCCTACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATg                                                                   <  1:948189/150‑1 (MQ=255)
                                                                                   cgcATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                                  <  1:699525/150‑1 (MQ=255)
                                                                                      aTCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGTGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCttt                                                               <  1:1666267/150‑1 (MQ=255)
                                                                                       tCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTc                                                              <  2:1517951/150‑1 (MQ=255)
                                                                                       tCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTc                                                              <  1:126765/150‑1 (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)
                                                                                             attgattgCTTAAGCGCCTTCATTCACCTGACGCCGCAGCAGGGAAAGCGCGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGCACCAGCACCAGATAACCGCCCTGGAACCCGAGGCTTTcatacat                                                        >  2:1576088/1‑150 (MQ=255)
                                                                                                    gCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTcatacata                                                       >  1:1671567/1‑144 (MQ=255)
                                                                                                    gCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGccc                                                 >  2:1087556/1‑150 (MQ=255)
                                                                                                    gCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGccc                                                 >  2:723272/1‑150 (MQ=255)
                                                                                                     cTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCg                                                >  1:1684241/1‑150 (MQ=255)
                                                                                                       tAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGCTGCTTTCCTACATATTGCCCGcc                                              >  1:2451937/1‑150 (MQ=255)
                                                                                                        aaGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCa                                             >  2:792760/1‑150 (MQ=255)
                                                                                                          gCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                                  <  2:80270/127‑1 (MQ=255)
                                                                                                          gCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGc                                                                  >  1:80270/1‑127 (MQ=255)
                                                                                                             aCTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACa                                        >  1:698499/1‑150 (MQ=255)
                                                                                                                  attcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACAt                                   <  2:774293/150‑1 (MQ=255)
                                                                                                                   ttcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATa                                  >  2:2075816/1‑150 (MQ=255)
                                                                                                                    tcaCCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATaa                                 >  1:764327/1‑150 (MQ=255)
                                                                                                                         tGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatc                            <  2:1855493/150‑1 (MQ=255)
                                                                                                                           acgacgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatcat                          <  1:1002119/150‑1 (MQ=255)
                                                                                                                            cgacgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAatcatc                         >  2:1672607/1‑150 (MQ=255)
                                                                                                                              acgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACAGGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGGCAGTACAGACATAAAAATCATCGc                       <  2:1661150/150‑1 (MQ=255)
                                                                                                                               cgCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCATCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGcc                      <  2:120446/150‑1 (MQ=255)
                                                                                                                                 cagcagGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAg                    <  2:1922119/150‑1 (MQ=255)
                                                                                                                                    cagGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCg                        <  1:2370314/143‑1 (MQ=255)
                                                                                                                                      gggAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCt               >  1:1425498/1‑150 (MQ=255)
                                                                                                                                      gggAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCt               >  2:552442/1‑150 (MQ=255)
                                                                                                                                      cgcttaGCGGCCCCGGCCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGCAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAACTCATCGCCAGTTGCt               >  1:368518/6‑150 (MQ=255)
                                                                                                                                       ggAAAGCGGGCCGGGGCCGATAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCCCCAGACCCAGCACCAGATAAGCGCCCTGGACACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGATGCtt              >  1:1560826/1‑150 (MQ=255)
                                                                                                                                          aaGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTaaa           >  1:355560/1‑150 (MQ=255)
                                                                                                                                              gggccggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAaagaag       <  1:194215/150‑1 (MQ=255)
                                                                                                                                               ggccggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGc      <  2:1666462/150‑1 (MQ=255)
                                                                                                                                                   ggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGGCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGCAGaa  <  1:282590/150‑1 (MQ=255)
                                                                                                                                                   ggggccgCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGCAGaa  <  1:274028/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGGTCGGATAAGGCGCTCGCGCCGCATCCGACATTGATTGCTTAAGCGACTTCATTCACCTGACGACGCAGCAGGGAAAGCGGGCCGGGGCCGCTAAGCGTGAACACGGAAATTAAGGTGAAGCCCAGCGCCACCAGACCCAGCACCAGATAAGCGCCCTGGAAACCGATGCTTTCATACATATTGCCCGCCAGTACAGACATAAAAATCATCGCCAGTTGCTTAAAGAAGCAGAA  >  NC_000913/361824‑362120

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