Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,591,131 C→G 6.7% E255Q (GAG→CAG)  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,1310CG 6.7% 195.5 / 3.9 60E255Q (GAG→CAG) PP_1395AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base G (2/2);  total (29/31)
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.42e-01

GCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGCCAAGT  >  NC_002947/1590986‑1591278
                                                                                                                                                  |                                                                                                                                                   
gCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGcctcct                                                                                                                                                   <  4:314100/148‑1 (MQ=255)
 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                                                                                           >  2:61434/1‑149 (MQ=255)
                                                       cACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAg                                                                                           >  6:226939/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                                                                                                                                  <  2:184176/94‑1 (MQ=255)
                                                                       cGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTc                                                                                                                                  >  1:184176/1‑94 (MQ=255)
                                                                       cGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTc                                                                            >  2:260387/1‑148 (MQ=255)
                                                                         ggCTTTCATACCCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAg                                                                         <  8:72753/149‑1 (MQ=255)
                                                                          gCTTTCATACCCCACCCGATACGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTg                                                                                                                                                  <  7:148526/75‑1 (MQ=255)
                                                                          gCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTg                                                                                                                                                  >  8:148526/1‑75 (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                                                                                                           <  8:4502/104‑1 (MQ=255)
                                                                                    cccacccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGa                                                                                                           >  7:4502/1‑104 (MQ=255)
                                                                                      cacccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGa                                                            <  5:226939/149‑1 (MQ=255)
                                                                                       acccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGAc                                                           <  5:202495/149‑1 (MQ=255)
                                                                                       acccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGAc                                                           >  7:2914/1‑149 (MQ=255)
                                                                                       acccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGAc                                                           <  1:135634/149‑1 (MQ=255)
                                                                                         ccGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGAcc                                                          >  6:299875/1‑148 (MQ=255)
                                                                                           gATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTg                                                                                                                                <  2:109150/76‑1 (MQ=255)
                                                                                           gATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTg                                                                                                                                >  1:109150/1‑76 (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                                                                         >  6:329859/1‑111 (MQ=255)
                                                                                                               tC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAg                                                                         <  5:329859/111‑1 (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)
                                                                                                                                                  |                                                                                                                                                   
GCCCGCAACCGAGCCACATCGCGAATCGGCGGCGCGCCGTAAAGCCGGCTGTATTCACGGCTGAACTGCGACGGGCTTTCATAACCCACCCGATAAGCCGCCACCGCCGCCTC‑CAGCCCGTCGTTCAACATCAGCCGCCGCGCCTCCTGCAGGCGCAGCTGCTTCTGGTACTGAAGCGGGCTCATCGAGGTCACGGCCTTGAAGCGGTGGTGCAAGGTCGAGGTGCTGAGGTTGACCTCACGGGCCAGGTCCTCGATACGCAGCGGTTGTTGGAAGTTGTTGTTGAGCCAAGT  >  NC_002947/1590986‑1591278

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