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.5% 503.7 / NA 204Y243Y (TAC→TAT) RPPX_08005integrase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (110/93);  total (111/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 209 total aligned reads displayed.

TAAAGCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTACTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGGCC  >  CP009974/1796000‑1796295
                                                                                                                                                    |                                                                                                                                                   
    gCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAAc                                                                                                                                               >  1:540809/1‑151 (MQ=14)
    gCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAAc                                                                                                                                               >  1:1261614/1‑151 (MQ=17)
     cAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAAcc                                                                                                                                              >  2:145370/1‑151 (MQ=17)
     cAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAAcc                                                                                                                                              >  1:3104808/1‑151 (MQ=17)
      agagCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCa                                                                                                                                             >  2:1093326/1‑151 (MQ=14)
      agagCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCa                                                                                                                                             <  2:1600080/151‑1 (MQ=17)
         gCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGa                                                                                                                                          <  1:2766931/151‑1 (MQ=17)
           aTGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATc                                                                                                                                        <  1:3724077/151‑1 (MQ=17)
           aTGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATc                                                                                                                                        <  2:2920177/151‑1 (MQ=17)
             gTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGc                                                                                                                                      >  1:3424028/1‑151 (MQ=17)
              tctcGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCa                                                                                                                                     <  1:3497791/151‑1 (MQ=17)
               ctcGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCAt                                                                                                                                    >  1:4881563/1‑151 (MQ=17)
                 cGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCa                                                                                                                                  <  2:71938/151‑1 (MQ=14)
                  gAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACCATTACAACCATGATCGCATCaa                                                                                                                                 <  2:4659205/151‑1 (MQ=17)
                   aaaaGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAg                                                                                                                                <  1:954511/151‑1 (MQ=17)
                   aaaaGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAg                                                                                                                                <  1:4703700/151‑1 (MQ=17)
                        gCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCAATATTACAACCACTATCGCATCAAGCTgag                                                                                                                           >  2:4269747/1‑151 (MQ=9)
                         cAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGATgagc                                                                                                                          >  2:142478/1‑150 (MQ=14)
                          aaTTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTGAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGAc                                                                                                                         >  1:3936043/1‑151 (MQ=14)
                          aaTTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTCCCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGCACATTCCCTATTACAACCATGATCGCATCAAGCTGAGAc                                                                                                                         >  1:2550666/1‑151 (MQ=9)
                          aaTTGCCTGGACAATGCGGCCATGGAAAGCTTCGTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGAc                                                                                                                         <  2:3836752/151‑1 (MQ=17)
                            ttGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTc                                                                                                                       >  2:3344903/1‑151 (MQ=17)
                            ttGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTc                                                                                                                       >  2:2470467/1‑151 (MQ=17)
                                cTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATg                                                                                                                   >  2:3795118/1‑151 (MQ=17)
                                  ggACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGc                                                                                                                 >  1:4057281/1‑151 (MQ=17)
                                    aCAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGCCTCAATGGCCt                                                                                                               >  1:4733187/1‑151 (MQ=14)
                                    aCAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCt                                                                                                               >  1:13462/1‑151 (MQ=17)
                                     cAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTg                                                                                                              >  1:2662600/1‑151 (MQ=17)
                                     cAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTg                                                                                                              <  2:1105129/151‑1 (MQ=17)
                                      aaTGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGa                                                                                                             >  2:933033/1‑151 (MQ=14)
                                      aaTGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGa                                                                                                             <  2:4744380/151‑1 (MQ=17)
                                              cATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCtgtg                                                                                                     >  2:2801193/1‑151 (MQ=17)
                                                  gAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAAt                                                                                                 >  2:2946058/1‑151 (MQ=17)
                                                  gAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAAt                                                                                                 >  2:1827690/1‑151 (MQ=17)
                                                  gAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAAt                                                                                                 <  1:246305/151‑1 (MQ=17)
                                                     aGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACa                                                                                              <  2:93143/151‑1 (MQ=17)
                                                      gCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAg                                                                                             <  2:3189844/151‑1 (MQ=17)
                                                         tcttcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGc                                                                                          >  1:2063644/1‑151 (MQ=17)
                                                          cttcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCt                                                                                         >  1:3634320/1‑151 (MQ=17)
                                                            tcGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCa                                                                                       <  2:2085706/151‑1 (MQ=17)
                                                                cACGCTCAAATCTGAGTTTTTCTACCTTAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATTGCCTGAGCCCTGTGGCATACAGGGCCCAGGCt                                                                                   >  2:1197105/1‑151 (MQ=255)
                                                                   gCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCa                                                                                <  1:5036929/151‑1 (MQ=17)
                                                                      cAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGGGCCCTGTGGAATGCAGGCCTCAGGCTGCCGCt                                                                             >  2:2622045/1‑151 (MQ=255)
                                                                       aaaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCtt                                                                            <  1:4275642/151‑1 (MQ=14)
                                                                       aaaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCtt                                                                            >  2:807222/1‑151 (MQ=14)
                                                                       aaaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCtt                                                                            <  1:4588606/151‑1 (MQ=14)
                                                                       aaaTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCtt                                                                            <  2:2246173/151‑1 (MQ=17)
                                                                         gtCTGAGGTGTTCTGGCTGGAGTGTGTGGGGCGCATTGATGGGGTGGGGGCGGGTCTGGATGAGTACATTGAGTATTACAACCATGATCGCATCAAGCTGAGGCTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTaa                                                                          <  2:4093589/150‑1 (MQ=11)
                                                                         aTCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTaa                                                                          >  2:4392494/1‑151 (MQ=14)
                                                                            tGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAAc                                                                       <  1:975918/151‑1 (MQ=17)
                                                                             gAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACt                                                                      >  1:184476/1‑151 (MQ=17)
                                                                             gAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACt                                                                      >  1:188151/1‑151 (MQ=17)
                                                                                  tttCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCa                                                                 >  2:872048/1‑151 (MQ=17)
                                                                                        ccTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACtttt                                                           <  2:2749239/151‑1 (MQ=14)
                                                                                         cTGAAGTGTTTCGAGAGCATTGGTGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTg                                                          <  2:496594/151‑1 (MQ=14)
                                                                                          tGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTgg                                                         >  1:3773258/1‑151 (MQ=17)
                                                                                          tGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTgg                                                         >  2:2239061/1‑151 (MQ=14)
                                                                                             aGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTggggg                                                      >  2:1924025/1‑151 (MQ=14)
                                                                                             aGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTggggg                                                      <  1:4462522/151‑1 (MQ=14)
                                                                                                       gagCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATa                                            >  1:1111575/1‑151 (MQ=11)
                                                                                                        agCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATAc                                           >  1:4390386/1‑151 (MQ=14)
                                                                                                           aTTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtg                                        >  2:2568384/1‑151 (MQ=14)
                                                                                                           aTTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtg                                        >  1:3385674/1‑151 (MQ=17)
                                                                                                           aTTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACCTTTGGGGGGCAGTCCATACCtg                                        >  2:3589290/1‑151 (MQ=14)
                                                                                                            ttGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtgt                                       <  1:5097244/151‑1 (MQ=14)
                                                                                                            ttGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCtgt                                       >  1:2927133/1‑151 (MQ=14)
                                                                                                               atgaGCTGAGAGCGGGTCTGAATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTggg                                    <  1:160549/151‑1 (MQ=17)
                                                                                                                tgaGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGa                                   >  2:632416/1‑151 (MQ=17)
                                                                                                                tgaGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGa                                   >  2:3244037/1‑151 (MQ=17)
                                                                                                                  aGCTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGc                                 >  1:2883749/1‑151 (MQ=14)
                                                                                                                    cggagaGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCTCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCgg                               <  2:4566526/149‑1 (MQ=12)
                                                                                                                    cTGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCgg                               <  1:1379842/151‑1 (MQ=17)
                                                                                                                     tGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCggg                              >  1:3291708/1‑151 (MQ=17)
                                                                                                                     tGAGAGCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCggg                              >  1:1759435/1‑151 (MQ=17)
                                                                                                                        gagCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCaa                           <  1:4206794/151‑1 (MQ=18)
                                                                                                                         agCGGGTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAg                          <  2:2063644/151‑1 (MQ=18)
                                                                                                                            gggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGccc                       >  2:3988648/1‑151 (MQ=37)
                                                                                                                             ggTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTTAATGGCCTGAGCCCTGTGGAATATAGGGCTCAGGCTGAAGCTTAAAACTGTTCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAgtctg                      >  2:303181/1‑147 (MQ=11)
                                                                                                                              gTCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCcgc                     <  1:83811/151‑1 (MQ=255)
                                                                                                                               tCTGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCcgcg                    <  2:533689/151‑1 (MQ=37)
                                                                                                                                 tGGATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGaa                  >  1:2265778/1‑151 (MQ=37)
                                                                                                                                  ggATGAGTACCTTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGAAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTACATAACGGTGGGAGCGGGCCAGCCCGCGAAg                 >  2:3575478/1‑151 (MQ=11)
                                                                                                                                  ggATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAg                 <  2:671036/151‑1 (MQ=38)
                                                                                                                                   gATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAgg                >  2:5211630/1‑151 (MQ=25)
                                                                                                                                   gATGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAgg                >  2:4704698/1‑151 (MQ=37)
                                                                                                                                    aTGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAggg               >  1:2341657/1‑151 (MQ=38)
                                                                                                                                    aTGAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAggg               >  1:4997290/1‑151 (MQ=255)
                                                                                                                                      gAGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGcc             >  1:1231131/1‑151 (MQ=255)
                                                                                                                                       aGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCg            >  1:923848/1‑151 (MQ=255)
                                                                                                                                       aGTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCg            >  2:4768610/1‑151 (MQ=255)
                                                                                                                                        gTACATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGc           >  1:3881942/1‑151 (MQ=255)
                                                                                                                                          aCATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAg         >  1:4506073/1‑151 (MQ=255)
                                                                                                                                           cATTCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgc        <  1:1012690/151‑1 (MQ=255)
                                                                                                                                             ttCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgcgc      <  1:263510/151‑1 (MQ=255)
                                                                                                                                             ttCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgcgc      <  2:3385674/151‑1 (MQ=255)
                                                                                                                                             ttCACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAgcgc      <  1:1104159/151‑1 (MQ=255)
                                                                                                                                               cACTATTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCgg    >  1:4636428/1‑151 (MQ=255)
                                                                                                                                                actattACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGGc   <  2:393405/151‑1 (MQ=255)
                                                                                                                                                actattACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGGc   <  1:645540/151‑1 (MQ=255)
                                                                                                                                                 ctattaCAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGGcc  <  2:51257/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TAAAGCAGAGCATGTCTCGAAAAGGCAATTGCCTGGACAATGCGGCCATGGAAAGCTTCTTCGGCACGCTCAAATCTGAGTTTTTCTACCTGAAGTGTTTCGAGAGCATTGATGAGCTGAGAGCGGGTCTGGATGAGTACATTCACTACTACAACCATGATCGCATCAAGCTGAGACTCAATGGCCTGAGCCCTGTGGAATACAGGGCTCAGGCTGCAGCTTAAAACTGTCCAACTTTTGGGGGGCAGTCCATACCTGTGGGAGCGGGCAAGCCCGCGAAGGGCCGCAGCGCGGCC  >  CP009974/1796000‑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: ‑