Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,168,668 T→A 9.4% T30S (ACT→TCT)  parE ← DNA topoisomerase IV, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,168,6680TA 9.4% 671.5 / 47.0 225T30S (ACT→TCT) parEDNA topoisomerase IV, subunit B
Reads supporting (aligned to +/- strand):  ref base T (100/103);  new base A (12/9);  total (112/113)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 228 total aligned reads displayed.

GAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGAC  >  CP009273/3168521‑3168813
                                                                                                                                                   |                                                                                                                                                 
taaTATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGt                                                                                                                                                   <  2:2706724/147‑1 (MQ=255)
gAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGt                                                                                                                                                   >  2:425448/1‑148 (MQ=255)
  atatCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTgg                                                                                                                                                 >  1:2261989/1‑148 (MQ=255)
     tCCCCCGGCATCCCGCGCCAATCGTAAATAACTACTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTAt                                                                                                                                              <  1:3056144/148‑1 (MQ=255)
     tCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTAt                                                                                                                                              <  1:1437040/148‑1 (MQ=255)
             caccccGCCCCCAGCGAAAATAACTACACACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATAc                                                                                                                                      <  1:1564876/145‑1 (MQ=255)
             cATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATAc                                                                                                                                      <  2:3522856/148‑1 (MQ=255)
              aTCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACa                                                                                                                                     <  2:2167252/148‑1 (MQ=255)
              aTCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACa                                                                                                                                     >  2:2803956/1‑148 (MQ=255)
               tCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACAt                                                                                                                                    >  2:2325958/1‑148 (MQ=255)
               tCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACAt                                                                                                                                    >  1:2049979/1‑148 (MQ=255)
               tCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACAt                                                                                                                                    >  2:2111157/1‑148 (MQ=255)
                       ccATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGAc                                                                                                                            >  1:2425565/1‑148 (MQ=255)
                         aTCGTCAATAACTACACACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGAcgg                                                                                                                          <  1:599096/148‑1 (MQ=255)
                              cAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGgcggcg                                                                                                                     >  1:1968623/1‑148 (MQ=255)
                                    cTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCgg                                                                                                               <  1:2038996/148‑1 (MQ=255)
                                     ttCTCACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGc                                                                                                              <  2:2316455/148‑1 (MQ=255)
                                     ttCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGc                                                                                                              >  2:1306703/1‑148 (MQ=255)
                                     ttCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGc                                                                                                              <  2:2481440/148‑1 (MQ=255)
                                     ttCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAAGGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGc                                                                                                              >  2:86158/1‑148 (MQ=255)
                                       cTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTc                                                                                                            >  2:94545/1‑148 (MQ=255)
                                          aCGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAg                                                                                                         >  2:2585801/1‑148 (MQ=255)
                                           cGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGc                                                                                                        >  1:2239049/1‑148 (MQ=255)
                                           cGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGc                                                                                                        >  2:82494/1‑148 (MQ=255)
                                              cTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCCAGCCCg                                                                                                     >  2:3528165/1‑148 (MQ=255)
                                              cTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCg                                                                                                     >  2:1181769/1‑148 (MQ=255)
                                              cTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCg                                                                                                     >  2:1885502/1‑148 (MQ=255)
                                                ggTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGt                                                                                                   <  2:3227566/148‑1 (MQ=255)
                                                ggTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGt                                                                                                   <  1:3418157/148‑1 (MQ=255)
                                                ggTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGt                                                                                                   <  2:1460699/148‑1 (MQ=255)
                                                  tCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGa                                                                                                 <  2:373339/148‑1 (MQ=255)
                                                   cAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAg                                                                                                >  1:713131/1‑148 (MQ=255)
                                                   cAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAg                                                                                                >  1:1459133/1‑148 (MQ=255)
                                                    aGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGt                                                                                               <  1:2353449/148‑1 (MQ=255)
                                                     gCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTa                                                                                              <  1:2611158/148‑1 (MQ=255)
                                                     gCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTa                                                                                              <  2:1881020/148‑1 (MQ=255)
                                                     gCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTa                                                                                              <  1:2651474/148‑1 (MQ=255)
                                                      cATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTAc                                                                                             >  1:1112773/1‑148 (MQ=255)
                                                      cATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTAc                                                                                             >  2:2041976/1‑148 (MQ=255)
                                                      cATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTAc                                                                                             >  1:1113240/1‑148 (MQ=255)
                                                          tAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCa                                                                                         <  1:188563/148‑1 (MQ=255)
                                                           aaaaTAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCaa                                                                                        >  1:1013835/1‑148 (MQ=255)
                                                            aaaTAACGTCCACGCGTTTTGCGGGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAAt                                                                                       <  1:1492937/148‑1 (MQ=255)
                                                             aaTAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATg                                                                                      >  1:643272/1‑148 (MQ=255)
                                                             aaTAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATg                                                                                      <  2:3251918/148‑1 (MQ=255)
                                                                aaCGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCa                                                                                   >  2:2853703/1‑148 (MQ=255)
                                                                aaCGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCa                                                                                   >  1:1245176/1‑148 (MQ=255)
                                                                 aCGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCAt                                                                                  <  2:564604/148‑1 (MQ=255)
                                                                 aCGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCAt                                                                                  <  1:421207/148‑1 (MQ=255)
                                                                    tCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAg                                                                               <  1:372703/148‑1 (MQ=255)
                                                                       aCGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGt                                                                            <  1:334979/148‑1 (MQ=255)
                                                                        cgcgTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGtt                                                                           >  2:901310/1‑148 (MQ=255)
                                                                         gcgTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTta                                                                          <  2:2398311/148‑1 (MQ=255)
                                                                            ttttGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATaa                                                                       >  2:2208349/1‑148 (MQ=255)
                                                                             tttGAGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAg                                                                      <  1:2529563/148‑1 (MQ=255)
                                                                               tGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGtt                                                                    <  2:2610193/148‑1 (MQ=255)
                                                                                    gACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGt                                                               <  2:3512727/148‑1 (MQ=255)
                                                                                     aCCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTc                                                              >  2:717140/1‑148 (MQ=255)
                                                                                      cccGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCa                                                             <  1:2336608/148‑1 (MQ=255)
                                                                                            aGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGttt                                                       >  2:2222797/1‑148 (MQ=255)
                                                                                             gTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTa                                                      <  1:3127248/148‑1 (MQ=255)
                                                                                             gTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTa                                                      >  2:459798/1‑148 (MQ=255)
                                                                                               gCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAg                                                    >  1:2201341/1‑148 (MQ=255)
                                                                                                cTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGt                                                   <  2:2261989/148‑1 (MQ=255)
                                                                                                cTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGt                                                   >  2:2136199/1‑148 (MQ=255)
                                                                                                    aTCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGt                                               <  2:447260/148‑1 (MQ=255)
                                                                                                      ccACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTCTAAGTTTGCGTCATGGTTTAAGTTAGTaa                                             >  2:3021045/1‑148 (MQ=255)
                                                                                                      ccACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTGTAAGTTAGTaa                                             >  2:1918200/1‑148 (MQ=255)
                                                                                                           cTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGa                                        <  1:2094416/148‑1 (MQ=255)
                                                                                                           cTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGa                                        >  1:1131425/1‑148 (MQ=255)
                                                                                                             gTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGt                                      <  1:861637/148‑1 (MQ=255)
                                                                                                                aTCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGa                                   <  2:2097674/148‑1 (MQ=255)
                                                                                                                aTCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGa                                   <  2:1696566/148‑1 (MQ=255)
                                                                                                                aTCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGa                                   <  2:2717219/148‑1 (MQ=255)
                                                                                                                  cAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATc                                 >  1:2482394/1‑148 (MQ=255)
                                                                                                                   aaggACTTCTTGCCCCAACTGGTTAGGGCGAGTGGTATCGGTAGACATGCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGGACCTCAATGGCATCAGCGTTATACGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCg                                <  2:2510347/145‑1 (MQ=255)
                                                                                                                   aaTGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCg                                <  2:1344311/148‑1 (MQ=255)
                                                                                                                    aTGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGt                               >  1:2703350/1‑148 (MQ=255)
                                                                                                                     tGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTc                              >  1:524471/1‑148 (MQ=255)
                                                                                                                      gACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCa                             >  2:2292833/1‑148 (MQ=255)
                                                                                                                      gACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGGTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCa                             >  2:2938908/1‑148 (MQ=255)
                                                                                                                       aCTTCTTGCCCCAAATGGTTAGGGCGAGAGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCag                            >  2:802039/1‑148 (MQ=255)
                                                                                                                        cttcttGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCaga                           >  1:3532514/1‑148 (MQ=255)
                                                                                                                            ttGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATg                       >  2:3376152/1‑148 (MQ=255)
                                                                                                                            ttGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATg                       >  2:1615233/1‑148 (MQ=255)
                                                                                                                                   aaaTGGTTAGGGCGAGTGGTATAGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGa                <  1:283144/148‑1 (MQ=255)
                                                                                                                                   aaaTGGTTAGGGCGAGAGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGa                >  2:756014/1‑148 (MQ=255)
                                                                                                                                    aaTGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGAc               >  2:1971522/1‑148 (MQ=255)
                                                                                                                                    aaTGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGAc               >  1:2771584/1‑148 (MQ=255)
                                                                                                                                     aTGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGAcc              >  2:1040559/1‑148 (MQ=255)
                                                                                                                                       ggTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCaa            >  2:2091212/1‑148 (MQ=255)
                                                                                                                                          tAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGaa         >  2:1039065/1‑148 (MQ=255)
                                                                                                                                           aGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGaaa        <  2:3190981/148‑1 (MQ=255)
                                                                                                                                           aGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGaaa        <  2:1512209/148‑1 (MQ=255)
                                                                                                                                           aGGGCGAGTGGTATCGGTATACATCCCAGTACGGCGGCGAACCGGATCAAGCCCTGTGAGTTGCTCAATGGCATCATCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGTTCGTCAGAGATGGTGCAGACATAGaaa        <  2:621944/148‑1 (MQ=255)
                                                                                                                                           aGGGCGAGCGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGaaa        <  2:1013835/148‑1 (MQ=255)
                                                                                                                                               cGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATcg    <  2:2601639/148‑1 (MQ=255)
                                                                                                                                                gAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATcga   <  2:1771299/148‑1 (MQ=255)
                                                                                                                                                gAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATcga   >  2:37605/1‑148 (MQ=255)
                                                                                                                                                 aGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGGTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATcgac  >  1:3328850/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAATATCCACCGGCATCCCGCGCCCATCGTCAATAACTTCTAACGACTGGTCAGCATGTAAAATAACGTCCACGCGTTTTGCGTGACCCGCCAGTGCTTCATCCACACTGTTATCAATGACTTCTTGCCCCAAATGGTTAGGGCGAGTGGTATCGGTATACATCCCCGGACGGCGGCGAACCGGCTCAAGCCCGGTGAGTACCTCAATGGCATCAGCGTTATAAGTTTGCGTCATGGTTTAAGTTAGTAATTCGAGTTGATCGTCAGAGATGGTGCAGACCAAGAAAATCGAC  >  CP009273/3168521‑3168813

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: