Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,591,135 C→G 6.5% L253L (CTG→CTC)  PP_1395 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,591,1350CG 6.5% 194.1 / 3.1 61L253L (CTG→CTC) PP_1395AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

CCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGCCAAGTGA  >  NC_002947/1590987‑1591280
                                                                                                                                                     |                                                                                                                                                 
cccGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGc                                                                                                                                                   <  6:165991/149‑1 (MQ=255)
  agCAACCGAGGCACATCGCGAATCGGCGGCGCGCAGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAg                                                                                                                                                 <  2:138691/148‑1 (MQ=255)
         gAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAgc                                                                                                                                          <  1:199420/149‑1 (MQ=255)
          aGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAgct                                                                                                                                         <  8:72677/149‑1 (MQ=255)
              acaTCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCtt                                                                                                                                     >  3:128552/1‑149 (MQ=255)
                   gcgAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTTGTTCAACATCAGCCGCTGCGCCTGCAGGAGgcgcg                                                                                                                                            <  1:48776/137‑2 (MQ=255)
                   gcgAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCTGCGCCTGCAGGAGgcgcg                                                                                                                                            >  2:48776/1‑136 (MQ=255)
                   gcgAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTgg                                                                                                                                 <  4:128552/148‑1 (MQ=255)
                    cgAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTa                                                                                                                               >  3:272553/1‑149 (MQ=255)
                           gcggcgCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCCCCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGc                                                                                                                        >  4:27227/1‑149 (MQ=255)
                             ggcgCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCgg                                                                                                                      >  2:135634/1‑149 (MQ=255)
                              gcgcgcCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCggg                                                                                                                     <  1:296203/149‑1 (MQ=255)
                                      tAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGa                                                                                                             <  7:83725/149‑1 (MQ=255)
                                        aaGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAg                                                                                                            <  2:298062/148‑1 (MQ=255)
                                                  tATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTGTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCtt                                                                                                 >  7:72753/1‑149 (MQ=255)
                                                    ttCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTg                                                                                                >  6:236668/1‑148 (MQ=255)
                                                      cACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAg                                                                                             >  6:226939/1‑149 (MQ=255)
                                                      cACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAg                                                                                             >  2:61434/1‑149 (MQ=255)
                                                                 tGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTCACAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGc                                                                                    >  1:73510/1‑148 (MQ=255)
                                                                 tGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGc                                                                                    >  7:188262/1‑147 (MQ=255)
                                                                   cGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAg                                                                                 <  6:206307/148‑1 (MQ=255)
                                                                   cGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAg                                                                                 <  2:73510/148‑1 (MQ=255)
                                                                    gACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGt                                                                               >  6:202495/1‑149 (MQ=255)
                                                                      cGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTc                                                                                                                                    >  1:184176/1‑94 (MQ=255)
                                                                      cGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTc                                                                                                                                    <  2:184176/94‑1 (MQ=255)
                                                                      cGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTc                                                                              >  2:260387/1‑148 (MQ=255)
                                                                        ggCTTTCATACCCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAg                                                                           <  8:72753/149‑1 (MQ=255)
                                                                                tAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCgg                                                                                          >  6:246671/1‑126 (MQ=255)
                                                                                tAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCgg                                                                                          <  5:246671/126‑1 (MQ=255)
                                                                                  acccacccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGt                                                                 <  3:1648/149‑1 (MQ=255)
                                                                                   cccacccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGa                                                                                                             >  7:4502/1‑104 (MQ=255)
                                                                                   cccacccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGa                                                                                                             <  8:4502/104‑1 (MQ=255)
                                                                                     cacccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGa                                                              <  5:226939/149‑1 (MQ=255)
                                                                                      acccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGAc                                                             >  7:2914/1‑149 (MQ=255)
                                                                                      acccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGAc                                                             <  1:135634/149‑1 (MQ=255)
                                                                                      acccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGAc                                                             <  5:202495/149‑1 (MQ=255)
                                                                                        ccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGAcc                                                            >  6:299875/1‑148 (MQ=255)
                                                                                          gATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTg                                                                                                                                  >  1:109150/1‑76 (MQ=255)
                                                                                          gATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTg                                                                                                                                  <  2:109150/76‑1 (MQ=255)
                                                                                                  gccACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAg                                                 <  5:236668/149‑1 (MQ=255)
                                                                                                      ccgccgCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGttc                                             >  5:69848/1‑147 (MQ=255)
                                                                                                         ccgccTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCGCGGGCCAGGTCCTCg                                          <  8:2914/149‑1 (MQ=255)
                                                                                                           gccTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACTGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGAt                                        <  3:27227/149‑1 (MQ=255)
                                                                                                             cTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATAc                                      >  2:95782/1‑149 (MQ=255)
                                                                                                              tC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAg                                                                           <  5:329859/111‑1 (MQ=255)
                                                                                                              tC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAg                                                                           >  6:329859/1‑111 (MQ=255)
                                                                                                                 cAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCa                                   >  8:155140/1‑149 (MQ=255)
                                                                                                                    cccGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACATCACGGGCCAGGTCCTCGATACGCAGCg                                >  6:92071/1‑149 (MQ=255)
                                                                                                                       gtcgtTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGgtt                             >  5:8827/1‑149 (MQ=255)
                                                                                                                          gtTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGgttgtt                          <  4:272553/149‑1 (MQ=255)
                                                                                                                              aaCATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTgg                        >  4:37291/1‑147 (MQ=255)
                                                                                                                               aCATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAg                     >  7:30785/1‑149 (MQ=255)
                                                                                                                                        gctgcgcCTGCAGGAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGttgttg                         >  3:138581/4‑136 (MQ=255)
                                                                                                                                        gctgcgcCTGCAGGAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGttgttg                         <  4:138581/133‑1 (MQ=255)
                                                                                                                                         ccgcGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGa           <  1:95782/149‑1 (MQ=255)
                                                                                                                                         ccgcGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGa           <  1:61434/149‑1 (MQ=255)
                                                                                                                                           gcgcCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGc         <  8:188262/149‑1 (MQ=255)
                                                                                                                                            cgcCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGcc        <  8:30785/149‑1 (MQ=255)
                                                                                                                                                tcctGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGGCCGCGAACGGCAGCGGTTGTGGGATGTTGTTGTTGAGCCAAGt    >  2:19168/1‑149 (MQ=255)
                                                                                                                                                  ctGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGCCAAGTGa  >  7:285367/1‑149 (MQ=255)
                                                                                                                                                    gCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGCCAAGt    <  2:292332/145‑1 (MQ=255)
                                                                                                                                                    gCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGCCAAGt    >  1:292332/1‑145 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGCCAAGTGA  >  NC_002947/1590987‑1591280

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