Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,796,148 C→T Y243Y (TAC→TAT)  RPPX_08005 → integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,796,1480CT99.6% 650.7 / ‑6.4 261Y243Y (TAC→TAT) RPPX_08005integrase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (145/115);  total (145/116)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 268 total aligned reads displayed.

AAAGCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTACTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGGCC  >  CP009974/1796001‑1796295
                                                                                                                                                   |                                                                                                                                                   
aaaGCAGAGCATGTCTCGAAAAGGCGATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACtattac                                                                                                                                                  >  1:3976210/1‑151 (MQ=14)
 aaGCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACa                                                                                                                                                 >  1:3698035/1‑151 (MQ=14)
  aGCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACaa                                                                                                                                                <  2:5512547/151‑1 (MQ=17)
     agagCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCa                                                                                                                                             >  1:5612463/1‑151 (MQ=17)
      gagCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCAt                                                                                                                                            >  2:6138277/1‑151 (MQ=17)
       agCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATg                                                                                                                                           <  2:4495547/151‑1 (MQ=17)
       agCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATg                                                                                                                                           >  2:1240424/1‑151 (MQ=14)
         cATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGTGTCTGGTTGAGTACATTCACTATTACAAccctcat                                                                                                                                         >  2:6481275/1‑146 (MQ=255)
         cATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGAt                                                                                                                                         >  1:1100702/1‑151 (MQ=14)
          aTGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATc                                                                                                                                        <  2:770091/151‑1 (MQ=14)
           tGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCg                                                                                                                                       >  2:116969/1‑151 (MQ=17)
             tctcGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCa                                                                                                                                     <  2:6343717/151‑1 (MQ=17)
                cGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCa                                                                                                                                  <  2:26468/151‑1 (MQ=14)
                 gAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCaa                                                                                                                                 <  1:3005821/151‑1 (MQ=14)
                       gCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTgag                                                                                                                           >  2:527819/1‑151 (MQ=17)
                       gCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTgag                                                                                                                           >  2:2861561/1‑151 (MQ=17)
                         aaTTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGCAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGAc                                                                                                                         >  2:5417853/1‑151 (MQ=14)
                              ccTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAAt                                                                                                                    >  2:2226122/1‑151 (MQ=17)
                               cTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATg                                                                                                                   <  1:1220772/151‑1 (MQ=17)
                               cTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATg                                                                                                                   <  1:1649478/151‑1 (MQ=17)
                                 ggACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGc                                                                                                                 >  2:3308206/1‑151 (MQ=14)
                                   aCAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCt                                                                                                               >  2:1855784/1‑151 (MQ=17)
                                      aTGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAg                                                                                                            <  1:3308206/151‑1 (MQ=17)
                                        gCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGc                                                                                                                               >  2:2013818/1‑130 (MQ=18)
                                        gCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGc                                                                                                                               <  1:2013818/130‑1 (MQ=18)
                                        gCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGcc                                                                                                          >  1:3081167/1‑151 (MQ=17)
                                          ggCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCt                                                                                                        >  2:36050/1‑151 (MQ=17)
                                             cATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCtgtg                                                                                                     <  1:489915/151‑1 (MQ=17)
                                                ggAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGGATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGaa                                                                                                  >  2:1064504/1‑151 (MQ=14)
                                                  aaaGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGa                                                                                                             <  2:57382/138‑1 (MQ=18)
                                                  aaaGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGa                                                                                                             >  1:57382/1‑138 (MQ=18)
                                                  aaaGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATa                                                                                                <  1:6322822/151‑1 (MQ=17)
                                                  aaaGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATa                                                                                                <  2:4533983/151‑1 (MQ=17)
                                                   aaGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATAc                                                                                               <  1:1912431/151‑1 (MQ=14)
                                                    aGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTgaga                                                                                                                          <  2:5240159/123‑1 (MQ=18)
                                                    aGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACa                                                                                              >  1:158307/1‑151 (MQ=14)
                                                      cttcttCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAgg                                                                                            >  1:3768039/1‑151 (MQ=14)
                                                      cttcttCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAgg                                                                                            <  1:361568/151‑1 (MQ=17)
                                                        tcttcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGc                                                                                          >  2:3442500/1‑151 (MQ=17)
                                                          ttcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTc                                                                                        >  2:1603179/1‑151 (MQ=17)
                                                           tcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCACTGGCCTGATCCCGGTGGAATACAGGGCTCa                                                                                       >  2:3018127/1‑151 (MQ=255)
                                                             ggCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAgg                                                                                     >  2:1788502/1‑151 (MQ=17)
                                                              gCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGc                                                                                    >  2:584729/1‑151 (MQ=17)
                                                              gCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGc                                                                                    >  2:4461091/1‑151 (MQ=17)
                                                               cACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCt                                                                                   >  2:1776736/1‑151 (MQ=17)
                                                                 cGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGc                                                                                 >  1:5490667/1‑151 (MQ=17)
                                                                     cAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCt                                                                             >  2:4696630/1‑151 (MQ=17)
                                                                      aaaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTTACTATTACAACCATGATCGCATCAAGCTAAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCtt                                                                            >  1:1379275/1‑151 (MQ=9)
                                                                      aaaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCtt                                                                            >  1:5379624/1‑151 (MQ=17)
                                                                       aaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTa                                                                           <  2:40030/151‑1 (MQ=14)
                                                                       aaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTa                                                                           <  1:6533933/151‑1 (MQ=14)
                                                                         tCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTaaa                                                                         >  1:5723416/1‑151 (MQ=17)
                                                                         tCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTaaa                                                                         >  1:1028911/1‑151 (MQ=17)
                                                                          cTGAGTTTTTCTACCTGACGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTaaaa                                                                        >  2:4329910/1‑151 (MQ=14)
                                                                           tGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAAc                                                                       <  2:2406500/151‑1 (MQ=17)
                                                                               tttttCTACCTGAGGTGTTTGGAGAGCATTGATGAGCTAAGAGCGGGTCTGGATAAGTACATTCACTATTACAACAATGAGCGCATCAAGCTGAGACTCAATGGCCTGAACCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTc                                                                   <  2:2438614/151‑1 (MQ=11)
                                                                               tttttCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTc                                                                   >  2:6190927/1‑151 (MQ=17)
                                                                               tttttCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTc                                                                   <  1:1753959/151‑1 (MQ=14)
                                                                                ttttCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTcc                                                                  >  1:5182582/1‑151 (MQ=17)
                                                                                     tACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACtt                                                             >  1:2859602/1‑151 (MQ=14)
                                                                                      aCCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACttt                                                            <  2:4070010/151‑1 (MQ=14)
                                                                                       ccTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACtttt                                                           >  1:4510779/1‑151 (MQ=17)
                                                                                           aaGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTgggg                                                       <  2:1670441/151‑1 (MQ=14)
                                                                                             gtgtTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTgggggg                                                     <  2:2697601/151‑1 (MQ=14)
                                                                                                        gCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATAcc                                          >  2:4808125/1‑151 (MQ=14)
                                                                                                        gCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATAcc                                          >  2:1449826/1‑151 (MQ=14)
                                                                                                        gCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATAcc                                          >  2:2150265/1‑151 (MQ=14)
                                                                                                          aTTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtc                                        >  2:6396010/1‑150 (MQ=14)
                                                                                                           ttGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtgt                                       <  2:502876/151‑1 (MQ=17)
                                                                                                           ttGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtgt                                       >  2:3670581/1‑151 (MQ=14)
                                                                                                            tgatgaGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtgtg                                      >  2:3341065/1‑151 (MQ=17)
                                                                                                              atgaGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTggg                                    <  1:5042680/151‑1 (MQ=17)
                                                                                                                 aGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGc                                 <  1:2665961/151‑1 (MQ=17)
                                                                                                                   cTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCgg                               >  1:2241907/1‑151 (MQ=17)
                                                                                                                     gagaGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGc                             >  2:6050566/1‑151 (MQ=14)
                                                                                                                      agagCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCGGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCa                            >  2:555362/1‑151 (MQ=17)
                                                                                                                       gagCGGCTCTCGAGCAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCaa                           <  2:3152568/151‑1 (MQ=255)
                                                                                                                        agCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGAGGGGCAGTCCATACCTGTGGGAGCGGGCAAg                          >  2:3936098/1‑151 (MQ=18)
                                                                                                                           gggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGccc                       >  1:1428571/1‑151 (MQ=37)
                                                                                                                           gggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGccc                       <  1:4829136/151‑1 (MQ=37)
                                                                                                                           gggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGccc                       <  1:3341065/151‑1 (MQ=37)
                                                                                                                            ggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCcg                      >  1:2274511/1‑151 (MQ=38)
                                                                                                                             gTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCcgc                     <  1:5936597/151‑1 (MQ=38)
                                                                                                                               cTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCGCGCGa                   >  2:4162456/1‑151 (MQ=255)
                                                                                                                               cTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGa                   <  1:5322577/151‑1 (MQ=37)
                                                                                                                                tGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGaa                  <  1:5396875/151‑1 (MQ=255)
                                                                                                                                 ggATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAg                 <  1:2344717/151‑1 (MQ=255)
                                                                                                                                  gATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGGGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAgg                <  1:47413/151‑1 (MQ=255)
                                                                                                                                      aGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCg            >  2:1540281/1‑151 (MQ=255)
                                                                                                                                      aGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCg            >  1:4572977/1‑151 (MQ=255)
                                                                                                                                        tACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCa          >  1:6472488/1‑151 (MQ=255)
                                                                                                                                         aCATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAg         <  1:6180296/151‑1 (MQ=255)
                                                                                                                                          cATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgc        <  1:452457/151‑1 (MQ=255)
                                                                                                                                          cATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgc        <  1:3184406/151‑1 (MQ=255)
                                                                                                                                           aTTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgcg       >  2:5252450/1‑151 (MQ=255)
                                                                                                                                           aTTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgcc       >  1:4082854/1‑150 (MQ=255)
                                                                                                                                            ttCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgcgc      >  2:4905035/1‑151 (MQ=255)
                                                                                                                                            ttCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgcgc      <  1:80776/151‑1 (MQ=255)
                                                                                                                                            ttCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgcgc      >  2:941154/1‑151 (MQ=255)
                                                                                                                                               actattACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGGc   >  2:4359469/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAAGCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTACTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGGCC  >  CP009974/1796001‑1796295

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑