Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,154,021 G→A 6.3% D719N (GAC→AAC)  mdtC → multidrug efflux system, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,154,0210GA 6.3% 827.9 / 29.5 254D719N (GAC→AAC) mdtCmultidrug efflux system, subunit C
Reads supporting (aligned to +/- strand):  ref base G (117/121);  new base A (4/12);  total (121/133)
Fisher's exact test for biased strand distribution p-value = 7.28e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 258 total aligned reads displayed.

ACCTGGCGGCACTGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTT  >  CP009273/2153872‑2154170
                                                                                                                                                     |                                                                                                                                                     
acCTGGCGGCACTGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCg                                                                                                                                                       <  2:3616950/150‑1 (MQ=255)
   tGGCGGCACTGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACg                                                                                                                                                    <  2:1403084/150‑1 (MQ=255)
   tGGCGGCACTGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACg                                                                                                                                                    <  2:5081561/150‑1 (MQ=255)
       ggcACTGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACa                                                                                                                                                >  1:1935639/1‑150 (MQ=255)
        gcACTGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACaa                                                                                                                                               >  1:1534565/1‑150 (MQ=255)
           cTGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAgcc                                                                                                                                            <  2:2023799/150‑1 (MQ=255)
            tGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCAACGTACAAgccg                                                                                                                                           <  1:4222129/150‑1 (MQ=255)
             gCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAgccgc                                                                                                                                          <  2:4234981/150‑1 (MQ=255)
               gagaATGGGAGCCGAAAATCCGCTATCAACTGGCGTCGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATGTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGCCGTACAAGCCGCCa                                                                                                                                        <  2:2383518/150‑1 (MQ=255)
               gagaATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCa                                                                                                                                        >  2:1865608/1‑150 (MQ=255)
                agaATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCaa                                                                                                                                       >  1:147710/1‑150 (MQ=255)
                agaATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCaa                                                                                                                                       <  1:2600118/150‑1 (MQ=255)
                   aTGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAg                                                                                                                                    <  1:2065834/150‑1 (MQ=255)
                    tGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGt                                                                                                                                   >  2:363311/1‑150 (MQ=255)
                    tGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGt                                                                                                                                   >  2:594702/1‑150 (MQ=255)
                     gggAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTc                                                                                                                                  >  1:4092163/1‑150 (MQ=255)
                     gggAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTc                                                                                                                                  >  2:842996/1‑150 (MQ=255)
                      ggAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCt                                                                                                                                 <  2:750314/150‑1 (MQ=255)
                      ggAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCt                                                                                                                                 <  1:4607782/150‑1 (MQ=255)
                      ggAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATAATCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCt                                                                                                                                 >  1:2702380/1‑150 (MQ=255)
                          ccGAACATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTa                                                                                                                             <  2:4480578/150‑1 (MQ=255)
                            gAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTaaa                                                                                                                           <  1:3295889/150‑1 (MQ=255)
                            gAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTaaa                                                                                                                           >  2:4960940/1‑150 (MQ=255)
                             aaaaTCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAAt                                                                                                                          <  1:4868739/150‑1 (MQ=255)
                             aaaaTCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAAt                                                                                                                          >  1:2781324/1‑150 (MQ=255)
                             aaaaTCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAAt                                                                                                                          <  1:4150294/150‑1 (MQ=255)
                              aaaTCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATa                                                                                                                         <  1:5069957/150‑1 (MQ=255)
                               aaTCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATaa                                                                                                                        >  1:1420823/1‑150 (MQ=255)
                                aTCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAAc                                                                                                                       >  1:1002221/1‑150 (MQ=255)
                                aTCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAAc                                                                                                                       <  2:117757/150‑1 (MQ=255)
                                 tCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACg                                                                                                                      <  2:1389553/150‑1 (MQ=255)
                                              ggCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCg                                                                                                         >  2:3633762/1‑150 (MQ=255)
                                                cGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGc                                                                                                       <  1:2674349/150‑1 (MQ=255)
                                                    gTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAAt                                                                                                   <  1:339141/150‑1 (MQ=255)
                                                        ccGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCg                                                                                               <  2:2200779/150‑1 (MQ=255)
                                                          ggAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAAACTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGAc                                                                                             <  2:869947/150‑1 (MQ=255)
                                                            aaCTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCAACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCa                                                                                           <  1:1342985/150‑1 (MQ=255)
                                                             aCTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCAt                                                                                          <  1:530178/150‑1 (MQ=255)
                                                               tGGCGGACATGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCAttt                                                                                        <  2:4480113/150‑1 (MQ=255)
                                                                ggCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa                                                                                       <  2:2573307/150‑1 (MQ=255)
                                                                ggCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTa                                                                                       >  2:1809122/1‑150 (MQ=255)
                                                                 gCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTAc                                                                                      <  1:3335531/150‑1 (MQ=255)
                                                                      cGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGcc                                                                                 >  2:2810034/1‑150 (MQ=255)
                                                                      cGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGcc                                                                                 >  1:3570290/1‑150 (MQ=255)
                                                                         gAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGAt                                                                              >  2:1628181/1‑150 (MQ=255)
                                                                           aCTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGa                                                                            >  1:2730687/1‑150 (MQ=255)
                                                                           aCTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGa                                                                            >  1:2284811/1‑150 (MQ=255)
                                                                                 aTCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGt                                                                      >  1:3052745/1‑150 (MQ=255)
                                                                                  tCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGta                                                                     <  2:3690177/150‑1 (MQ=255)
                                                                                    agcagGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGtata                                                                   <  2:959315/150‑1 (MQ=255)
                                                                                    agcagGATAACGGCGCGGAGATGAATCGGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGtata                                                                   <  2:3846078/150‑1 (MQ=255)
                                                                                       agGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAg                                                                <  2:299184/150‑1 (MQ=255)
                                                                                       agGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAg                                                                >  2:2712995/1‑150 (MQ=255)
                                                                                       agGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAg                                                                >  2:3918463/1‑150 (MQ=255)
                                                                                       agGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAg                                                                >  2:4894492/1‑150 (MQ=255)
                                                                                          aTAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgtgg                                                             >  2:2042890/1‑150 (MQ=255)
                                                                                           tAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgtggt                                                            <  2:147710/150‑1 (MQ=255)
                                                                                           tAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgtggt                                                            >  1:3113418/1‑150 (MQ=255)
                                                                                           tAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgtggt                                                            >  1:505408/1‑150 (MQ=255)
                                                                                           tAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCAACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgtggt                                                            <  2:2123259/150‑1 (MQ=255)
                                                                                            aaCGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgtggtg                                                           >  1:1204249/1‑150 (MQ=255)
                                                                                            aaCGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgtggtg                                                           >  2:3928767/1‑150 (MQ=255)
                                                                                            aaCGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAgtggtg                                                           >  1:2583708/1‑150 (MQ=255)
                                                                                             aCGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGa                                                          >  2:1870732/1‑150 (MQ=255)
                                                                                                  gcgGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGaa                                                     <  1:601971/150‑1 (MQ=255)
                                                                                                    ggAGATGACTCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGt                                                   <  2:1792588/150‑1 (MQ=255)
                                                                                                    ggAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCAACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGt                                                   <  1:3867655/150‑1 (MQ=255)
                                                                                                    ggAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCAACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGt                                                   <  2:1689961/150‑1 (MQ=255)
                                                                                                     gagaTGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTg                                                  <  1:371856/150‑1 (MQ=255)
                                                                                                     gagaTGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCAACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTg                                                  >  1:3617701/1‑150 (MQ=255)
                                                                                                      agaTGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCTAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTgg                                                 <  2:899200/150‑1 (MQ=255)
                                                                                                       gaTGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGa                                                >  1:219148/1‑150 (MQ=255)
                                                                                                        aTGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGAt                                               <  2:2419056/150‑1 (MQ=255)
                                                                                                        aTGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGAt                                               <  1:3970913/150‑1 (MQ=255)
                                                                                                         tGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATc                                              <  2:4107074/150‑1 (MQ=255)
                                                                                                         tGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATc                                              <  1:5135873/150‑1 (MQ=255)
                                                                                                         tGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCAACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATc                                              >  2:3709387/1‑150 (MQ=255)
                                                                                                          gAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATcc                                             <  2:4843019/150‑1 (MQ=255)
                                                                                                          gAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATcc                                             >  1:696666/1‑150 (MQ=255)
                                                                                                           aaTCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCcg                                            >  2:157540/1‑150 (MQ=255)
                                                                                                                ggTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCta                                       >  2:1871858/1‑150 (MQ=255)
                                                                                                                 gTTTACGACCGCGACACCATGGCACGGCTGGGAATCAACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCtat                                      <  1:4036361/150‑1 (MQ=255)
                                                                                                                    tACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATAcc                                   <  1:1089803/150‑1 (MQ=255)
                                                                                                                           gcgACACCATGGCTCGGCTGGGAATCGACTTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGAGCTATACCCAGGACa                            <  1:1763255/150‑1 (MQ=255)
                                                                                                                            cgACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACAt                           >  2:3411337/1‑150 (MQ=255)
                                                                                                                            cgACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACAt                           <  1:4521210/150‑1 (MQ=255)
                                                                                                                              acacCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCa                         <  1:2536085/150‑1 (MQ=255)
                                                                                                                               cacCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAg                        >  2:2596984/1‑150 (MQ=255)
                                                                                                                                acCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGt                       >  1:1670880/1‑150 (MQ=255)
                                                                                                                                 ccATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTg                      <  1:1899401/150‑1 (MQ=255)
                                                                                                                                  cATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGACGTGGATCCGCGCTATACCCAGGACATCAGTgc                     >  2:3596243/1‑150 (MQ=255)
                                                                                                                                     ggCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCt                  <  2:4129910/150‑1 (MQ=255)
                                                                                                                                      gCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTg                 <  2:3174386/150‑1 (MQ=255)
                                                                                                                                        aCGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGa               >  1:4889956/1‑150 (MQ=255)
                                                                                                                                         cGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGaa              >  1:4410268/1‑150 (MQ=255)
                                                                                                                                            cTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGaaaaa           >  1:5022962/1‑150 (MQ=255)
                                                                                                                                               ggAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATg        <  2:1023047/150‑1 (MQ=255)
                                                                                                                                                  aTCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTc     >  2:2228266/1‑150 (MQ=255)
                                                                                                                                                     gACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAACAAATGTTCGtt  >  2:2469165/1‑150 (MQ=255)
                                                                                                                                                     gACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGtt  <  1:877559/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCTGGCGGCACTGCGAGAATGGGAGCCGAAAATCCGCAAAAAACTGGCGACGTTGCCGGAACTGGCGGACGTGAACTCCGATCAGCAGGATAACGGCGCGGAGATGAATCTGGTTTACGACCGCGACACCATGGCACGGCTGGGAATCGACGTACAAGCCGCCAACAGTCTGTTAAATAACGCCTTCGGTCAGCGGCAAATCTCGACCATTTACCAGCCGATGAACCAGTATAAAGTGGTGATGGAAGTGGATCCGCGCTATACCCAGGACATCAGTGCGCTGGAAAAAATGTTCGTT  >  CP009273/2153872‑2154170

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: