Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,885,881 G→T 5.6% A303S (GCC→TCC)  mdtL → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,885,8810GT 5.6% 1269.5 / 29.0 329A303S (GCC→TCC) mdtLmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (184/126);  new base T (5/14);  total (189/140)
Fisher's exact test for biased strand distribution p-value = 7.16e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

CCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCGGTGGTTGGTATCGGC  >  CP009273/3885732‑3886030
                                                                                                                                                     |                                                                                                                                                     
      cATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGtt                                                                                                                                                 <  1:340905/150‑1 (MQ=255)
      cATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGtt                                                                                                                                                 <  2:2636533/150‑1 (MQ=255)
        ttCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTtc                                                                                                                                               >  1:1521731/1‑150 (MQ=255)
         tctcCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTtct                                                                                                                                              >  2:2110298/1‑150 (MQ=255)
          ctcCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTtctc                                                                                                                                             >  2:996543/1‑150 (MQ=255)
          ctcCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTtctc                                                                                                                                             <  1:3867948/150‑1 (MQ=255)
            ccACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCgg                                                                                                                                           >  2:3432731/1‑150 (MQ=255)
            ccACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCgg                                                                                                                                           >  2:2664880/1‑150 (MQ=255)
                gCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAgg                                                                                                                                       >  2:4428323/1‑150 (MQ=255)
                 ccATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGt                                                                                                                                      >  2:4149899/1‑150 (MQ=255)
                 ccATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTACACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGt                                                                                                                                      >  2:3749712/1‑150 (MQ=255)
                   aTTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGttt                                                                                                                                    >  1:4567527/1‑150 (MQ=255)
                     ttGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTg                                                                                                                                  <  2:144918/150‑1 (MQ=255)
                        cgcTGGGAATTTTTAAGCCACGTACGTTGATGATCACCACGCAGGTGTTATTCCTGGCGGCGGGTATCACTCTAGCCGTTTCACCTTCCCAAGCGGTATCTCTGTTAGGAATCACGCAGAAATGCGCCGGGGACACGGTAGGTTTGGgtg                                                                                                                               >  2:2210874/1‑150 (MQ=255)
                          cTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGtgtg                                                                                                                             >  1:2763999/1‑150 (MQ=255)
                             ggAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCg                                                                                                                          <  2:3378283/150‑1 (MQ=255)
                             ggAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCg                                                                                                                          >  1:1157984/1‑150 (MQ=255)
                             ggAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCg                                                                                                                          <  1:2608762/150‑1 (MQ=255)
                             ggAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCg                                                                                                                          <  1:2051620/150‑1 (MQ=255)
                                   tttAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGt                                                                                                                    <  2:1765420/150‑1 (MQ=255)
                                    ttAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTc                                                                                                                   <  1:2926515/150‑1 (MQ=255)
                                       aGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAgg                                                                                                                <  1:2209560/150‑1 (MQ=255)
                                       aGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAgg                                                                                                                >  1:888234/1‑150 (MQ=255)
                                         ccACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCg                                                                                                              >  2:1296647/1‑150 (MQ=255)
                                          cACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGt                                                                                                             >  2:546596/1‑150 (MQ=255)
                                          cACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGt                                                                                                             >  1:4059695/1‑150 (MQ=255)
                                                    gatgatCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGtt                                                                                                   <  2:1686049/150‑1 (MQ=255)
                                                    gatgatCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGtt                                                                                                   >  1:739318/1‑150 (MQ=255)
                                                    gatgatCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGtt                                                                                                   >  1:3843888/1‑150 (MQ=255)
                                                     atgatCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGtgt                                                                                                  >  1:4398988/1‑148 (MQ=255)
                                                     aggatCCCCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGttt                                                                                                  <  2:3659098/148‑1 (MQ=255)
                                                       gatCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTc                                                                                                <  1:253621/150‑1 (MQ=255)
                                                        atCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCa                                                                                               >  1:3030885/1‑150 (MQ=255)
                                                        atCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCa                                                                                               >  2:3337968/1‑150 (MQ=255)
                                                         tCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCCCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCAt                                                                                              >  1:1814233/1‑150 (MQ=255)
                                                          cACCTCGCAGGTGTTATTCCTGGCTGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTATGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCAtt                                                                                             <  1:3930626/150‑1 (MQ=255)
                                                           aCCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTa                                                                                            <  1:4049521/150‑1 (MQ=255)
                                                            ccTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGTTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAc                                                                                           >  1:3014603/1‑150 (MQ=255)
                                                                 cAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGAGGTTTCTCTGTTTGGTATCACTCTGATTTGCTACGGTTTATCGGGAGGTTTTGGTGTGGCGATGAGTAAGGCGTTAGGGCCGTTTTAATTAcgagcg                                                                                      >  2:3597423/1‑150 (MQ=255)
                                                                  aGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCgg                                                                                     >  2:4623401/1‑150 (MQ=255)
                                                                  aGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCgg                                                                                     >  1:3250059/1‑150 (MQ=255)
                                                                  aGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCgg                                                                                     >  1:1492080/1‑150 (MQ=255)
                                                                   ggTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCggg                                                                                    >  2:4317334/1‑150 (MQ=255)
                                                                     tgtTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCg                                                                                  <  1:1577993/150‑1 (MQ=255)
                                                                        tATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAg                                                                               >  1:2691173/1‑150 (MQ=255)
                                                                        tATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAg                                                                               >  2:4027160/1‑150 (MQ=255)
                                                                         aTTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGc                                                                              >  2:3487642/1‑150 (MQ=255)
                                                                         aTTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGc                                                                              >  2:1045117/1‑150 (MQ=255)
                                                                           tCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCa                                                                            <  2:3927487/150‑1 (MQ=255)
                                                                              tGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCTCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCt                                                                         <  2:4588101/150‑1 (MQ=255)
                                                                              tGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCAATACGCGCGGGCGTATCCAGCt                                                                         >  1:369579/1‑150 (MQ=255)
                                                                                gcggcgGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCg                                                                       <  1:4392216/150‑1 (MQ=255)
                                                                                gcggcgGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCg                                                                       <  1:1003700/150‑1 (MQ=255)
                                                                                 cggcggGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGa                                                                      >  2:1534269/1‑150 (MQ=255)
                                                                                  ggcggGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGAc                                                                     >  2:4549502/1‑150 (MQ=255)
                                                                                    cggGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCTCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCt                                                                   <  2:3771202/150‑1 (MQ=255)
                                                                                    cggGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCTCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCt                                                                   <  2:735387/150‑1 (MQ=255)
                                                                                    cggGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCt                                                                   >  1:4588980/1‑150 (MQ=255)
                                                                                    cggGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCt                                                                   <  2:639077/150‑1 (MQ=255)
                                                                                      gggATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTa                                                                 >  1:636084/1‑150 (MQ=255)
                                                                                      gggATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTa                                                                 >  1:2975079/1‑150 (MQ=255)
                                                                                      gggATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTGTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTTTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAACCGCGGGCGTAGCCAGCCCGACCtga                                                                 >  2:3110866/1‑148 (MQ=255)
                                                                                      gggATCACTCTTGCCGTTACACCTTCCCATGCGGTTTCTCAGTTTGGTAACACGCTGAATTGAGCAGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCATGCGTAAGGGCCTTGTTCATTACGAGCGGGCGAAGCCAGCACGAcccaa                                                                 >  1:1418829/1‑147 (MQ=255)
                                                                                       ggATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAg                                                                <  1:3724957/150‑1 (MQ=255)
                                                                                        gATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAgg                                                               <  2:1981778/150‑1 (MQ=255)
                                                                                         aTCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGt                                                              >  2:1582393/1‑150 (MQ=255)
                                                                                          tCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTa                                                             <  1:16373/150‑1 (MQ=255)
                                                                                            aCTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTATGGCCGTTTTAATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTAtt                                                           >  1:1320608/1‑150 (MQ=255)
                                                                                             ctctTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCAATACGCGCGGGCGTAGCCACCACGACCCTAggggatg                                                          >  1:1991634/1‑145 (MQ=255)
                                                                                                ttGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCCGCTCGACCGTAGGTATTgcgc                                                       >  1:2374859/1‑150 (MQ=255)
                                                                                                 tGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCTCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCa                                                      <  1:3985659/150‑1 (MQ=255)
                                                                                                 tGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCa                                                      <  1:3575174/150‑1 (MQ=255)
                                                                                                    cGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGGATTGCGCAGGt                                                   >  2:2249590/1‑150 (MQ=255)
                                                                                                      tttCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGttt                                                 >  1:1827193/1‑150 (MQ=255)
                                                                                                          aCCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCgg                                             <  1:2958925/150‑1 (MQ=255)
                                                                                                           ccTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGt                                            >  2:3113648/1‑150 (MQ=255)
                                                                                                            cTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGtt                                           <  2:314585/150‑1 (MQ=255)
                                                                                                            cTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGtt                                           <  2:87557/150‑1 (MQ=255)
                                                                                                             ttCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTATCCACCTCGACCTTAGGTATTTCGCAGGTTTGCGGTTc                                          >  1:1741018/1‑150 (MQ=255)
                                                                                                             ttCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTc                                          >  2:3271083/1‑150 (MQ=255)
                                                                                                              tCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCGTAGGTATTGCGCAGGTTTGCGGTTCg                                         >  1:1655522/1‑150 (MQ=255)
                                                                                                                   tGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACt                                    >  1:171425/1‑150 (MQ=255)
                                                                                                                   tGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCGTAGGTATTGCGCAGGTTTGCGGTTCGTCACt                                    >  2:2820885/1‑150 (MQ=255)
                                                                                                                   tGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTAAGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTTAGCAGGTTTGCGGTTAGTCAca                                    >  1:3603734/1‑149 (MQ=255)
                                                                                                                     cGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACtgt                                  <  2:3323909/150‑1 (MQ=255)
                                                                                                                     cGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGGTAGGGCCGTTTTCATTAAGCGCGGGCGTAGCCAGCTCGAACTTAGGTATTGCGCAGGTTTGCGGTTCGTCACtgt                                  >  2:2136680/1‑150 (MQ=255)
                                                                                                                      ggTTTCTCTGTTTGGTATCACGCTGATTTGCTCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACtgtg                                 <  2:3851245/150‑1 (MQ=255)
                                                                                                                      ggTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACtgtg                                 <  1:514778/150‑1 (MQ=255)
                                                                                                                        tttCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGa                               >  2:1095323/1‑150 (MQ=255)
                                                                                                                         ttCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGAt                              <  2:3036348/150‑1 (MQ=255)
                                                                                                                           ctctGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGAttt                            <  1:426840/150‑1 (MQ=255)
                                                                                                                               gTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTtggct                        <  2:2977707/150‑1 (MQ=255)
                                                                                                                                  tGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTtggctggc                     >  1:4445310/1‑150 (MQ=255)
                                                                                                                                      aTCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCg                 >  2:4454955/1‑150 (MQ=255)
                                                                                                                                       tCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCgg                >  1:2279121/1‑150 (MQ=255)
                                                                                                                                         aCGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCggtg              <  2:2523905/150‑1 (MQ=255)
                                                                                                                                            cTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCGGtggtt           <  1:4293139/150‑1 (MQ=255)
                                                                                                                                              gATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCGGtggttgg         <  1:4160165/150‑1 (MQ=255)
                                                                                                                                                tttGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCGGTGGTTGGTa       <  1:2110298/150‑1 (MQ=255)
                                                                                                                                                  tGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCGGTGGTTGGTATc     >  2:2098494/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCACCTCGCAGGTGTTATTCCTGGCGGCGGGGATCACTCTTGCCGTTTCACCTTCCCATGCGGTTTCTCTGTTTGGTATCACGCTGATTTGCGCCGGTTTCTCGGTAGGTTTTGGTGTGGCGATGAGTCAGGCGTTAGGGCCGTTTTCATTACGCGCGGGCGTAGCCAGCTCGACCTTAGGTATTGCGCAGGTTTGCGGTTCGTCACTGTGGATTTGGCTGGCAGCGGTGGTTGGTATCGGC  >  CP009273/3885732‑3886030

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