Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,885,645 C→A 8.3% T224K (ACG→AAG)  mdtL → multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,885,6450CA 8.3% 1046.4 / 46.7 284T224K (ACG→AAG) mdtLmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base C (117/143);  new base A (9/15);  total (126/158)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 290 total aligned reads displayed.

TGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCAC  >  CP009273/3885498‑3885792
                                                                                                                                                   |                                                                                                                                                   
tGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGt                                                                                                                                                   >  2:1932164/1‑150 (MQ=255)
  tttATTTGAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCa                                                                                                                                                 <  2:153733/150‑1 (MQ=255)
     aTTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCg                                                                                                                                              >  1:3982097/1‑150 (MQ=255)
     aTTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCg                                                                                                                                              <  1:1840045/150‑1 (MQ=255)
        ttAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAACCCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTa                                                                                                                                           <  2:2859872/150‑1 (MQ=255)
        ttAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTa                                                                                                                                           >  1:188516/1‑150 (MQ=255)
         tAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTAt                                                                                                                                          <  1:2451142/150‑1 (MQ=255)
            aaGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGc                                                                                                                                       >  2:4570231/1‑150 (MQ=255)
            aaGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGc                                                                                                                                       >  1:3860731/1‑150 (MQ=255)
              gAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTg                                                                                                                                     >  2:3748021/1‑150 (MQ=255)
               aaaCGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGa                                                                                                                                    <  1:1689554/150‑1 (MQ=255)
               aaaCGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGa                                                                                                                                    <  1:1717212/150‑1 (MQ=255)
               aaaCGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGa                                                                                                                                    <  1:4939049/150‑1 (MQ=255)
                aaCGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGAt                                                                                                                                   <  1:1569092/150‑1 (MQ=255)
                aaCGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATg                                                                                                                                  >  1:1305790/1‑150 (MQ=255)
                 aCGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATg                                                                                                                                  >  1:833066/1‑150 (MQ=255)
                 aCGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATg                                                                                                                                  <  2:5060062/150‑1 (MQ=255)
                      cccAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAAt                                                                                                                             >  2:1790501/1‑150 (MQ=255)
                       ccAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATc                                                                                                                            >  2:1240524/1‑150 (MQ=255)
                        cAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCa                                                                                                                           <  1:283460/150‑1 (MQ=255)
                        cAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCa                                                                                                                           <  1:3845958/150‑1 (MQ=255)
                           cGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATgg                                                                                                                        <  2:4766082/150‑1 (MQ=255)
                           cGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATgg                                                                                                                        <  2:483853/150‑1 (MQ=255)
                           cGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACAAGTCACCGGTATTGCTGATGGAAATCATgg                                                                                                                        <  2:4244147/150‑1 (MQ=255)
                            ggCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATggg                                                                                                                       <  2:2521146/150‑1 (MQ=255)
                            ggCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATggg                                                                                                                       <  1:493135/150‑1 (MQ=255)
                              cccccGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTCTTGCTGATGGAAATCATggggg                                                                                                                     >  2:3145456/1‑149 (MQ=255)
                              cccccGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGt                                                                                                                     <  2:4366169/150‑1 (MQ=255)
                               ccccGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGtt                                                                                                                    <  2:924397/150‑1 (MQ=255)
                               ccccGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTGTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATAATGGGGtt                                                                                                                    <  2:3711044/150‑1 (MQ=255)
                                 ccGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTg                                                                                                                  <  2:2663643/150‑1 (MQ=255)
                                    cAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAgc                                                                                                               <  1:2898181/150‑1 (MQ=255)
                                     aGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGTgcg                                                                                                              >  1:597226/1‑150 (MQ=255)
                                       cTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGcgcg                                                                                                            >  2:2444317/1‑150 (MQ=255)
                                        ttCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCgg                                                                                                           <  2:1093484/150‑1 (MQ=255)
                                         tCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGt                                                                                                          <  2:1999067/150‑1 (MQ=255)
                                         tCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGt                                                                                                          >  2:3532446/1‑150 (MQ=255)
                                         tCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGt                                                                                                          >  2:3641340/1‑150 (MQ=255)
                                             aTAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAAt                                                                                                      >  2:3267803/1‑150 (MQ=255)
                                             aTAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGGATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAAt                                                                                                      >  1:766740/1‑150 (MQ=255)
                                             aTAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATAACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAAt                                                                                                      >  1:299634/1‑150 (MQ=255)
                                              tAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATa                                                                                                     >  2:2113813/1‑150 (MQ=255)
                                              tAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATa                                                                                                     <  1:4379441/150‑1 (MQ=255)
                                               aaaCCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATAc                                                                                                    <  1:3828023/150‑1 (MQ=255)
                                               aaaCCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACAAGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATAc                                                                                                    <  2:622877/150‑1 (MQ=255)
                                                aaCCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACg                                                                                                   <  1:4825589/150‑1 (MQ=255)
                                                aaCCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACg                                                                                                   >  2:2944141/1‑150 (MQ=255)
                                                      gagaAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGccacca                                                                                             >  1:1434148/1‑150 (MQ=255)
                                                        gaAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCAtt                                                                                           <  1:455442/150‑1 (MQ=255)
                                                         aaaaTAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACAAGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTa                                                                                          <  1:2712556/150‑1 (MQ=255)
                                                            aTAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATgg                                                                                       <  1:231073/150‑1 (MQ=255)
                                                            aTAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATgg                                                                                       >  1:1184511/1‑150 (MQ=255)
                                                            aTAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGGTTGAGCGCGGTGAATACGCCACCATTATgg                                                                                       >  2:18320/1‑150 (MQ=255)
                                                             tAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGgc                                                                                      <  2:2294703/150‑1 (MQ=255)
                                                             gcgCGAGTCGCTGCTTAACCGTGTTTTCCTCAGCCGTGGTGTGAGCACCCCCCTCAGCGTTTCGGTGCTCCTCCCTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCTTGGGGTTTGAGCGCGGTGAATCCGCCACCATTAGGgc                                                                                      <  2:1593787/148‑1 (MQ=255)
                                                                 gAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTg                                                                                  <  2:3819585/150‑1 (MQ=255)
                                                                    tCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGAcc                                                                               <  1:588552/150‑1 (MQ=255)
                                                                        tgctTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTg                                                                           <  1:4327405/150‑1 (MQ=255)
                                                                          ctTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGc                                                                         >  2:4719063/1‑150 (MQ=255)
                                                                           ttAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCg                                                                        <  2:4555112/150‑1 (MQ=255)
                                                                             aaCCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTc                                                                      >  2:1043019/1‑150 (MQ=255)
                                                                             aaCCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTc                                                                      <  2:4808867/150‑1 (MQ=255)
                                                                               ccGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACAAGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAg                                                                    <  1:1207835/150‑1 (MQ=255)
                                                                               ccGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACAAGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAg                                                                    <  1:4328124/150‑1 (MQ=255)
                                                                                  ttttttCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCAt                                                                 >  2:1989283/1‑150 (MQ=255)
                                                                                   tttttCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATg                                                                >  2:798021/1‑150 (MQ=255)
                                                                                    ttttCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGa                                                               >  2:1700203/1‑150 (MQ=255)
                                                                                         cTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGtt                                                          <  1:3058096/150‑1 (MQ=255)
                                                                                         cTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACAAGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGtt                                                          >  1:4299413/1‑150 (MQ=255)
                                                                                           cAGCCGCGTGGCGATCACCCCCCCCCGCGGTTCGGTGATCCTCACTATCGTCAACACGTCACCGGTATTGCAGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTc                                                        <  2:4483306/150‑1 (MQ=255)
                                                                                               cGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATtc                                                    <  1:4279436/150‑1 (MQ=255)
                                                                                                gtgtTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATtct                                                   >  2:5095835/1‑150 (MQ=255)
                                                                                                 tgttgtTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACAAGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATtctc                                                  <  2:3801443/150‑1 (MQ=255)
                                                                                                       tATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCAGGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGcc                                            >  1:4903080/1‑150 (MQ=255)
                                                                                                        aTCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCa                                           >  2:4761759/1‑150 (MQ=255)
                                                                                                            ccaccCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTg                                       <  2:1230605/150‑1 (MQ=255)
                                                                                                             caccCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTgc                                      <  1:4386147/150‑1 (MQ=255)
                                                                                                              accCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTgcg                                     >  1:4151796/1‑150 (MQ=255)
                                                                                                                ccTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCt                                   <  1:3359613/150‑1 (MQ=255)
                                                                                                                 cTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTg                                  <  1:2047289/150‑1 (MQ=255)
                                                                                                                 cTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTg                                  <  1:781981/150‑1 (MQ=255)
                                                                                                                    aGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGa                               <  2:1822350/150‑1 (MQ=255)
                                                                                                                    aGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGa                               >  2:4659765/1‑150 (MQ=255)
                                                                                                                     gCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGaa                              <  1:1039308/150‑1 (MQ=255)
                                                                                                                       gTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAAtt                            >  1:3700138/1‑150 (MQ=255)
                                                                                                                        tttCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAAttt                           <  1:4634253/150‑1 (MQ=255)
                                                                                                                        tttCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAAttt                           <  2:4281797/150‑1 (MQ=255)
                                                                                                                           cGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTa                        >  2:2103192/1‑150 (MQ=255)
                                                                                                                             gTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAg                      <  1:508811/150‑1 (MQ=255)
                                                                                                                             gTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAg                      <  1:2555072/150‑1 (MQ=255)
                                                                                                                              tGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGc                     <  2:1913001/150‑1 (MQ=255)
                                                                                                                               gATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGcc                    >  1:3165314/1‑150 (MQ=255)
                                                                                                                                   cTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCacgt                >  2:100967/1‑150 (MQ=255)
                                                                                                                                       cTTTCGTCAACACGTCACCGGTAGTGCTGATGGACATCATGGGGTTTGAGCGCGGTGAATACGCCACCCGGATGGCGCTGACCGCTGGCGTCAGCATGACCGTGTCATTCTCCACGCCCGTGGCGCTGGGCCTTTTGCAGCCccgtacgg            >  2:192895/1‑149 (MQ=255)
                                                                                                                                         ttCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTtg          <  1:3135412/150‑1 (MQ=255)
                                                                                                                                         ttCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTtg          <  1:4732148/150‑1 (MQ=255)
                                                                                                                                           cGTCAACACGTCACCTGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGTCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTtgat        <  1:1929984/150‑1 (MQ=255)
                                                                                                                                              cAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTgatgat     <  2:671695/150‑1 (MQ=255)
                                                                                                                                               aaCACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATc    <  2:1145553/150‑1 (MQ=255)
                                                                                                                                               aaCACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATc    <  2:1146010/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGTTTATTTTAAAAGAAACGCGCCCAGCGGCCCCCGCAGCTTCGGATAAACCACGAGAAAATAGCGAGTCGCTGCTTAACCGTTTTTTCCTCAGCCGTGTTGTTATCACCACCCTCAGCGTTTCGGTGATCCTCACTTTCGTCAACACGTCACCGGTATTGCTGATGGAAATCATGGGGTTTGAGCGCGGTGAATACGCCACCATTATGGCGCTGACCGCTGGCGTCAGCATGACCGTTTCATTCTCCACGCCATTTGCGCTGGGAATTTTTAAGCCACGTACGTTGATGATCAC  >  CP009273/3885498‑3885792

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