Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 481,730 A→G Y892H (TAC→CAC)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913481,7300AG98.0% 362.1 / ‑1.2 98Y892H (TAC→CAC) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (52/44);  total (53/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TCTTCGCCGACAACCCAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAG  >  NC_000913/481583‑481877
                                                                                                                                                   |                                                                                                                                                   
tcttcGCCGACAACCCAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCa                                                                                                                                                   <  1:1670370/150‑1 (MQ=255)
tcttcGCCGACAACCCAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCa                                                                                                                                                   <  2:1302422/150‑1 (MQ=255)
 cttcGCCGACAACCCAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAgc                                                                                                                                                 >  2:685763/1‑151 (MQ=255)
       cGACAACCCAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCgccgcc                                                                                                                                           >  1:930440/1‑151 (MQ=255)
              ccAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGacac                                                                                                                                     >  2:329375/1‑150 (MQ=255)
               cAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGacac                                                                                                                                     >  1:212238/1‑149 (MQ=255)
               cAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGacac                                                                                                                                     >  2:897978/1‑149 (MQ=255)
                aaTGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAgg                                                                                                                                  <  2:1143651/151‑1 (MQ=255)
                  tGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGa                                                                                                                                 >  1:670406/1‑150 (MQ=255)
                    gTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAaca                                                                                                                              <  1:689735/151‑1 (MQ=255)
                    gTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAaca                                                                                                                              <  1:1523298/151‑1 (MQ=255)
                    gTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAac                                                                                                                               >  1:1690431/1‑150 (MQ=255)
                     ttGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAacac                                                                                                                             >  1:1249018/1‑151 (MQ=255)
                     ttGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAacac                                                                                                                             >  2:1357901/1‑151 (MQ=255)
                       gtgAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACg                                                                                                                            >  1:1250299/1‑150 (MQ=255)
                       gtgAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGa                                                                                                                           <  2:1643270/151‑1 (MQ=255)
                         gAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGAc                                                                                                                          <  2:135421/150‑1 (MQ=255)
                          aGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACaa                                                                                                                        <  2:809036/151‑1 (MQ=255)
                           gCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAAt                                                                                                                       <  1:248244/151‑1 (MQ=255)
                            cAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATc                                                                                                                      >  2:726580/1‑151 (MQ=255)
                             aGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATc                                                                                                                      >  1:250380/1‑150 (MQ=255)
                             aGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCa                                                                                                                     <  2:1575310/151‑1 (MQ=255)
                             aGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCa                                                                                                                     <  1:474356/151‑1 (MQ=255)
                                   aCCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGaaa                                                                                                               <  2:1671698/151‑1 (MQ=255)
                                   aCCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGaa                                                                                                                <  1:1409997/150‑1 (MQ=255)
                                    ccTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAAt                                                                                                              <  2:883049/151‑1 (MQ=255)
                                      tGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATc                                                                                                             >  2:1504256/1‑150 (MQ=255)
                                       ggAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATc                                                                                                             <  1:1194744/149‑1 (MQ=255)
                                        gAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATcgcg                                                                                                          >  1:1022357/1‑151 (MQ=255)
                                          aGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGt                                                                                                         >  1:1516083/1‑150 (MQ=255)
                                           gTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTAc                                                                                                       <  1:723151/151‑1 (MQ=255)
                                            tAAACGTCATTGGTCAGGCCACGGAAGGCGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTAc                                                                                                       >  1:153684/1‑150 (MQ=255)
                                             aaaCGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTAc                                                                                                       >  1:1455897/1‑149 (MQ=255)
                                              aaCGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGt                                                                                                    <  1:379577/151‑1 (MQ=255)
                                              aaCGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGt                                                                                                    <  2:1298780/151‑1 (MQ=255)
                                                   cATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAgg                                                                                               >  1:491543/1‑151 (MQ=255)
                                                   cATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAgg                                                                                               >  2:589045/1‑151 (MQ=255)
                                                   cATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAgg                                                                                               <  2:302798/151‑1 (MQ=255)
                                                    aTTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGt                                                                                              <  1:1714694/151‑1 (MQ=255)
                                                         tCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTg                                                                                         >  2:986976/1‑151 (MQ=255)
                                                          cAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTgg                                                                                        >  1:1619232/1‑151 (MQ=255)
                                                           aGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGt                                                                                       <  1:1430741/151‑1 (MQ=255)
                                                            ggCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGtt                                                                                      <  2:1116070/151‑1 (MQ=255)
                                                              ccACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGc                                                                                    <  1:336918/151‑1 (MQ=255)
                                                                 cGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCgg                                                                                 >  2:1013973/1‑151 (MQ=255)
                                                                 cGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCgg                                                                                 <  1:1638501/151‑1 (MQ=255)
                                                                  ggAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCgg                                                                                 >  1:1707072/1‑150 (MQ=255)
                                                                   gAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCgg                                                                                 >  1:1045202/1‑149 (MQ=255)
                                                                   gAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGgag                                                                               <  2:411154/151‑1 (MQ=255)
                                                                     aGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGgagag                                                                             <  2:376266/151‑1 (MQ=255)
                                                                      ggTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGgag                                                                               >  2:270873/1‑148 (MQ=255)
                                                                      ggTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGgag                                                                               >  1:1135636/1‑148 (MQ=255)
                                                                       gTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGAc                                                                           <  1:827728/151‑1 (MQ=255)
                                                                       gTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGAc                                                                           >  2:1357972/1‑151 (MQ=255)
                                                                         ggCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACg                                                                          <  1:1478867/150‑1 (MQ=255)
                                                                           cagccagcAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTc                                                                       >  2:1117502/1‑151 (MQ=255)
                                                                           cagccagcAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTc                                                                       >  1:1470081/1‑151 (MQ=255)
                                                                            agccagcAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTcc                                                                      >  2:1444836/1‑151 (MQ=255)
                                                                                agcAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGAt                                                                  >  2:1016231/1‑151 (MQ=255)
                                                                                agcAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGAt                                                                  <  1:811501/151‑1 (MQ=255)
                                                                                  cAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGAt                                                                  >  2:1268492/1‑149 (MQ=255)
                                                                                   aaCGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAgg                                                               >  1:206838/1‑151 (MQ=255)
                                                                                   aaCGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAgg                                                               >  1:674009/1‑151 (MQ=255)
                                                                                     cGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGAc                                                             <  2:1012374/151‑1 (MQ=255)
                                                                                     cGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGAc                                                             >  2:527116/1‑151 (MQ=255)
                                                                                        aCCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATc                                                          >  1:945141/1‑151 (MQ=255)
                                                                                        aCCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATc                                                          >  2:1459969/1‑151 (MQ=255)
                                                                                        aCCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATc                                                          >  2:6799/1‑151 (MQ=255)
                                                                                            aTAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGAGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCg                                                      >  1:975550/1‑151 (MQ=255)
                                                                                              aaCCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTc                                                    <  1:530518/151‑1 (MQ=255)
                                                                                                cccccAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCa                                                  <  2:795717/151‑1 (MQ=255)
                                                                                                   ccAGCGGAACGACCAGCATAACGGAGCACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTc                                               >  1:768663/1‑151 (MQ=255)
                                                                                                   ccAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTc                                               <  2:988630/151‑1 (MQ=255)
                                                                                                   ccAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTc                                               >  1:1027367/1‑151 (MQ=255)
                                                                                                    cAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCa                                              <  2:1652275/151‑1 (MQ=255)
                                                                                                     aGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCa                                              >  2:1329593/1‑150 (MQ=255)
                                                                                                        ggAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGc                                          >  1:1008115/1‑151 (MQ=255)
                                                                                                          aaCGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCa                                        <  2:250403/151‑1 (MQ=255)
                                                                                                           aCGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCa                                        >  2:97842/1‑150 (MQ=255)
                                                                                                                   cATAACGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTa                               <  1:160503/151‑1 (MQ=255)
                                                                                                                       aCGGAGAACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGc                            >  2:691253/1‑150 (MQ=255)
                                                                                                                           agaACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGttt                       >  2:1149450/1‑151 (MQ=255)
                                                                                                                           agaACGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGttt                       >  1:1273996/1‑151 (MQ=255)
                                                                                                                             aaCGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGc                     <  2:1690546/151‑1 (MQ=255)
                                                                                                                              aCGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGc                     >  1:1511433/1‑150 (MQ=255)
                                                                                                                               cGGAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTc                   <  1:685705/151‑1 (MQ=255)
                                                                                                                                 gAATCGACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGc                 >  1:1378030/1‑151 (MQ=255)
                                                                                                                                   aTCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGcc                >  2:1535348/1‑150 (MQ=255)
                                                                                                                                      gACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgtt            <  1:949146/151‑1 (MQ=255)
                                                                                                                                      gACCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgtt            <  2:1210572/151‑1 (MQ=255)
                                                                                                                                       aCCAGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAgttg           <  1:1187311/151‑1 (MQ=255)
                                                                                                                                          aGCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTc        >  2:90163/1‑151 (MQ=255)
                                                                                                                                           gCTCTCGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTcc       >  2:1102001/1‑151 (MQ=255)
                                                                                                                                             tctcGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCAt     >  1:190503/1‑151 (MQ=255)
                                                                                                                                                cGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATc    <  2:1249115/149‑1 (MQ=255)
                                                                                                                                                cGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCa   <  1:277699/150‑1 (MQ=255)
                                                                                                                                                cGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCa   <  1:492744/150‑1 (MQ=255)
                                                                                                                                                cGTGCAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAg  <  2:894253/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TCTTCGCCGACAACCCAATGGTTGTGAGCAGGCCTACCTGGAAGTAAACGTCATTGGTCAGGCCACGGAAGGTGGCAGCCAGCAACGCACCGATAACCCCCAGCGGAACGACCAGCATAACGGAGAACGGAATCGACCAGCTCTCGTACAGCGCCGCCAGACACAGGAACACGACAATCAACGAAATCGCGTACAGTGAAGGTGCCTGGTTGCCGGAGAGACGTTCCTGATAGGACATCCCCGTCCAGTCATAGCCAACACCGGTAGGCAGTTTGCTCGCCAGTTGTTCCATCAG  >  NC_000913/481583‑481877

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