Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,165,107 T→C intergenic (‑301/‑174) 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,1070TC88.2% 127.5 / 12.7 68intergenic (‑301/‑174)RPPX_18525/RPPX_18530transporter/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (5/3);  new base C (32/28);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 7.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGGCCATGC  >  CP009974/4164962‑4165212
                                                                                                                                                 |                                                                                                         
cTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTa                                                                                                      >  2:390869/1‑151 (MQ=255)
cTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTa                                                                                                      <  2:1359570/151‑1 (MQ=255)
cTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTa                                                                                                      <  2:1359252/151‑1 (MQ=255)
  ccAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAAc                                                                                                    <  1:354048/151‑1 (MQ=255)
  ccAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGGTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAAc                                                                                                    >  2:187869/1‑151 (MQ=255)
    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                                                                                             >  2:1272166/1‑147 (MQ=17)
         ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                                                                                             >  1:1247991/1‑147 (MQ=17)
         ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                                                                                             >  2:3001703/1‑147 (MQ=255)
            aGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCtt                                                                                         >  1:2878604/1‑151 (MQ=14)
            aGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCCTAACCAACGCTCtt                                                                                         >  2:395646/1‑151 (MQ=14)
             gCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCttt                                                                                        >  2:4964708/1‑151 (MQ=255)
             gCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCttt                                                                                        <  1:2910484/151‑1 (MQ=14)
             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:5373572/1‑151 (MQ=14)
                aaGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAAc                                                                                     >  1:2751961/1‑151 (MQ=14)
                 aGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACa                                                                                    >  2:1656483/1‑151 (MQ=14)
                 aGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGAGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACa                                                                                    <  2:1704970/151‑1 (MQ=14)
                 aGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGAGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACa                                                                                    <  2:1703497/151‑1 (MQ=14)
                  gCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACaa                                                                                   >  2:1431345/1‑151 (MQ=14)
                   cAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAACGCTCTTAACCAACGCTCTTTAACaaa                                                                                  >  2:1171278/1‑151 (MQ=255)
                    aaaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacg                                                                                               >  2:3157018/1‑136 (MQ=255)
                    aaaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacg                                                                                               <  1:3157018/138‑3 (MQ=255)
                     aaGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATc                                                                                <  1:2008053/151‑1 (MQ=14)
                       gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:1292099/151‑1 (MQ=14)
                       gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:5288161/151‑1 (MQ=11)
                       gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  2:2696719/151‑1 (MQ=11)
                       gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:3948886/151‑1 (MQ=11)
                       gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  2:2696310/151‑1 (MQ=14)
                         tAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAAt                                                                            <  2:2904996/151‑1 (MQ=11)
                         tAGCCACGGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAAt                                                                            <  2:2336375/151‑1 (MQ=11)
                          aGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATc                                                                           <  2:3999496/151‑1 (MQ=11)
                          aGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATc                                                                           <  1:4241568/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:2338108/1‑151 (MQ=11)
                             cACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAg                                                                        <  1:4696314/151‑1 (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:2506306/1‑151 (MQ=11)
                                 gTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAt                                                                    >  2:1340150/1‑151 (MQ=11)
                                 gTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAt                                                                    >  2:1120419/1‑151 (MQ=1)
                                  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:2747741/4‑151 (MQ=1)
                                                     agttAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAg                                                 >  2:5007280/4‑151 (MQ=1)
                                                                               ttaaaaaTAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTAt                       <  1:1930250/148‑1 (MQ=1)
                                                                                                    gccctGACGAAAGCGTAGCATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCAAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGGCCATGc  <  2:4636855/148‑1 (MQ=0)
                                                                                                                                                 |                                                                                                         
CTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGGCCATGC  >  CP009974/4164962‑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: