Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,150,896 C→A 5.0% A718E (GCG→GAG)  mdtB → multidrug efflux system, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,150,8960CA 5.0% 687.1 / 9.7 177A718E (GCG→GAG) mdtBmultidrug efflux system, subunit B
Reads supporting (aligned to +/- strand):  ref base C (86/82);  new base A (4/5);  total (90/87)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 177 total aligned reads displayed.

GGTCAGCCGCACCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACA  >  CP009273/2150750‑2151042
                                                                                                                                                  |                                                                                                                                                  
   cAGCCGCACCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTAt                                                                                                                                                <  2:2380980/148‑1 (MQ=255)
       cGCACCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCa                                                                                                                                            <  1:2243641/148‑1 (MQ=255)
         cACCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAAt                                                                                                                                          >  2:1822896/1‑148 (MQ=255)
         cACCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAAt                                                                                                                                          >  1:372029/1‑148 (MQ=255)
         cACCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAAt                                                                                                                                          >  1:372569/1‑148 (MQ=255)
         cACCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAAt                                                                                                                                          >  1:1002205/1‑148 (MQ=255)
          aCCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATg                                                                                                                                         <  2:2990474/148‑1 (MQ=255)
           cccAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGt                                                                                                                                        <  1:248707/148‑1 (MQ=255)
             cAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTg                                                                                                                                      <  1:1017348/148‑1 (MQ=255)
             cAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTg                                                                                                                                      <  1:3179727/148‑1 (MQ=255)
             cAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTg                                                                                                                                      >  2:1429392/1‑148 (MQ=255)
             cAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTg                                                                                                                                      >  2:2245426/1‑148 (MQ=255)
              aGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGa                                                                                                                                     <  1:2884897/148‑1 (MQ=255)
               gTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGAt                                                                                                                                    <  2:621054/148‑1 (MQ=255)
                    aGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGa                                                                                                                               <  1:2696506/148‑1 (MQ=255)
                       ttACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAg                                                                                                                            >  1:833315/1‑148 (MQ=255)
                             tGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAg                                                                                                                      >  2:635052/1‑148 (MQ=255)
                             tGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAg                                                                                                                      >  2:2622594/1‑148 (MQ=255)
                               cAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGcc                                                                                                                    <  2:2816284/148‑1 (MQ=255)
                               cAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGcc                                                                                                                    <  2:2488899/148‑1 (MQ=255)
                                  gCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTc                                                                                                                 <  2:188895/148‑1 (MQ=255)
                                  gCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTc                                                                                                                 <  2:189405/148‑1 (MQ=255)
                                    cACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTg                                                                                                               >  2:2071478/1‑148 (MQ=255)
                                      cGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTggg                                                                                                             >  1:2576138/1‑148 (MQ=255)
                                      cGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGAGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTggg                                                                                                             >  1:1787825/1‑148 (MQ=255)
                                         cACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGAt                                                                                                          >  1:3496823/1‑148 (MQ=255)
                                            tGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAg                                                                                                       >  1:903363/1‑148 (MQ=255)
                                             ggATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGc                                                                                                      <  2:41999/148‑1 (MQ=255)
                                             ggATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGc                                                                                                      <  2:297187/148‑1 (MQ=255)
                                              gATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCa                                                                                                     <  1:2432286/148‑1 (MQ=255)
                                               aTGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCAt                                                                                                    >  2:2155896/1‑148 (MQ=255)
                                               aTGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGAGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCAt                                                                                                    >  2:2263626/1‑148 (MQ=255)
                                                tGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATg                                                                                                   <  1:3034040/148‑1 (MQ=255)
                                                     tCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGa                                                                                              >  2:1946315/1‑148 (MQ=255)
                                                        gTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGt                                                                                           <  1:374869/148‑1 (MQ=255)
                                                        gTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGt                                                                                           <  1:838693/148‑1 (MQ=255)
                                                        gTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGt                                                                                           >  1:2856042/1‑148 (MQ=255)
                                                        gTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGt                                                                                           >  2:2700291/1‑148 (MQ=255)
                                                         tACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTc                                                                                          >  1:496522/1‑148 (MQ=255)
                                                          aCCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCg                                                                                         <  2:1476010/148‑1 (MQ=255)
                                                          aCCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTAATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCg                                                                                         <  1:2146862/148‑1 (MQ=255)
                                                                   ccACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGccc                                                                                <  2:2120814/148‑1 (MQ=255)
                                                                     acaGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTg                                                                              >  1:1170916/1‑148 (MQ=255)
                                                                         ttGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACTCCCTGTACa                                                                          >  1:2542199/1‑148 (MQ=255)
                                                                         ttGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACa                                                                          <  1:3178994/148‑1 (MQ=255)
                                                                           gATCGAAAAACACCAGCAACAGCCACACCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGCTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAAc                                                                        <  2:2166583/148‑1 (MQ=255)
                                                                            aTGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAAcg                                                                       >  1:804527/1‑148 (MQ=255)
                                                                              ggAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAAcgcg                                                                     <  1:2910143/148‑1 (MQ=255)
                                                                                aaaaaCTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGtt                                                                   >  1:89531/1‑148 (MQ=255)
                                                                                 aaaaCTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGttt                                                                  <  1:1417058/148‑1 (MQ=255)
                                                                                  aaaCTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTg                                                                 <  1:794612/148‑1 (MQ=255)
                                                                                  aaaCTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTg                                                                 <  1:749407/148‑1 (MQ=255)
                                                                                   aaCTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTgg                                                                >  1:265360/1‑148 (MQ=255)
                                                                                    aCTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGt                                                               >  1:1639581/1‑148 (MQ=255)
                                                                                     cTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTc                                                              <  2:326090/148‑1 (MQ=255)
                                                                                      tCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCa                                                             >  2:1709698/1‑148 (MQ=255)
                                                                                      tCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCa                                                             >  2:1710238/1‑148 (MQ=255)
                                                                                         aGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCg                                                          <  2:3072108/148‑1 (MQ=255)
                                                                                          gCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCgg                                                         <  1:2442812/148‑1 (MQ=255)
                                                                                          gCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCgg                                                         >  2:2949681/1‑148 (MQ=255)
                                                                                          gCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCgg                                                         <  2:1006836/148‑1 (MQ=255)
                                                                                           cAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGc                                                        >  2:2654683/1‑148 (MQ=255)
                                                                                             aCTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTg                                                      <  1:1822896/148‑1 (MQ=255)
                                                                                             aCTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTg                                                      <  1:1076015/148‑1 (MQ=255)
                                                                                              cTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGa                                                     >  2:665137/1‑148 (MQ=255)
                                                                                                gCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGAtt                                                   >  2:2976281/1‑148 (MQ=255)
                                                                                                gCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGAtt                                                   >  2:1199046/1‑148 (MQ=255)
                                                                                                 ccACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGAttt                                                  <  2:1017393/148‑1 (MQ=255)
                                                                                                 ccACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGAttt                                                  <  1:88574/148‑1 (MQ=255)
                                                                                                    caGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCa                                               >  1:1946254/1‑148 (MQ=255)
                                                                                                    caGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGAGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCa                                               <  1:1159956/148‑1 (MQ=255)
                                                                                                      gCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACt                                             <  1:2226025/148‑1 (MQ=255)
                                                                                                       cTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGAtt                                                   >  1:2733618/1‑141 (MQ=255)
                                                                                                        tttCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTAt                                           <  2:833307/148‑1 (MQ=255)
                                                                                                        tttCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTAt                                           <  1:1549204/148‑1 (MQ=255)
                                                                                                          tCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTAttt                                         <  2:3385575/148‑1 (MQ=255)
                                                                                                           cTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGAGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTta                                        <  2:2241119/148‑1 (MQ=255)
                                                                                                            tGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTtat                                       >  1:203801/1‑148 (MQ=255)
                                                                                                            tGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTtat                                       >  1:1078362/1‑148 (MQ=255)
                                                                                                             gATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTtata                                      >  2:69396/1‑148 (MQ=255)
                                                                                                             gATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTtata                                      <  2:1582310/148‑1 (MQ=255)
                                                                                                                gTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGGTTGGTCAGCGGCTGATTTCCACTATTTATACTc                                   >  2:1454975/1‑148 (MQ=255)
                                                                                                                gTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGCCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTc                                   <  1:1491997/148‑1 (MQ=255)
                                                                                                                 tctcCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCa                                  <  2:2739792/148‑1 (MQ=255)
                                                                                                                 tctcCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTGTGGTCAGCGGCTGATTTCCACTATTTATACTCa                                  >  2:1549382/1‑148 (MQ=255)
                                                                                                                    ccAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGc                               >  2:2852949/1‑148 (MQ=255)
                                                                                                                    ccAGCGACTGGCAGGACAAAGGGCTGGTGGAGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGc                               <  1:3364874/148‑1 (MQ=255)
                                                                                                                       gCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCaa                            >  2:2212799/1‑148 (MQ=255)
                                                                                                                           cTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAg                        <  1:120254/148‑1 (MQ=255)
                                                                                                                            tGGCAGGACAAAGGGCTGGTGGAGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGt                       <  2:1524273/148‑1 (MQ=255)
                                                                                                                               cAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATc                    >  1:2355318/1‑148 (MQ=255)
                                                                                                                                aGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATcg                   <  2:1078548/148‑1 (MQ=255)
                                                                                                                                  gACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAGCCAGTATcgcg                 >  1:1101330/1‑148 (MQ=255)
                                                                                                                                  gACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATcgcg                 <  1:3283373/148‑1 (MQ=255)
                                                                                                                                  gACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATcgcg                 >  1:705553/1‑148 (MQ=255)
                                                                                                                                  gACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATcgcg                 >  1:2567921/1‑148 (MQ=255)
                                                                                                                                             ggtggCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAg      <  2:903363/148‑1 (MQ=255)
                                                                                                                                              gtggCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGc     >  2:3546257/1‑148 (MQ=255)
                                                                                                                                              gtggCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGc     >  2:3546442/1‑148 (MQ=255)
                                                                                                                                                 gCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGcaca  >  1:1522758/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
GGTCAGCCGCACCCAGTACCAGTTTACCTTGCAGGCCACGTCACTGGATGCGCTCAGTACCTGGGTGCCACAGTTGATGGAAAAACTCCAGCAACTGCCACAGCTTTCTGATGTCTCCAGCGACTGGCAGGACAAAGGGCTGGTGGCGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACA  >  CP009273/2150750‑2151042

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: