Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,475,329 T→G 5.6% intergenic (‑85/+72) tal ← / ← dusA transaldolase/tRNA dihydrouridine synthase DusA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,475,3290TG 5.6% 195.4 / 3.8 73intergenic (‑85/+72)tal/dusAtransaldolase/tRNA dihydrouridine synthase DusA
Reads supporting (aligned to +/- strand):  ref base T (40/28);  new base G (2/2);  total (42/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

GTCCCCGGTGTCGGCGACCACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATGTCTGCCGACAG  >  NC_002947/2475185‑2475478
                                                                                                                                                |                                                                                                                                                      
ggccccGGGGGCGGCGACCACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGccc                                                                                                                                                   <  1:665351/148‑1 (MQ=255)
agacattcacaaggcggctgaacccttggcggtgtttgaccaggccgtggagttgttgcaggggcggtgacctcagagggagggaacctcgggcttgcttgcgactcgaatcttcacctttcgtacagccctgTGCGGACCCGCTGGccc                                                                                                                                                   <  2:1229168‑M1/17‑1 (MQ=255)
     cGGTGTCGGCGACCACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGcaca                                                                                                                                              <  2:757397/150‑1 (MQ=255)
             gCGACCACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTa                                                                                                                                      >  1:461115/1‑150 (MQ=255)
             gCGACCACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTa                                                                                                                                      >  1:1116989/1‑150 (MQ=255)
             gCGACCACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTa                                                                                                                                      >  1:461105/1‑150 (MQ=255)
               gACCACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCGGTACg                                                                                                                                    >  1:756174/1‑150 (MQ=255)
                 ccACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGaa                                                                                                                                  >  2:1262989/1‑150 (MQ=255)
                  cACGGTGGTGAACTGCTTGAGTTGTTCCCGCTTGGAGGTCATGGGTGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGaaa                                                                                                                                 <  1:547838/150‑1 (MQ=255)
                  cACGGTGGTGAACTGCTTGAGTTGTCCCAGCTTCGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACACATCGCTTCCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGaaa                                                                                                                                 <  1:125153/150‑1 (MQ=255)
                   aCGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAg                                                                                                                                >  1:147325/1‑150 (MQ=255)
                     ggtggtGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGt                                                                                                                              <  1:1639485/150‑1 (MQ=255)
                     ggtggtGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGt                                                                                                                              <  1:1639606/150‑1 (MQ=255)
                     ggtggtGAACTGCTTGAGATGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTAAAGGGCCGCCGAAACGCTGCGGGCCTGCTTGCCCGCACAGGTCTGTACGAAAGGt                                                                                                                              >  2:1753720/1‑150 (MQ=255)
                      gtggtgAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTg                                                                                                                             <  2:1515210/150‑1 (MQ=255)
                       tggtgAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGa                                                                                                                            <  1:905181/150‑1 (MQ=255)
                        ggtgAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGaa                                                                                                                           <  2:1253372/150‑1 (MQ=255)
                              cTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCg                                                                                                                     >  1:569404/1‑150 (MQ=255)
                                  ttGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTGCCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTc                                                                                                                 <  2:824717/150‑1 (MQ=255)
                                  ttGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCACCCGCTCAAGGGCCGCCGAAACCCTGCGGACACGCTGGCCCGCACAAGGCCGTACGAAAGGTGAAGATTCGAGTc                                                                                                                 >  1:1866317/1‑150 (MQ=255)
                                       ttgttCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGGCCAGCGGGTCCGcagg                                                                                                                                              >  1:243616/1‑114 (MQ=255)
                                          ttCCAGCTTGGAGGTCATGGGCGTGCTCTGGCATGTGCATTTACTCGGCATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGCCgcaagcaa                                                                                                        >  1:365917/1‑150 (MQ=255)
                                             cAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGcc                                                                                                     >  2:278158/1‑150 (MQ=255)
                                             cAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGcc                                                                                                     >  1:1189623/1‑150 (MQ=255)
                                             cAGCTTGGAGGTCATGCGCTTGCTCTGTCCTGTGCATTTACTCGCCATTACCCGAGCCCCGGCAGCCCCTCAAGCGCCGCCGCCACCCTGCGGCCCCCCTGCCCCCCCCAGGGCTCTACGACCGGTGCACATTCGAGTCtcacgccaccc                                                                                                      >  2:839127/1‑145 (MQ=255)
                                               gCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTa                                                                                                                                      >  1:99291/1‑116 (MQ=255)
                                               gCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTa                                                                                                                                      <  2:99291/116‑1 (MQ=255)
                                                 ttGGAGGTCATGGGCGTGCTCTGTCCTGTGGATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCCGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAg                                                                                                 >  2:1119624/1‑150 (MQ=255)
                                                 ttGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAg                                                                                                 >  1:1049389/1‑150 (MQ=255)
                                                    gAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGtt                                                                                              <  1:632266/150‑1 (MQ=255)
                                                    gAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGtt                                                                                              >  2:604281/1‑150 (MQ=255)
                                                     aGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGCACCGCCGTCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTtc                                                                                             <  2:259748/150‑1 (MQ=255)
                                                          aTGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTtccctcc                                                                                        >  1:57007/1‑150 (MQ=255)
                                                          aTGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTtccctcc                                                                                        >  1:271523/1‑150 (MQ=255)
                                                          aTGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTtccc                                                                                           <  1:1024399/147‑1 (MQ=255)
                                                          aTGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTtccc                                                                                           >  2:1024399/1‑147 (MQ=255)
                                                            gggCGTGCCCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTccctccct                                                                                      <  2:1042000/150‑1 (MQ=255)
                                                                gTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGa                                                                                  <  2:147325/150‑1 (MQ=255)
                                                                gTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGa                                                                                  <  1:1262989/150‑1 (MQ=255)
                                                                gTGCTATGTCCTGTGCATTTACTTGACATTATCCGAGACCTGGCAGCCGCGCAAGGACCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAGAGGGGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGa                                                                                  >  2:420245/1‑150 (MQ=255)
                                                                    tctGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTc                                                                              >  1:979763/1‑150 (MQ=255)
                                                                                aTTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGca                                                                  >  2:1489048/1‑150 (MQ=255)
                                                                                aTTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGca                                                                  >  1:1815262/1‑150 (MQ=255)
                                                                                aTTTACTCGACATTACCCGAGCCCAGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGca                                                                  >  1:258135/1‑150 (MQ=255)
                                                                                      tCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCaacaacc                                                            >  2:1765056/1‑149 (MQ=255)
                                                                                         aCATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACACCTCCa                                                         >  2:22931/1‑150 (MQ=255)
                                                                                                gttgctgccGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACGGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTCCAACCACTCCACGGCCTg                                                  >  2:769690/8‑150 (MQ=255)
                                                                                                ccGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGTTTAAGATTCGAGTCGCAAGAAACCCCGACTTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTg                                                  >  1:1869686/1‑150 (MQ=255)
                                                                                                ccGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTg                                                  >  1:250355/1‑150 (MQ=255)
                                                                                                ccGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCCCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTg                                                  >  2:1363988/1‑150 (MQ=255)
                                                                                                ccGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACACTGCGGACCCGCTGGCCCGCACAGGGCTGTAAGAAAGGGGAAGATTCGCGACGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTg                                                  >  1:134037/1‑150 (MQ=255)
                                                                                                     ccccGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCaa                                             >  2:1339559/1‑150 (MQ=255)
                                                                                                               cttcacctttcgtacagccctgTGCGGGCCAGCGGGTCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCaa                                   <  1:1229168‑M1/128‑1 (MQ=255)
                                                                                                               cGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCTAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCaa                                   <  1:278158/150‑1 (MQ=255)
                                                                                                                    aaGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAg                                                                                                 <  1:158977/83‑1 (MQ=255)
                                                                                                                    aaGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAg                                                                                                 >  2:158977/1‑83 (MQ=255)
                                                                                                                      gggCCGCCGAAACCCTCCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCa                            <  1:1602968/150‑1 (MQ=255)
                                                                                                                        gccgccGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGCTTCGAGGCGCAAGCAAGCACCAGGGTCCCTCCCTCTGAGGTCAGCGCCTCTGCAACAACTACACGGCCTGGGCAAACCCCGCCAAGGGTTCAgc                          >  2:893776/1‑150 (MQ=255)
                                                                                                                        gccgccGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAgg                                                                                                <  1:1356489/80‑1 (MQ=255)
                                                                                                                        gccgccGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAgg                                                                                                >  2:1356489/1‑80 (MQ=255)
                                                                                                                                ccctgTGCGGGCCAGCGGGTCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGa                                                                                  <  2:283259‑M1/81‑1 (MQ=255)
                                                                                                                                ccctgTGCGGGCCAGCGGGTCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGa                                                                                  >  1:283259‑M1/6‑86 (MQ=255)
                                                                                                                                aaCCCTGCCGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTtgt                  <  1:893776/150‑1 (MQ=255)
                                                                                                                                aaCCATGCGGCCCCGCTGGCACGCACAGGGCTGCACTAACGGTGTAGACTCGAGCCGAAAGTATGCCCTTTGTTCCCTCCCTCTTAGGTCACCGCCCCTGCAACACCCCCACGGCCTGGTCAAACTCCGCCAAGGGTTCAGCCGCCTtgt                  <  1:894876/150‑1 (MQ=255)
                                                                                                                                   ccTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCAACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGaa               >  1:1337955/1‑150 (MQ=255)
                                                                                                                                    cTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGtt                                                                                              <  1:921831/70‑1 (MQ=255)
                                                                                                                                    cTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGtt                                                                                              >  2:921831/1‑70 (MQ=255)
                                                                                                                                     TGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATg             <  2:1889378‑M1/150‑1 (MQ=255)
                                                                                                                                          aCCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCACGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCTCCGCCCCTGCAACAACTCCACCGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATGTCTGc        <  2:979763/150‑1 (MQ=255)
                                                                                                                                          aCCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCCCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATGTCTGc        <  1:1363988/150‑1 (MQ=255)
                                                                                                                                           cccGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATGTCTGcc       <  2:271523/150‑1 (MQ=255)
                                                                                                                                            ccGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGCCTAGAACAGACGTACCGGTCCAGGATTTCTCCGTTGCATTTCACCGCCCCTGCAACAACTCCACGGCCTGTTCAAACACCGCCAAGGTGGCAGCCGCCTTGTGAATGTCTGCCg      >  1:35751/1‑150 (MQ=255)
                                                                                                                                             cGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATGTCTGCCGa     >  2:1738849/1‑150 (MQ=255)
                                                                                                                                               atGGCCCGAAAACGGTTTTACTAAATGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATGTCTGCCGACa   <  2:1049389/149‑1 (MQ=255)
                                                                                                                                                tGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCCCCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATGTCTGCCGACAg  <  1:604281/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GTCCCCGGTGTCGGCGACCACGGTGGTGAACTGCTTGAGTTGTTCCAGCTTGGAGGTCATGGGCGTGCTCTGTCCTGTGCATTTACTCGACATTACCCGAGCCCCGGCAGCCGCTCAAGGGCCGCCGAAACCCTGCGGACCCGCTGGCCCGCACAGGGCTGTACGAAAGGTGAAGATTCGAGTCGCAAGCAAGCCCGAGGTTCCCTCCCTCTGAGGTCACCGCCCCTGCAACAACTCCACGGCCTGGTCAAACACCGCCAAGGGTTCAGCCGCCTTGTGAATGTCTGCCGACAG  >  NC_002947/2475185‑2475478

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: