Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,039,241 A→C 8.1% intergenic (‑261/+33) 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,2410AC 8.1% 181.0 / 5.3 59intergenic (‑261/+33)PP_5290/pyrEtransposase/orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base C (3/2);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGGAGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGG  >  NC_002947/6039094‑6039387
                                                                                                                                                   |                                                                                                                                                  
gaggaggGGTGATTTCTCCTACGGTATGTACTGCGCGAACAGTCAGAAGCGAGTTAAGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTgaa                                                                                                                                                   <  5:86572/149‑1 (MQ=255)
        gTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAACCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCCACACACAGAAGAAGGACg                                                                                                                                           >  4:32185/1‑149 (MQ=17)
                  cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCa                                                                                                                                  <  2:74111/148‑1 (MQ=255)
                      ggTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTGTCATTCCt                                                                                                                             >  4:150755/1‑149 (MQ=255)
                         cTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTTCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTggg                                                                                                                          >  2:90166/1‑149 (MQ=255)
                            tACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTggg                                                                                                                          <  3:31674/146‑1 (MQ=255)
                            tACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCa                                                                                                                       <  3:115995/149‑1 (MQ=255)
                              cTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGATGAAGGGCGCCTTTTTCATTCCTGGGCCAgg                                                                                                                     >  2:219694/1‑149 (MQ=255)
                              cTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAgg                                                                                                                     <  6:31405/149‑1 (MQ=255)
                                    gAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATc                                                                                                                  >  6:89832/1‑146 (MQ=255)
                                    gAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATc                                                                                                                  <  5:89832/146‑1 (MQ=255)
                                    gAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGa                                                                                                               <  2:58704/149‑1 (MQ=255)
                                       cAGGCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCTAGTCCTCACATACAAGCGGGTTTACCCGCGAAAGGGCCTGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATcc                                                                                                            <  3:32185/149‑1 (MQ=37)
                                        aGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATccc                                                                                                           >  7:179795/1‑149 (MQ=255)
                                             gAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGAGCCGCTGCAGGAGGGCGCCGTTTTGCTTACGGGGGCAGGCTGAGATGCCGTaag                                                                                                      >  4:153542/1‑147 (MQ=18)
                                             gAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACt                                                                                                      >  1:188774/1‑149 (MQ=255)
                                                   aGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCa                                                                                                                 >  6:20420/1‑132 (MQ=255)
                                                   aGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCa                                                                                                                 <  5:20420/132‑1 (MQ=255)
                                                            cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                                       >  2:56645/1‑149 (MQ=255)
                                                            cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                                       >  8:79170/1‑149 (MQ=255)
                                                                   ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGc                                                                                >  4:72513/1‑149 (MQ=255)
                                                                    tccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGcc                                                                               >  3:176266/1‑149 (MQ=255)
                                                                     cccTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCg                                                                              <  7:79170/149‑1 (MQ=255)
                                                                      ccTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGTGCCCGGTAGGCCTCGACAGcctg                                                                             >  4:151791/1‑147 (MQ=255)
                                                                      ccTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGTGCCCGGTAGGCCTCGACAGACgt                                                                             >  8:98197/1‑148 (MQ=255)
                                                                       cTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCgg                                                                             >  8:15636/1‑148 (MQ=255)
                                                                          cAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCgg                                                                             <  5:266743/145‑1 (MQ=255)
                                                                          cAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCgg                                                                             >  6:266743/1‑145 (MQ=255)
                                                                              tCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCt                                                                     >  2:181248/1‑149 (MQ=255)
                                                                                cTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGt                                                                   <  5:185066/149‑1 (MQ=255)
                                                                                  aCCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                                 >  3:148949/1‑149 (MQ=255)
                                                                                      tACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCa                                                                                                                 >  2:135263/1‑97 (MQ=255)
                                                                                      tACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCa                                                                                                                 <  1:135263/97‑1 (MQ=255)
                                                                                      tACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGc                                                             <  1:234582/149‑1 (MQ=255)
                                                                                           gCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgtg                                                         >  2:30405/1‑148 (MQ=255)
                                                                                              ggTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATc                                                     <  3:72513/149‑1 (MQ=255)
                                                                                               gTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                                 >  6:122477/1‑136 (MQ=255)
                                                                                               gTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGttt                                                                 <  5:122477/136‑1 (MQ=255)
                                                                                                   aCCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAg                                                >  4:91877/1‑149 (MQ=255)
                                                                                                    cccGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGc                                               <  8:179795/149‑1 (MQ=255)
                                                                                                    cccGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGc                                               <  1:181248/149‑1 (MQ=255)
                                                                                                       gcgAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAg                                            >  3:19217/1‑149 (MQ=255)
                                                                                                              ggCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTcc                                     <  4:176266/149‑1 (MQ=255)
                                                                                                                ccAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGATTGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAg                                   <  1:90166/149‑1 (MQ=255)
                                                                                                                ccAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAg                                   <  4:19217/149‑1 (MQ=255)
                                                                                                                    aTCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCAcc                               >  5:105021/1‑149 (MQ=255)
                                                                                                                       aGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTgg                            <  4:148949/149‑1 (MQ=255)
                                                                                                                                 aaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCATGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaaa                  <  3:150755/149‑1 (MQ=255)
                                                                                                                                   aaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCt                                                            <  4:61713/105‑1 (MQ=38)
                                                                                                                                   aaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCt                                                            >  3:61713/1‑105 (MQ=38)
                                                                                                                                   aaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTACTGGGCCAGGATCAGATCCCGTACTGCGCCCCGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCTTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaaaa                 >  4:190438/1‑147 (MQ=255)
                                                                                                                                     aaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAAt              <  7:15636/149‑1 (MQ=255)
                                                                                                                                       ggCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCa                                                                           >  2:4494/1‑86 (MQ=38)
                                                                                                                                       ggCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCa                                                                           <  1:4494/86‑1 (MQ=38)
                                                                                                                                           aCCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGac        <  1:30405/149‑1 (MQ=255)
                                                                                                                                             ccTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTc                                      >  7:54777/1‑117 (MQ=255)
                                                                                                                                             ccTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTc                                      <  8:54777/117‑1 (MQ=255)
                                                                                                                                               tctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGCTGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTGCAGCACCTGGGTCAGCGAAACAATGCTGaccacc    >  4:24233/1‑149 (MQ=255)
                                                                                                                                               tctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt           <  6:56828/142‑1 (MQ=255)
                                                                                                                                               tctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt           >  5:56828/1‑142 (MQ=255)
                                                                                                                                                   aagaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCgg  <  2:188774/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGGAGGGGTGATTTCTCCTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGG  >  NC_002947/6039094‑6039387

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: ‑