Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,165,118 Δ1 bp intergenic (‑312/‑163) RPPX_18525 ← / → RPPX_18530 transporter/16S ribosomal RNA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099744,165,1160A.84.0% 76.6 / 35.0 50intergenic (‑310/‑165)RPPX_18525/RPPX_18530transporter/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base . (23/19);  total (28/22)
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

AACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGGCCATGC  >  CP009974/4164966‑4165212
                                                                                                                                                      |                                                                                                
aaCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCa                                                                                                  >  2:4296990/1‑151 (MQ=255)
aaCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCa                                                                                                  >  1:3055713/1‑151 (MQ=255)
 aCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaa                                                                                                 <  2:3467071/151‑1 (MQ=255)
 aCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaa                                                                                                 <  2:3283800/151‑1 (MQ=255)
 aCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaa                                                                                                 <  1:764557/151‑1 (MQ=255)
 aCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaa                                                                                                 >  1:4956311/1‑151 (MQ=255)
    gggTAGGGGAGCAAAGGTAGCCACTGTATTGGGTCAAGTGTTTGAGGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgc                                                                                              <  2:4178266/150‑4 (MQ=14)
     ggTAGCGAAGCAAAGGTAGCCGCGGTATTGAGTCAAGTGTTTGAAGCTAAACAAAAAAAAAAAAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                                                                                             <  2:996101/151‑5 (MQ=14)
     ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                                                                                             >  1:1247991/1‑147 (MQ=17)
     ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                                                                                             >  2:3001703/1‑147 (MQ=255)
     ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                                                                                             >  2:1272166/1‑147 (MQ=17)
        aGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCtt                                                                                         >  1:2878604/1‑151 (MQ=14)
        aGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCCTAACCAACGCTCtt                                                                                         >  2:395646/1‑151 (MQ=14)
         gCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCttt                                                                                        <  1:2910484/151‑1 (MQ=14)
         gCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCttt                                                                                        >  2:4964708/1‑151 (MQ=255)
         gCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCttt                                                                                        <  1:466794/151‑1 (MQ=255)
          cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGTCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                                                                                       >  2:5199298/1‑151 (MQ=14)
          cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                                                                                       <  2:2718078/151‑1 (MQ=14)
          cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                                                                                       >  1:2611951/1‑151 (MQ=14)
          cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCCGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                                                                                       >  2:5199299/1‑151 (MQ=14)
            aaGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAAc                                                                                     >  2:1980312/1‑151 (MQ=14)
            aaGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAAc                                                                                     >  1:2751961/1‑151 (MQ=14)
            aaGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAAc                                                                                     >  1:5373572/1‑151 (MQ=14)
             aGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACa                                                                                    >  2:1656483/1‑151 (MQ=14)
             aGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGAGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACa                                                                                    <  2:1703497/151‑1 (MQ=14)
             aGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGAGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACa                                                                                    <  2:1704970/151‑1 (MQ=14)
              gCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACaa                                                                                   >  2:1431345/1‑151 (MQ=14)
               cAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAACGCTCTTAACCAACGCTCTTTAACaaa                                                                                  >  2:1171278/1‑151 (MQ=255)
                aaaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacg                                                                                               <  1:3157018/138‑3 (MQ=255)
                aaaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacg                                                                                               >  2:3157018/1‑136 (MQ=255)
                 aaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATc                                                                                <  1:2008053/151‑1 (MQ=14)
                   gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:5288161/151‑1 (MQ=11)
                   gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  2:2696719/151‑1 (MQ=11)
                   gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  2:2696310/151‑1 (MQ=14)
                   gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:3948886/151‑1 (MQ=11)
                   gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:1292099/151‑1 (MQ=14)
                     tAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAAt                                                                            <  2:2904996/151‑1 (MQ=11)
                     tAGCCACGGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAAt                                                                            <  2:2336375/151‑1 (MQ=11)
                      aGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATc                                                                           <  1:4241568/151‑1 (MQ=11)
                      aGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATc                                                                           <  2:3999496/151‑1 (MQ=11)
                      aGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCCTTAACAAATCGaacc                                                                           >  1:4522224/1‑149 (MQ=11)
                       gCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGTGTAGAATGCGCGCGACGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCa                                                                          >  2:5211270/1‑151 (MQ=11)
                       gCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCa                                                                          <  1:3910626/151‑1 (MQ=11)
                       gCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGAGTAGAATGCGCGCAACGGTTGAGATGAAAAGCTCTTAACCAACGATCTTTAACAAATCGAATCa                                                                          >  2:5211315/1‑151 (MQ=11)
                        ccACTGTATTGAGTCAAGTGTTTGAAGCTAAAAGAGAAAAAAAAAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCaa                                                                         <  2:820138/151‑1 (MQ=11)
                         cACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAg                                                                        <  1:4696314/151‑1 (MQ=11)
                         cACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAg                                                                        >  1:2338108/1‑151 (MQ=11)
                          aCTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGc                                                                       >  2:1328922/1‑151 (MQ=255)
                           cTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTGTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCa                                                                      >  2:3723316/1‑151 (MQ=11)
                            tGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCaa                                                                     >  1:3431595/1‑151 (MQ=11)
                             gaagcgAGTTAAGTGACTGAAGTTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCaa                                                                     >  1:3462081/6‑151 (MQ=1)
                             gTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAt                                                                    >  2:1340150/1‑151 (MQ=11)
                             gTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAt                                                                    >  2:1120419/1‑151 (MQ=1)
                             gTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAt                                                                    >  2:2506306/1‑151 (MQ=11)
                              tATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAtt                                                                   >  1:1854148/1‑151 (MQ=1)
                              tATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAtt                                                                   <  1:1182131/151‑1 (MQ=255)
                               agagAGTTAAGTGACTGAAGTTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAAtt                                                                   <  1:3859624/148‑1 (MQ=11)
                                 cgAGTTAAGTGACTGAAGTTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCCCTTTAAACAATCGAATCAAGCAATTCg                                                                 >  2:4366049/2‑151 (MQ=11)
                                  gAGTTAAGTGACTGAAGTTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCgt                                                                >  1:4646021/1‑151 (MQ=255)
                                                 agttAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAg                                                 >  2:5007280/4‑151 (MQ=1)
                                                 agttAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAg                                                 >  2:2747741/4‑151 (MQ=1)
                                                                           ttaaaaaTAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTAt                       <  1:1930250/148‑1 (MQ=1)
                                                                                                gccctGACGAAAGCGTAGCATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCAAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGGCCATGc  <  2:4636855/148‑1 (MQ=0)
                                                                                                                                                      |                                                                                                
AACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGGCCATGC  >  CP009974/4164966‑4165212

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: