Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,147,980 C→A 9.4% R162S (CGC→AGC)  mdtA → multidrug efflux system, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,147,9800CA 9.4% 756.1 / 34.7 212R162S (CGC→AGC) mdtAmultidrug efflux system, subunit A
Reads supporting (aligned to +/- strand):  ref base C (121/71);  new base A (9/11);  total (130/82)
Fisher's exact test for biased strand distribution p-value = 1.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 214 total aligned reads displayed.

TTTACTGGCAGAAATTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGG  >  CP009273/2147835‑2148127
                                                                                                                                                 |                                                                                                                                                   
 ttACTGGCAGAAATTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGcc                                                                                                                                                  >  1:2692693/1‑148 (MQ=255)
  tACTGGCAGAAATTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCa                                                                                                                                                 <  1:13623/148‑1 (MQ=255)
      ggCAGAAATTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAg                                                                                                                                             >  2:2737579/1‑148 (MQ=255)
        cAGAAATTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCt                                                                                                                                           <  1:2345905/148‑1 (MQ=255)
          gAAATTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTgg                                                                                                                                         >  1:1161518/1‑148 (MQ=255)
             aTTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATg                                                                                                                                      >  1:1834194/1‑148 (MQ=255)
             aTTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATg                                                                                                                                      >  1:1291037/1‑148 (MQ=255)
                  ccccAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCaa                                                                                                                                 >  2:3163997/1‑148 (MQ=255)
                   cccAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCAGCCAGGAGCTGGATGCCCAAc                                                                                                                                <  2:2460371/148‑1 (MQ=255)
                     cagccagTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCCAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAg                                                                                                                              >  2:605260/1‑148 (MQ=255)
                      agccagTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAgg                                                                                                                             >  1:1195480/1‑148 (MQ=255)
                        ccagTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCCGGAGCTGGATGCCCAACAGgcg                                                                                                                           >  2:548757/1‑148 (MQ=255)
                        ccagTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGgcg                                                                                                                           <  1:1136019/148‑1 (MQ=255)
                         cagTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGgcgc                                                                                                                          <  1:3105172/148‑1 (MQ=255)
                            ttCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTgg                                                                                                                       <  1:707333/148‑1 (MQ=255)
                            ttCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTgg                                                                                                                       >  1:1055939/1‑148 (MQ=255)
                             tCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGt                                                                                                                      >  1:2318573/1‑148 (MQ=255)
                              cAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAAAAACTGGCAAAAACAAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTc                                                                                                                     >  1:1034990/1‑148 (MQ=255)
                                 aGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGt                                                                                                                  >  2:2272165/1‑148 (MQ=255)
                                 aGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGt                                                                                                                  >  1:3204175/1‑148 (MQ=255)
                                     gCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGaaa                                                                                                              >  1:1882848/1‑148 (MQ=255)
                                        ttAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCg                                                                                                           >  2:2495507/1‑148 (MQ=255)
                                        ttAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCg                                                                                                           >  2:2127682/1‑148 (MQ=255)
                                         tAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGa                                                                                                          <  1:1295220/148‑1 (MQ=255)
                                         tAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGa                                                                                                          <  2:683082/148‑1 (MQ=255)
                                         tAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGa                                                                                                          <  1:569314/148‑1 (MQ=255)
                                         tAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGa                                                                                                          <  1:1066610/148‑1 (MQ=255)
                                         tAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCAGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGa                                                                                                          >  2:397330/1‑148 (MQ=255)
                                          aGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGaa                                                                                                         <  2:3061239/148‑1 (MQ=255)
                                           gCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCAGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAg                                                                                                        <  2:2073010/148‑1 (MQ=255)
                                            cacaAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAgg                                                                                                       >  2:908015/1‑148 (MQ=255)
                                                 gCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTGTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCa                                                                                                  <  2:3303005/148‑1 (MQ=255)
                                                 gCCCAGGGCCAACTGGAAAAAGATAAAGCCACGCTTGCCAACGCACGCCGAGACCTGGCGCGTTATCAACAACTGGCAAAAACCAAACTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCa                                                                                                  >  2:1045751/1‑148 (MQ=255)
                                                    cAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTa                                                                                               <  1:1109010/148‑1 (MQ=255)
                                                    cAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTa                                                                                               <  1:2949613/148‑1 (MQ=255)
                                                     aGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTaa                                                                                              <  1:1656530/148‑1 (MQ=255)
                                                       ggCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAgg                                                                                            >  1:2826784/1‑148 (MQ=255)
                                                       ggCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAgg                                                                                            <  2:1386381/148‑1 (MQ=255)
                                                         ccAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCt                                                                                          <  2:2752702/148‑1 (MQ=255)
                                                              tGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACATCTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCtgatga                                                                                     <  2:796829/148‑1 (MQ=255)
                                                              tGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCAGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCtgatga                                                                                     >  1:159841/1‑148 (MQ=255)
                                                                gCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGaag                                                                                   >  1:1813898/1‑148 (MQ=255)
                                                                    aaGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGaagcaag                                                                               <  2:844120/148‑1 (MQ=255)
                                                                    aaGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGaagcaag                                                                               <  2:842938/148‑1 (MQ=255)
                                                                    aaGATAAAGCCACGCTTGCAAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGaagcaag                                                                               >  2:779133/1‑148 (MQ=255)
                                                                      gATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCACTCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTTGTCAGTGAAACCGAAGGCACCATTAAGGCTTATTAAGCAAGCg                                                                             <  2:1027355/148‑1 (MQ=255)
                                                                      gATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCg                                                                             <  2:2957508/148‑1 (MQ=255)
                                                                      gATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCg                                                                             <  1:1395370/148‑1 (MQ=255)
                                                                        tAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGtt                                                                           >  1:1745482/1‑148 (MQ=255)
                                                                         aaaGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGTCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTg                                                                          >  1:2479291/1‑148 (MQ=255)
                                                                           aGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGcc                                                                        >  1:551933/1‑148 (MQ=255)
                                                                           aGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGcc                                                                        >  2:2472044/1‑148 (MQ=255)
                                                                             ccACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAg                                                                      >  2:3017947/1‑148 (MQ=255)
                                                                                        aaCGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAATAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCCGCGCGCAGCTGc                                                           >  2:3233082/1‑148 (MQ=255)
                                                                                               gccgTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGa                                                    <  2:1688688/148‑1 (MQ=255)
                                                                                                ccgTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGAc                                                   >  1:3042936/1‑148 (MQ=255)
                                                                                                  gTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTg                                                 <  1:2272165/148‑1 (MQ=255)
                                                                                                    gACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGa                                               >  1:2591242/1‑148 (MQ=255)
                                                                                                    gACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCAGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGa                                               <  2:2385064/148‑1 (MQ=255)
                                                                                                      ccTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGc                                             >  2:2119285/1‑148 (MQ=255)
                                                                                                          gcgcGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGa                                         <  2:1809557/148‑1 (MQ=255)
                                                                                                          gcgcGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGa                                         >  2:529653/1‑148 (MQ=255)
                                                                                                          gcgaGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGa                                         >  2:417355/1‑148 (MQ=255)
                                                                                                           cgcgTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGAt                                        <  1:800282/148‑1 (MQ=255)
                                                                                                             cgTTATCAACAACTGGCAAAAACCAATCTCGTTTCCAGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTa                                      >  2:1181289/1‑148 (MQ=255)
                                                                                                               ttATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTAcc                                    >  2:1736867/1‑148 (MQ=255)
                                                                                                               ttATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTAcc                                    >  2:891784/1‑148 (MQ=255)
                                                                                                               ttATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTAcc                                    >  1:2415411/1‑148 (MQ=255)
                                                                                                               ttATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTAcc                                    >  2:1642326/1‑148 (MQ=255)
                                                                                                                  tCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCa                                 >  2:2219576/1‑148 (MQ=255)
                                                                                                                  tCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCa                                 >  2:2697934/1‑148 (MQ=255)
                                                                                                                   caacaaCTGGCAAAAACCACTCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAc                                <  1:1707054/148‑1 (MQ=255)
                                                                                                                   caacaaCTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAc                                >  2:1149508/1‑148 (MQ=255)
                                                                                                                    aacaacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAcc                               >  1:1523419/1‑148 (MQ=255)
                                                                                                                    aacaacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCAcc                               >  2:1664502/1‑148 (MQ=255)
                                                                                                                     acaacTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCa                              >  1:1637285/1‑148 (MQ=255)
                                                                                                                        acTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTc                           >  2:1058111/1‑148 (MQ=255)
                                                                                                                         cTGGCAAAACCCACGCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCg                          <  2:2941629/148‑1 (MQ=255)
                                                                                                                         cTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCg                          >  1:261435/1‑148 (MQ=255)
                                                                                                                          tGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGa                         <  2:2736474/148‑1 (MQ=255)
                                                                                                                             cAAAAACCAATCTCGTTTCCAGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATgg                      <  1:397330/148‑1 (MQ=255)
                                                                                                                                  aCCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTcgcg                 <  2:3204175/148‑1 (MQ=255)
                                                                                                                                  aCCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGAGGGTcgcg                 >  1:50156/1‑148 (MQ=255)
                                                                                                                                   ccAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGt                <  2:2672095/148‑1 (MQ=255)
                                                                                                                                      aTCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTgg             >  2:1470558/1‑148 (MQ=255)
                                                                                                                                      aTCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTgg             >  2:3279342/1‑148 (MQ=255)
                                                                                                                                      aTCTCGTTTCCAGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTgg             <  1:3023164/148‑1 (MQ=255)
                                                                                                                                       tctcGTTTCCCGCCAGGAGCTGGCTGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGt            <  2:1576280/148‑1 (MQ=255)
                                                                                                                                       tctcGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGt            <  1:3311114/148‑1 (MQ=255)
                                                                                                                                       tctcGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGt            >  1:2168236/1‑148 (MQ=255)
                                                                                                                                       tctcGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGt            >  2:2640262/1‑148 (MQ=255)
                                                                                                                                         tcGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGtct          <  1:1000846/148‑1 (MQ=255)
                                                                                                                                          cGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGtctc         >  2:2424588/1‑148 (MQ=255)
                                                                                                                                          cGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGtctc         >  1:935256/1‑148 (MQ=255)
                                                                                                                                          cGTTTACCGCCAGGAGCTGGATGAAAAACACGCGCTGGTCAGTGAAACCGAAGGCACCAGTAAAGGTGATGAAGCAAGCGTTGCCATCGCTCAGCTGCAACACGACTGGAGACGGATTACCGCAACAGTCGATGTTCGCGTTGGtctc         >  2:2849962/1‑148 (MQ=255)
                                                                                                                                            tttCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCaa       >  2:1920314/1‑148 (MQ=255)
                                                                                                                                            tttCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCaa       >  2:1058736/1‑148 (MQ=255)
                                                                                                                                             ttCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAg      >  2:830872/1‑148 (MQ=255)
                                                                                                                                              tCCAGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGc     >  2:2997338/1‑148 (MQ=255)
                                                                                                                                                ccGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCAGTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAg   <  2:487304/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
TTTACTGGCAGAAATTGACCCCAGCCAGTTCAAAGTTGCATTAGCACAAGCCCAGGGCCAACTGGCAAAAGATAAAGCCACGCTTGCCAACGCCCGCCGTGACCTGGCGCGTTATCAACAACTGGCAAAAACCAATCTCGTTTCCCGCCAGGAGCTGGATGCCCAACAGGCGCTGGTCAGTGAAACCGAAGGCACCATTAAGGCTGATGAAGCAAGCGTTGCCAGCGCGCAGCTGCAACTCGACTGGAGCCGGATTACCGCACCAGTCGATGGTCGCGTTGGTCTCAAGCAGG  >  CP009273/2147835‑2148127

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: