Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 838,010:1 +C 9.5% intergenic (‑43/+37) 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,0101.C 9.5% 241.4 / 19.3 63intergenic (‑43/+37)pth/rplYpeptidyl‑tRNA hydrolase/50S ribosomal protein L25/general stress protein Ctc
Reads supporting (aligned to +/- strand):  ref base . (24/33);  new base C (4/2);  total (28/35)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGAGCTTTC  >  NC_002947/837868‑838158
                                                                                                                                                |                                                                                                                                                    
tGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGt                                                                                                                                               >  2:1213700/1‑150 (MQ=255)
tGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGt                                                                                                                                               <  2:1002913/150‑1 (MQ=255)
  gCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTaa                                                                                                                                             <  1:264569/150‑1 (MQ=255)
  gCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTaa                                                                                                                                             <  1:40365/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                                                                                        <  1:109627/146‑1 (MQ=255)
                                                           ccGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCcttctgcttctg                                                                                        >  2:109627/1‑146 (MQ=255)
                                                                    aCCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTg                                                                            <  1:819006/150‑1 (MQ=255)
                                                                     cccAGGCCAACGATCAACTGGATGGCGGTCACGTGAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTgg                                                                           >  1:353465/1‑150 (MQ=255)
                                                                     cccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTgg                                                                           >  1:412331/1‑150 (MQ=255)
                                                                      ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa                                                                                      <  2:810207/138‑1 (MQ=255)
                                                                      ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGa                                                                                      >  1:810207/1‑138 (MQ=255)
                                                                      ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGa                                                                          >  2:1031661/1‑150 (MQ=255)
                                                                      ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGa                                                                          >  1:1030113/1‑150 (MQ=255)
                                                                      ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGa                                                                          <  1:442773/150‑1 (MQ=255)
                                                                      ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCCGCGCGCACTGTaa                                                                                                                                             <  2:20735/83‑1 (MQ=38)
                                                                      ccAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCCGCGCGCACTGTaa                                                                                                                                             >  1:20735/1‑83 (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                                                                                                                   >  2:182251/1‑78 (MQ=255)
                                                                                                     gTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTgg                                                                                                                   <  1:182251/78‑1 (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                        <  1:757401/142‑1 (MQ=255)
                                                                                                                                ggTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGaa                        >  2:757401/1‑142 (MQ=255)
                                                                                                                                ggTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGc                >  1:777314/1‑150 (MQ=255)
                                                                                                                                 gTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGcc               >  1:558106/1‑150 (MQ=255)
                                                                                                                                 gTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGt                                        >  1:746901/1‑125 (MQ=255)
                                                                                                                                 gTTACAGTGCCGCGCGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGt                                        <  2:746901/125‑1 (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)
                                                                                                                                              cgGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGAGCTTTc  <  2:558106/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGGCAAGGCGTTCAACGAAAAGAGCCCCTGCGTTATGCCGGGTCTGTTCGTATTCGGGGCCGGGGTTACCCAGGCCAACGATCAACTGGATGGCGGTCACGTCAGGGGCTCTTCCTTGGAGTTGGTGGGTTACAGTGCGCGCGGCACTGTAACCCGATCAACGCCGGCGTGCGAGTGGATTACTCGGCAGCGCCTTCTGCTTCTGGAGCAACGCGTGGAGCGTGAACGTTGGCAACAGCTTTGTCATCACCGTGGGCCAGAGCGACGAACTCTACGCCTTTCGGAGCTTTC  >  NC_002947/837868‑838158

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: