Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 4,165,100 C→T intergenic (‑294/‑181) 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,1000CT90.1% 176.7 / 12.7 81intergenic (‑294/‑181)RPPX_18525/RPPX_18530transporter/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (37/36);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

ATGCGCACCACTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGGCCATGC  >  CP009974/4164952‑4165212
                                                                                                                                                    |                                                                                                                
aTGCGCACCACTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGa                                                                                                                >  1:4848499/1‑151 (MQ=255)
   ggggcgacTCCAACGGGGTAGCGAAGCAAAGCTAGCCACTGGAGTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaa                                                                                                             <  2:3185266/145‑2 (MQ=17)
   cgcACCACTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGCGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaa                                                                                                             <  2:5253371/151‑2 (MQ=255)
   cgcACCACTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaa                                                                                                             >  1:140019/1‑150 (MQ=255)
    gcACCACTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaag                                                                                                            <  1:1863242/151‑3 (MQ=255)
    gcACCACTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaag                                                                                                            >  1:779713/1‑149 (MQ=255)
    gcACCACTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaag                                                                                                            <  2:2791774/151‑3 (MQ=255)
     caccagTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTGTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaagc                                                                                                           <  2:3522327/151‑4 (MQ=255)
     caccacTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaagc                                                                                                           <  1:2002608/151‑4 (MQ=255)
       ccacTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaagctc                                                                                                         >  1:1410620/1‑146 (MQ=255)
        cacTTCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaagctct                                                                                                        >  2:4474252/1‑145 (MQ=255)
        cacTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaagctct                                                                                                        <  2:661067/151‑7 (MQ=255)
        cacTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGaaaagctct                                                                                                        <  1:3203593/151‑7 (MQ=255)
          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                                                                                                  >  1:3055713/1‑151 (MQ=255)
              aaCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCa                                                                                                  >  2:4296990/1‑151 (MQ=255)
               aCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaa                                                                                                 <  2:3283800/151‑1 (MQ=255)
               aCGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaa                                                                                                 <  2:3467071/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:3001703/1‑147 (MQ=255)
                   ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                                                                                             >  1:1247991/1‑147 (MQ=17)
                   ggTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCaacgct                                                                                             >  2:1272166/1‑147 (MQ=17)
                      aGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCtt                                                                                         >  1:2878604/1‑151 (MQ=14)
                      aGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCCTAACCAACGCTCtt                                                                                         >  2:395646/1‑151 (MQ=14)
                       gCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCttt                                                                                        <  1:466794/151‑1 (MQ=255)
                       gCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCttt                                                                                        >  2:4964708/1‑151 (MQ=255)
                       gCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCttt                                                                                        <  1:2910484/151‑1 (MQ=14)
                        cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGTCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                                                                                       >  2:5199298/1‑151 (MQ=14)
                        cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                                                                                       >  1:2611951/1‑151 (MQ=14)
                        cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                                                                                       <  2:2718078/151‑1 (MQ=14)
                        cGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCCGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTa                                                                                       >  2:5199299/1‑151 (MQ=14)
                          aaGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAAc                                                                                     >  1:2751961/1‑151 (MQ=14)
                          aaGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAAc                                                                                     >  1:5373572/1‑151 (MQ=14)
                          aaGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAAc                                                                                     >  2:1980312/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                                                                              <  2:2696310/151‑1 (MQ=14)
                                 gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  2:2696719/151‑1 (MQ=11)
                                 gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:1292099/151‑1 (MQ=14)
                                 gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:3948886/151‑1 (MQ=11)
                                 gCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGa                                                                              <  1:5288161/151‑1 (MQ=11)
                                   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: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:1182131/151‑1 (MQ=255)
                                            tATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGATGAAAAGCTCTTAACCAACGCTCTTTAACAAATCGAATCAAGCAAtt                                                                   >  1:1854148/1‑151 (MQ=1)
                                             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)
                                                                                                                                                    |                                                                                                                
ATGCGCACCACTCCAACGCGGTAGCGAAGCAAAGCTAGCCACTGTATTGAGTCAAGTGTTTGAAGCTAAACGAGTTTCTCGCAAAAAACTTCAAAATAAACGCTTGACAGGCTCTGAGGAAAGCGTAGAATGCGCGCCTCGGTTGAGACGAAAGGTTCTTAACCAAACGCTCTTTAACAAATCGAATCAAGCAATTCGTGTGGGTGCTTGTGAGTACGGACTGATAGTCAAAAAGATTATCAGCATCACAAGTGGCCATGC  >  CP009974/4164952‑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: