Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,039,301 A→G 8.9% intergenic (‑265/+29) 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,3010AG 8.9% 159.7 / 5.2 56intergenic (‑265/+29)PP_5290/pyrEtransposase/orotate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base G (3/2);  total (29/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 = 7.05e-01

CTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGATA  >  NC_002947/6039168‑6039447
                                                                                                                                      |                                                                                                                                                  
cTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAtt                                                                                                                                    <  1:447329/150‑1 (MQ=255)
     gTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                               <  1:417996/150‑1 (MQ=255)
     gTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTAAACGCGAAAGGGCAAGATTAGGGGCGAG‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTgg                                                                                                                               <  1:692568/150‑1 (MQ=18)
         gTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCTTCTTCAGAGGGTCGc                                                                                                                                              <  2:451569/131‑1 (MQ=18)
         gTACTGCGCGAACAGTCAGAAGCGAGGTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCTTCTTCAGAGGGTCGc                                                                                                                                              >  1:451569/1‑131 (MQ=18)
             tGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGAt                                                                                                                       >  2:945892/1‑150 (MQ=255)
               cgcgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCACAAAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATc                                                                                                                      <  2:857623/150‑1 (MQ=255)
                 cgAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAAGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCATGAGCCAGGATCAGa                                                                                                                   >  1:334783/1‑150 (MQ=255)
                        taaGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTa                                                                                                            >  2:1556740/3‑150 (MQ=255)
                                gAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGccc                                                                                                    >  1:439451/1‑150 (MQ=255)
                                  gTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCgg                                                                                                  >  2:333057/1‑150 (MQ=255)
                                     aTGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAg                                                                                               <  1:385033/150‑1 (MQ=255)
                                         cccTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCt                                                                                           <  1:1002488/150‑1 (MQ=255)
                                          cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTc                                                                                          >  1:1543299/1‑150 (MQ=255)
                                          cctcctcctCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTc                                                                                          >  2:979118/1‑150 (MQ=255)
                                            tcctcctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGa                                                                                        >  2:460856/1‑150 (MQ=255)
                                                 ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTg                                                                                                                                <  2:1047875/105‑1 (MQ=255)
                                                 ctccCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTg                                                                                                                                >  1:1047875/1‑105 (MQ=255)
                                                     cTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTGTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCGCGGTAGGCCTCGACAGCCGGCa                                                                               >  1:21826/1‑150 (MQ=255)
                                                          aGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCCGAAACGGCGCTTG‑AAAGGAAAGCGAGCTTCTGAATAAGGTCGCCTTTTTCATTCCTGGGCAAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGc                                                                          >  1:31022/1‑150 (MQ=255)
                                                          aGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGc                                                                          >  1:1220626/1‑150 (MQ=255)
                                                            tCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAt                                                                                                                                     >  2:1007732/1‑89 (MQ=255)
                                                            tCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCAt                                                                                                                                     <  1:1007732/89‑1 (MQ=255)
                                                             ccTACCATACAAGCGGGTTAAGCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGt                                                                       <  2:1613712/150‑1 (MQ=255)
                                                               tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGaagaag                                                                                                                                                   <  1:15128/72‑1 (MQ=255)
                                                               tACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGaagaag                                                                                                                                                   >  2:15128/1‑72 (MQ=255)
                                                                aCCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCttttt                                                                                                                                        <  1:480659/82‑1 (MQ=255)
                                                                aCCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCttttt                                                                                                                                        >  2:480659/1‑82 (MQ=255)
                                                                       aaGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCtgtg                                                             >  1:1540781/1‑150 (MQ=255)
                                                                        aGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttt                                                                                                                                          >  1:274440/1‑72 (MQ=21)
                                                                        aGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCTTCTTCAGAGGGTCGCCttt                                                                                                                                          <  2:274440/72‑1 (MQ=25)
                                                                          cGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGAt                                                          <  1:945892/150‑1 (MQ=255)
                                                                           gggTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATc                                                         <  2:850438/150‑1 (MQ=255)
                                                                                 aCCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGc                                                   <  1:460856/150‑1 (MQ=255)
                                                                                  cccGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGcc                                                  >  2:46447/1‑150 (MQ=255)
                                                                                     gcgAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAgg                                               >  2:45433/1‑150 (MQ=255)
                                                                                     gcgAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAgg                                               >  2:1339557/1‑150 (MQ=255)
                                                                                       gAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGaa                                             <  2:1543299/150‑1 (MQ=255)
                                                                                       gAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGaa                                             <  1:358521/150‑1 (MQ=255)
                                                                                          aGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTAATGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCGGTTTCAGCTGTGGATCGTCAGGCAGGAACTc                                          <  2:495349/150‑1 (MQ=255)
                                                                                             gCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAg                                       >  2:1536347/1‑150 (MQ=255)
                                                                                                    cAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTgg                                <  1:333057/150‑1 (MQ=255)
                                                                                                          gcAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGAATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGc                          <  2:1220626/150‑1 (MQ=255)
                                                                                                            gaatgaaaaaGGCGACCTTCTTCAGAGGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCg                         >  2:1355422/6‑150 (MQ=38)
                                                                                                            ac‑aAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGa                        >  1:643444/1‑150 (MQ=255)
                                                                                                            ac‑aAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTACAGCCCCTGGATcggaga                        >  2:1168859/1‑145 (MQ=255)
                                                                                                               aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaaa                      <  1:1556740/150‑1 (MQ=255)
                                                                                                               aaaaaaaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGaaa                      <  1:979118/150‑1 (MQ=255)
                                                                                                                    aaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATg                 <  2:111641/150‑1 (MQ=255)
                                                                                                                    aaGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATg                 <  2:439451/150‑1 (MQ=255)
                                                                                                                      ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTTGGCCTCGACAGCCGGCAGATGCTGTTTCa                                                                   >  1:715359/1‑98 (MQ=255)
                                                                                                                      ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCa                                                                   <  2:715359/98‑1 (MQ=255)
                                                                                                                      ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt               <  1:1536347/150‑1 (MQ=255)
                                                                                                                      ggCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCt               >  1:782761/1‑150 (MQ=255)
                                                                                                                             ctctGAAGAAGGTCTCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGc                <  2:1326883/142‑1 (MQ=255)
                                                                                                                             ctctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGc                >  1:1326883/1‑142 (MQ=255)
                                                                                                                              tctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCg       >  1:193256/1‑150 (MQ=255)
                                                                                                                              tctGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCg       >  1:1500451/1‑150 (MQ=255)
                                                                                                                                 gaagaaGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGa    <  1:1339557/150‑1 (MQ=255)
                                                                                                                                   agaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGCCCACCGGGATa  <  2:21826/150‑1 (MQ=255)
                                                                                                                                   agaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGATa  >  2:1636969/1‑150 (MQ=255)
                                                                                                                                    gaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCt                                  <  2:295227/117‑1 (MQ=255)
                                                                                                                                    gaagGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCt                                  >  1:295227/1‑117 (MQ=255)
                                                                                                                                      |                                                                                                                                                  
CTACGGTCTGTACTGCGCGAACAGTCAGAAGCGAGTTATGGCCCTCCTCCTCCCTTCAAGTCCTACCATACAAGCGGGTTTACCCGCGAAAGGGCCAGATCAGGCGGCAC‑AAAAAAAGGCGACCCTCTGAAGAAGGTCGCCTTTTTCATTCCTGGGCCAGGATCAGATCCCGTACTGCGCCCGGTAGGCCTCGACAGCCGGCAGATGCTGTTTCAGCTGTGGATCGTCAGCCAGGAACTCCAGCACCTGGGTCAGCGAAACAATGCTGACCACCGGGATA  >  NC_002947/6039168‑6039447

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