Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 1,920,803 T→A intergenic (‑432/+123) RPPX_08560 ← / ← RPPX_08565 16S ribosomal RNA/CoA‑transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099741,920,8030TA97.4% 239.6 / ‑4.1 114intergenic (‑432/+123)RPPX_08560/RPPX_0856516S ribosomal RNA/CoA‑transferase
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base A (54/56);  total (58/56)
Fisher's exact test for biased strand distribution p-value = 2.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 114 total aligned reads displayed.

CTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTTGATCGTAGCCCTTACAGC‑‑A‑‑CCTTCCGAATCGTCTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAGG  >  CP009974/1920691‑1920950
                                                                                                                    |                                                                                                                                                   
cTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgag                                                                                                               >  1:3927399/1‑151 (MQ=2)
 ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGCTGTCGCTGGAAGTGGGGCGCACTATAATGGGATTCgaga                                                                                                              >  2:1465645/1‑151 (MQ=2)
 ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga                                                                                                              <  2:2855223/151‑1 (MQ=2)
 ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga                                                                                                              <  2:6126431/151‑1 (MQ=2)
 ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga                                                                                                              <  1:4873599/151‑1 (MQ=2)
 ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga                                                                                                              <  1:5821978/151‑1 (MQ=2)
 ttGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCgaga                                                                                                              <  2:1913893/151‑1 (MQ=2)
   ggAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCTTTCCGATGTCGCTGGAAGTGGGGCGCATTATACGGCGCTTCGAGAgg                                                                                                            >  2:232995/1‑151 (MQ=2)
    gAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAggg                                                                                                           >  1:1748310/1‑151 (MQ=2)
    gAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAggg                                                                                                           <  1:6005577/151‑1 (MQ=2)
     aaGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGc                                                                                                          >  1:638951/1‑151 (MQ=2)
     aaGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGc                                                                                                          >  1:2643515/1‑151 (MQ=2)
       gcgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGt                                                                                                        <  2:2918939/151‑1 (MQ=11)
       gcgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGt                                                                                                        <  2:4032767/151‑1 (MQ=2)
        cgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTc                                                                                                       <  2:4428118/151‑1 (MQ=2)
        cgcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATGATAAGGGGATTCGAGAGGGCGTc                                                                                                       >  2:2514959/1‑151 (MQ=2)
         gcTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCa                                                                                                      <  2:1201368/151‑1 (MQ=2)
          cTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCTATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa                                                                                                     >  2:1739571/1‑151 (MQ=2)
          cTGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa                                                                                                     >  1:611301/1‑151 (MQ=2)
           tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa                                                                                                     >  2:2521693/1‑150 (MQ=2)
           tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCaa                                                                                                     <  1:2521693/150‑1 (MQ=2)
           tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAAc                                                                                                    >  1:4154449/1‑151 (MQ=2)
           tGTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAAc                                                                                                    <  1:5317769/151‑1 (MQ=2)
            gTCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAAcc                                                                                                   <  2:638550/151‑1 (MQ=11)
             tCAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCt                                                                                                  >  1:173209/1‑151 (MQ=2)
              cAACCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCtt                                                                                                 <  2:1673337/151‑1 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                                <  2:425745/151‑1 (MQ=2)
               aaCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCttt                                                                                                >  1:5513620/1‑151 (MQ=2)
                aCCATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTa                                                                                               >  1:1515934/1‑151 (MQ=2)
                 ccATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTaa                                                                                              <  1:3548520/151‑1 (MQ=2)
                 ccATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTaa                                                                                              <  1:81832/151‑1 (MQ=2)
                 ccATCATCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTaa                                                                                              <  2:2910332/151‑1 (MQ=2)
                  catcatCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAt                                                                                             >  2:684853/1‑151 (MQ=2)
                  catcatCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAt                                                                                             >  1:3280911/1‑151 (MQ=2)
                  catcatCTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAt                                                                                             >  2:3542561/1‑151 (MQ=2)
                   atcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAtt                                                                                            <  1:406450/151‑1 (MQ=2)
                   atcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAtt                                                                                            <  2:2715844/151‑1 (MQ=2)
                    tcatcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAAttt                                                                                           <  2:5443525/151‑1 (MQ=2)
                      atcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCa                                                                                         <  2:613995/151‑1 (MQ=2)
                      atcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCa                                                                                         <  1:5552615/151‑1 (MQ=2)
                      atcTCAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGC‑‑TTTCCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCa                                                                                         <  2:3325031/151‑1 (MQ=2)
                          cAACCGCTTTCGATCATTCGACCGTAGCCCCTCCAGCTTT‑‑CCTTTCACTACTTAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATaa                                                                                     <  2:905783/151‑1 (MQ=11)
                                                                            actacttAACTCATTGAATCTCAAGGAGTTTTTCTTTTCTTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTc                                       >  2:3716136/7‑151 (MQ=21)
                                                                            actacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTc                                       >  1:5476146/7‑151 (MQ=255)
                                                                            actacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTc                                       >  1:1907592/7‑151 (MQ=255)
                                                                            actacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTc                                       <  2:549432/145‑1 (MQ=255)
                                                                              ttcttAACTCATTGAATCTCAAGGAGTTTTTCGCTCTGATGCCGCTGGAAGTGGGCCGTATTTTAAGGGGATTTTAGAGCGCGTCAACCTTTAATTTCAATAACCGATGATATCGCAGAAAAAACAAAGCGGGGACGCCTGACTGcccccc                                     >  2:735182/5‑151 (MQ=18)
                                                                              tacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTccc                                     >  1:5431657/5‑151 (MQ=255)
                                                                              tacttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTccc                                     >  2:1468401/5‑151 (MQ=255)
                                                                               tcttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcccc                                    <  1:6413016/148‑1 (MQ=255)
                                                                               acttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcccc                                    <  1:2389049/148‑1 (MQ=255)
                                                                               acttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTcccc                                    <  1:4772594/148‑1 (MQ=255)
                                                                                cttAACTCATTGAATCTCAAGGAGTTTTTCGTTTCGATTTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCg                                   >  2:5100858/3‑151 (MQ=255)
                                                                                cttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCg                                   >  1:336414/3‑151 (MQ=255)
                                                                                 ttAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGc                                  >  2:2734525/2‑151 (MQ=255)
                                                                                  tAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCt                                 <  1:1065415/151‑1 (MQ=255)
                                                                                  tAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCt                                 >  2:605087/1‑151 (MQ=255)
                                                                                  tAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCt                                 >  1:4472161/1‑151 (MQ=255)
                                                                                  tAACTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCt                                 >  1:5420264/1‑151 (MQ=255)
                                                                                   aaCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTa                                >  1:2899026/1‑151 (MQ=255)
                                                                                    aCTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTAt                               >  1:4062873/1‑151 (MQ=255)
                                                                                     cTCCTTGAATCTCAAGGAGATTTTCGTTCCGAGGTCGGTGGAAGTGGGGCGCATTTTAAGGGGGTTGGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTAtt                              >  2:4141782/1‑151 (MQ=25)
                                                                                     cTCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTAtt                              <  1:5303340/151‑1 (MQ=255)
                                                                                      tCATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTc                             <  2:2423018/151‑1 (MQ=255)
                                                                                       cATTGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCa                            <  1:2100678/151‑1 (MQ=255)
                                                                                         ttGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATc                          <  1:5568355/151‑1 (MQ=255)
                                                                                         ttGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATc                          >  2:3918924/1‑151 (MQ=255)
                                                                                          tGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCt                         <  2:1095158/151‑1 (MQ=255)
                                                                                          tGAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGAGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCt                         <  1:5591947/151‑1 (MQ=255)
                                                                                           gAATCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTa                        >  1:4379242/1‑151 (MQ=255)
                                                                                            aaTCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTag                       >  1:1939930/1‑151 (MQ=255)
                                                                                            aaTCTCAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTag                       >  2:5572217/1‑151 (MQ=255)
                                                                                             aTCTCAAGGAGTTTTTTTTTTTTTTTTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGTGTGGCCTGACGGCCTCCCCGCTATTCATCTaga                      >  2:2345114/1‑151 (MQ=25)
                                                                                              tctcAAGGAGTTTTTTGTTCTTTTTTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTagag                     >  2:1496180/1‑151 (MQ=25)
                                                                                              tctcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTagag                     <  2:962713/151‑1 (MQ=255)
                                                                                               ctcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGc                    >  1:226244/1‑151 (MQ=255)
                                                                                               ctcAAGGAGCTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGTGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGAGGAGGCCTTACGGCCACCCCGCTATTCACCTAGAGc                    >  1:2806457/1‑151 (MQ=21)
                                                                                                tcAAGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGca                   <  1:4127146/151‑1 (MQ=255)
                                                                                                  aaGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATGGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGcaca                 >  2:6219545/1‑151 (MQ=255)
                                                                                                  aaGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGcaca                 >  2:2199350/1‑151 (MQ=255)
                                                                                                   aGGAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCacac                <  1:4727509/151‑1 (MQ=255)
                                                                                                    ggAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACAcc               <  1:1338081/151‑1 (MQ=255)
                                                                                                     gAGTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGCGCACACCa              >  2:698284/1‑151 (MQ=255)
                                                                                                     gAGTTTTTCGGTCCGATGTCGCTGGAAGTGGGGCGCAGTATAAGGTGATTCGAGAGGGCGCCAACCTTTAATTTCAATAAACGATGATAGCGCAGTAAAAACAACGCTGAGAGGCCTGACGGCCTCCCCGCGATCCCTCTAGGGCCCAccc              >  2:3822304/1‑150 (MQ=18)
                                                                                                       gTTTTTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGTTTCGAGAGGGCGTCAACCTTTAATCTCAATAACCGTTGTTAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGACCTCCCCGCTATTCCGCTAGAGCACACCAGc            >  2:2456359/1‑151 (MQ=25)
                                                                                                        ttttTCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGca           <  1:3323965/151‑1 (MQ=255)
                                                                                                         ttttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcac          <  1:667322/148‑1 (MQ=255)
                                                                                                          tttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcacc         <  2:395217/149‑1 (MQ=255)
                                                                                                          tttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGAAAGAGAAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcacc         <  2:3711316/149‑1 (MQ=255)
                                                                                                           ttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcacca        <  1:6219545/150‑1 (MQ=255)
                                                                                                           ttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcacca        <  1:5125674/150‑1 (MQ=255)
                                                                                                           ttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcacca        >  1:4402353/2‑151 (MQ=255)
                                                                                                           ttCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcacca        >  2:727181/2‑151 (MQ=255)
                                                                                                            tCGTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGcaccac       <  1:3223256/151‑1 (MQ=255)
                                                                                                              gTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGa     <  1:3430898/151‑1 (MQ=255)
                                                                                                              gTTCCGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGa     >  2:3430898/1‑151 (MQ=255)
                                                                                                                 ccGATGTGGGTGTAAGTGGGGCGCGTTATAAGGGGATCCGAGAGGGCGTCAACCTTTACTTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAgg  <  2:5741072/151‑1 (MQ=39)
                                                                                                                 ccGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAgg  >  2:461764/1‑151 (MQ=255)
                                                                                                                 ccGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAgg  >  2:2656409/1‑151 (MQ=255)
                                                                                                                 ccGATGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAgg  >  1:1559919/1‑151 (MQ=255)
                                                                                                                    |                                                                                                                                                   
CTTGGAAGCGCTGTCAACCATCATCTCAACCGCTTTCGATCATTTGATCGTAGCCCTTACAGC‑‑A‑‑CCTTCCGAATCGTCTAACTCATTGAATCTCAAGGAGTTTGTCGTTCCGTTGTCGCTGGAAGTGGGGCGCATTATAAGGGGATTCGAGAGGGCGTCAACCTTTAATTTCAATAACCGTTGATAGCGCAGAAAAAACAAAGCGGGGAGGCCTGACGGCCTCCCCGCTATTCATCTAGAGCACACCAGCACCACGAAGG  >  CP009974/1920691‑1920950

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: