Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 838,006 Δ1 bp 9.5% intergenic (‑39/+41) pth ← / ← rplY peptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947838,0060G. 9.5% 191.6 / 15.8 63intergenic (‑39/+41)pth/rplYpeptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc
Reads supporting (aligned to +/- strand):  ref base G (24/32);  new base . (4/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

ACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCG  >  NC_002947/837857‑838150
                                                                                                                                                      |                                                                                                                                                 
acGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTgcg                                                                                                                                                   >  1:703616/1‑150 (MQ=255)
  gCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTgcgcg                                                                                                                                                 >  1:580075/1‑150 (MQ=255)
           tGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGt                                                                                                                                       <  2:1002913/150‑1 (MQ=255)
           tGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGt                                                                                                                                       >  2:1213700/1‑150 (MQ=255)
             gCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTaa                                                                                                                                     <  1:40365/150‑1 (MQ=255)
             gCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTaa                                                                                                                                     <  1:264569/150‑1 (MQ=255)
               aaGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAAcc                                                                                                                                   <  1:115675/150‑1 (MQ=255)
                aGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAAccc                                                                                                                                  >  1:1033393/1‑150 (MQ=255)
                        aaCGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACg                                                                                                                          >  1:537285/1‑150 (MQ=255)
                           gAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCg                                                                                                                       <  2:525108/150‑1 (MQ=255)
                             aaaGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGc                                                                                                                     >  2:1007357/1‑150 (MQ=255)
                             aaaGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGc                                                                                                                     >  1:55592/1‑150 (MQ=255)
                               agagCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTCTATTGATCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGt                                                                                                                   <  1:1043012/150‑1 (MQ=255)
                                    cccTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAg                                                                                                              <  2:722455/150‑1 (MQ=255)
                                       tGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTgg                                                                                                           >  2:442773/1‑150 (MQ=255)
                                        gCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATGAACGCCGGCGTGCGAGTGGa                                                                                                          <  2:172162/150‑1 (MQ=255)
                                        gCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGa                                                                                                          <  2:834466/150‑1 (MQ=255)
                                         cGTCATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTCGGAGTTGGTGGGTTACAGTGCGCGCGGAACTGTAACCCGATCAACGCCGGCGTGCGAGTGGAt                                                                                                         <  1:340533/150‑1 (MQ=255)
                                               gCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGGTCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGGGTTGGTGGGTTACAGTGCGCGCGGCACTGTGTCGCGATCAACGCCGGCGTGCGAGTGGATTACTCg                                                                                                   <  2:250911/150‑1 (MQ=255)
                                                       tGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGcc                                                                                           <  2:537285/150‑1 (MQ=255)
                                                           cGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttct                                                                                       >  1:451914/1‑150 (MQ=255)
                                                             tATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgc                                                                                     <  2:55592/150‑1 (MQ=255)
                                                                  ggggCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttctg                                                                                >  2:297486/1‑150 (MQ=255)
                                                                      ccGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttctg                                                                                >  2:109627/1‑146 (MQ=255)
                                                                      ccGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttctg                                                                                <  1:109627/146‑1 (MQ=255)
                                                                               aCCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTg                                                                    <  1:819006/150‑1 (MQ=255)
                                                                                cccAGGCCAACGATCAACTGGATGGCGGTCACGTGAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTgg                                                                   >  1:353465/1‑150 (MQ=255)
                                                                                cccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTgg                                                                   >  1:412331/1‑150 (MQ=255)
                                                                                 ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa                                                                              >  1:810207/1‑138 (MQ=255)
                                                                                 ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa                                                                              <  2:810207/138‑1 (MQ=255)
                                                                                 ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGa                                                                  >  1:1030113/1‑150 (MQ=255)
                                                                                 ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGa                                                                  <  1:442773/150‑1 (MQ=255)
                                                                                 ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGa                                                                  >  2:1031661/1‑150 (MQ=255)
                                                                                 ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCCGCGCGCACTGTaa                                                                                                                                     >  1:20735/1‑83 (MQ=38)
                                                                                 ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCCGCGCGCACTGTaa                                                                                                                                     <  2:20735/83‑1 (MQ=38)
                                                                                     gCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGt                                                              <  2:580075/150‑1 (MQ=255)
                                                                                         aCGAGCAACGGGATGGCGGCCACGTAAGGGGCTCTTCCTTGGCGGTGGTGGGTTACAGTGCCCGCGGCACTGTAACGCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAAc                                                          <  2:122302/150‑1 (MQ=255)
                                                                                            aTCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGtt                                                       <  1:1007357/150‑1 (MQ=255)
                                                                                             tCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTg                                                      <  2:1030113/150‑1 (MQ=255)
                                                                                             tCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTg                                                      <  1:110277/150‑1 (MQ=255)
                                                                                                        gCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCttt                                           <  1:1031661/150‑1 (MQ=255)
                                                                                                            tCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa                                                                              >  2:526830/1‑111 (MQ=255)
                                                                                                            tCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa                                                                              <  1:526830/111‑1 (MQ=255)
                                                                                                             cACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGtcat                                      <  2:1121883/150‑1 (MQ=255)
                                                                                                                gTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTgg                                                                                                           <  1:182251/78‑1 (MQ=255)
                                                                                                                gTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTgg                                                                                                           >  2:182251/1‑78 (MQ=255)
                                                                                                                gTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCAc                                   >  2:789278/1‑150 (MQ=255)
                                                                                                                         tctTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCa                          >  2:687137/1‑150 (MQ=255)
                                                                                                                           ttCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCaga                        <  2:412331/150‑1 (MQ=255)
                                                                                                                           ttCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACACGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCaga                        >  1:159127/1‑150 (MQ=255)
                                                                                                                                   aGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGaa                <  2:1033393/150‑1 (MQ=255)
                                                                                                                                     ttGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTACGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAAct              >  2:1089859/1‑150 (MQ=255)
                                                                                                                                         tggGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTAc          <  1:297486/150‑1 (MQ=255)
                                                                                                                                           ggTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGaa                >  2:757401/1‑142 (MQ=255)
                                                                                                                                           ggTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGaa                <  1:757401/142‑1 (MQ=255)
                                                                                                                                           ggTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGc        >  1:777314/1‑150 (MQ=255)
                                                                                                                                            gTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGcc       >  1:558106/1‑150 (MQ=255)
                                                                                                                                            gTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGt                                <  2:746901/125‑1 (MQ=255)
                                                                                                                                            gTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGt                                >  1:746901/1‑125 (MQ=255)
                                                                                                                                            gTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGcc       >  2:397483/1‑150 (MQ=255)
                                                                                                                                             ttACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCt      >  1:344900/1‑150 (MQ=255)
                                                                                                                                                cAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGc                      >  1:264511/1‑131 (MQ=255)
                                                                                                                                                cAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGc                      <  2:264511/131‑1 (MQ=255)
                                                                                                                                                cAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTc   <  1:776893/150‑1 (MQ=255)
                                                                                                                                                 aGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCg  <  1:1122650/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
ACGCTGGGCGCTGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCG  >  NC_002947/837857‑838150

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: