Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,039,236 C→T 13.8% intergenic (‑256/+38) PP_5290 ← / ← pyrE transposase/orotate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,039,2360CT13.8% 139.9 / 12.7 58intergenic (‑256/+38)PP_5290/pyrEtransposase/orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (27/23);  new base T (4/4);  total (31/27)
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.94e-01

GAAAGGAGGAGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCT  >  NC_002947/6039089‑6039378
                                                                                                                                                   |                                                                                                                                              
gAAAGGAGGAGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCct                                                                                                                                               >  2:56909/1‑149 (MQ=255)
  aaGGAGGAGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCctct                                                                                                                                             <  1:51927/149‑1 (MQ=255)
  aaGGAGGAGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCctct                                                                                                                                             <  4:45372/149‑1 (MQ=255)
     gaggaggGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTg                                                                                                                                            >  7:38990/1‑147 (MQ=255)
        gaggGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTgaagaa                                                                                                                                       <  6:63244/149‑1 (MQ=255)
             gTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCg                                                                                                                                  >  4:251401/1‑149 (MQ=255)
             gTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCg                                                                                                                                  >  2:286051/1‑149 (MQ=255)
             gTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCg                                                                                                                                  >  2:131735/1‑149 (MQ=255)
              tGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCGc                                                                                                                                 <  3:251401/149‑1 (MQ=255)
                       cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttttttt                                                                                                                         <  6:297558/148‑3 (MQ=255)
                       cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttttt                                                                                                                           <  5:277631/146‑1 (MQ=255)
                          cGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTc                                                                                                                      <  1:286051/148‑1 (MQ=255)
                                      cgcgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATc                                                                                                         <  3:114758/149‑1 (MQ=255)
                                       gcgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGAt                                                                                                          >  1:48743/1‑147 (MQ=255)
                                                     gCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAgg                                                                                                            <  8:223979/131‑1 (MQ=255)
                                                     gCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAgg                                                                                                            >  7:223979/1‑131 (MQ=255)
                                                     gCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAgg                                                                                                            >  7:228203/1‑131 (MQ=255)
                                                     gCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAgg                                                                                                            <  8:228203/131‑1 (MQ=255)
                                                     gCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACGgcg                                                                                          >  2:273516/1‑149 (MQ=255)
                                                         gTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCg                                                                                      <  2:83534/149‑1 (MQ=255)
                                                             tGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCAGGTAg                                                                                  >  8:176078/1‑149 (MQ=255)
                                                             tGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTACTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAg                                                                                  >  8:180739/1‑149 (MQ=255)
                                                               gCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGCAGGTAGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGc                                                                                >  5:224728/1‑149 (MQ=255)
                                                                 cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                              >  2:56431/1‑149 (MQ=255)
                                                                 cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                              >  8:84690/1‑149 (MQ=255)
                                                                 cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTACTGGGACAGCATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                              >  6:168728/1‑149 (MQ=255)
                                                                  ctcctcctcCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTc                                                                             <  8:192032/149‑1 (MQ=255)
                                                                         tccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGcc                                                                      >  2:86794/1‑149 (MQ=255)
                                                                         tccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGcc                                                                      >  8:153858/1‑149 (MQ=255)
                                                                         tccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGcc                                                                      >  7:141248/1‑149 (MQ=255)
                                                                           ccTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGcc                                                                      <  1:131735/147‑1 (MQ=255)
                                                                            cTTCAAGTCCAACCATACAAGCGGGCTAACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAACGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGc                                                                   <  1:302552/149‑1 (MQ=255)
                                                                                  gTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACGAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGa                                                                           >  6:291120/1‑135 (MQ=255)
                                                                                     cTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTg                                                           <  1:56909/148‑1 (MQ=255)
                                                                                     cTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGt                                                          <  2:48743/149‑1 (MQ=255)
                                                                                       aCCATACAAGCGGGTTTAGCCGCGAAAGGGCCAGATCAGGCGGCACAAAAACAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                        <  6:198706/149‑1 (MQ=255)
                                                                                        tcATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTc                                                       <  7:127403/148‑1 (MQ=255)
                                                                                        ccATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAgg                                                                                 <  2:293664/123‑1 (MQ=255)
                                                                                        ccATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAgg                                                                                 >  1:293664/1‑123 (MQ=255)
                                                                                        ccATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                        >  8:225266/1‑148 (MQ=255)
                                                                                        ccATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTc                                                       <  1:86794/149‑1 (MQ=255)
                                                                                         cATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCa                                                      <  8:38990/149‑1 (MQ=255)
                                                                                          aTACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTc                                                       >  8:127403/1‑147 (MQ=255)
                                                                                           tACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGc                                                    >  5:295889/1‑149 (MQ=255)
                                                                                              aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                 <  7:180739/149‑1 (MQ=255)
                                                                                              aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                 <  7:176078/149‑1 (MQ=255)
                                                                                                      ttACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                 >  6:94456/1‑141 (MQ=255)
                                                                                                      ttACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                 <  5:94456/141‑1 (MQ=255)
                                                                                                          ccGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                 >  1:69029/1‑137 (MQ=255)
                                                                                                          ccGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgt                                                 <  2:69029/137‑1 (MQ=255)
                                                                                                                    gCCAGATCAGGCTGCAGAAAAAAAGGCGACCCTCTGAAGAAGGTCGGATTTTTCATTCCTGGTCCAGGATCAGATCCCGTACTGCGCCCGGTAGGACTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCa                           <  5:168728/149‑1 (MQ=255)
                                                                                                                    gCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCa                           <  1:56431/149‑1 (MQ=255)
                                                                                                                        gATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGTCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCAc                       >  7:120152/1‑149 (MQ=255)
                                                                                                                             ggcggcACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTggg                  >  8:205962/1‑149 (MQ=255)
                                                                                                                                   acaAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGc             >  4:78115/1‑148 (MQ=255)
                                                                                                                                    atgaaaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTgc                                                                                           <  7:267743/66‑1 (MQ=18)
                                                                                                                                    atgaaaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTgc                                                                                           >  8:267743/4‑69 (MQ=18)
                                                                                                                                     aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGa           >  5:255641/1‑148 (MQ=255)
                                                                                                                                       aaaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaac         >  1:207410/1‑147 (MQ=255)
                                                                                                                                         aaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACaa      <  7:225266/149‑1 (MQ=255)
                                                                                                                                             gCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt  >  8:197992/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
GAAAGGAGGAGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCT  >  NC_002947/6039089‑6039378

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: