Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,572,249 T→A 6.7% L124* (TTG→TAG) ‡ pcaB → 3‑carboxy‑cis,cis‑muconate cycloisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,572,2490TA 6.7% 169.2 / 2.7 60L124* (TTG→TAG) ‡pcaB3‑carboxy‑cis,cis‑muconate cycloisomerase
Reads supporting (aligned to +/- strand):  ref base T (28/26);  new base A (2/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGT  >  NC_002947/1572101‑1572391
                                                                                                                                                    |                                                                                                                                              
ggCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGAttt                                                                                                                                                <  6:108744/149‑1 (MQ=255)
  cAACTCGGCCTTTCCGCTTGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGa                                                                                                                                              <  2:204647/149‑1 (MQ=255)
   aaCTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGAt                                                                                                                                             >  7:225601/1‑149 (MQ=255)
         gCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATCAAATCGAggg                                                                                                                                       <  1:121000/149‑2 (MQ=255)
         gCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGAGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGc                                                                                                                                       >  3:121138/1‑149 (MQ=255)
                   tGGTGAAAGCGTTGGGCAAGGTTATCGCCACCGGCGTGCCAGAGGCTGCGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTCGTTATGCAGCTGCGCGATCCCCTCGATTTGATCTAGGCACACCTCgg                                                                                                                              >  7:129918/1‑148 (MQ=255)
                     gTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGc                                                                                                                             >  3:81869/1‑147 (MQ=255)
                            cGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGcc                                                                                                                    <  3:168860/149‑1 (MQ=255)
                             gTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGcc                                                                                                                    >  5:43302/1‑148 (MQ=255)
                                 ggCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCCAGCTTGCCGaggt                                                                                                               <  4:243723/149‑4 (MQ=255)
                                 ggCAAGGTGATCGCCACCGGAGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAc                                                                                                               >  8:129082/1‑149 (MQ=255)
                                   cAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAccc                                                                                                             >  3:311901/1‑149 (MQ=255)
                                   cAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAccc                                                                                                             >  2:63150/1‑149 (MQ=255)
                                     aGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTg                                                                                                           >  1:79976/1‑149 (MQ=255)
                                     aGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTg                                                                                                           <  4:81869/149‑1 (MQ=255)
                                       gTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAcc                                                                                                              >  1:234255/1‑144 (MQ=255)
                                       gTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATAcc                                                                                                              <  2:234255/144‑1 (MQ=255)
                                       gTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGt                                                                                                          <  4:190204/148‑1 (MQ=255)
                                          aTCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCgc                                                                                                       <  2:258598/148‑1 (MQ=255)
                                             gCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCccag                                                                                                     >  7:45786/1‑147 (MQ=255)
                                              ccACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGcagcag                                                                                                  <  2:79976/149‑1 (MQ=255)
                                               cACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAgg                                                                                                 >  1:153214/1‑149 (MQ=255)
                                                     cGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGa                                                                                           <  4:311901/149‑1 (MQ=255)
                                                      gTGCCAGAGGCTGAGCGCTATGTGCACATTGGGGCCACCCTCCAGGATGCGATGGATACCGGTATGGTTATGCAGCTGCGAGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGACGATACCCTGTCGCAGCCGGCCTTg                                                                                            >  5:173423/1‑147 (MQ=255)
                                                         ccAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTCGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCa                                                                                       >  4:304050/1‑149 (MQ=255)
                                                            gagGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGc                                                                                    <  4:317804/149‑1 (MQ=255)
                                                                   aGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTCGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATa                                                                               <  3:304050/147‑1 (MQ=255)
                                                                    gcgcTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGc                                                                            <  2:153214/149‑1 (MQ=255)
                                                                    gcgcTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGc                                                                            <  2:144380/149‑1 (MQ=255)
                                                                          tgtgCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTtgg                                                                      <  6:43302/149‑1 (MQ=255)
                                                                           gtgCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCTGCAGGCCTTGAAGCACGCCGATACGCCCTtg                                                                       >  7:228672/1‑147 (MQ=255)
                                                                           gtgCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCtt                                                                        >  3:220143/1‑146 (MQ=255)
                                                                           gtgCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCtt                                                                        <  4:220143/146‑1 (MQ=255)
                                                                            tgCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCCTCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCTCAGCAGGCCTTGAAGCACGCCGATACGCCCTtggt                                                                     <  7:129082/148‑1 (MQ=255)
                                                                             gCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTtggtgg                                                                   >  6:184592/1‑149 (MQ=255)
                                                                              cACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCGTTGAAGCACGCCGATACGCCCTtggtgg                                                                   <  1:183209/148‑1 (MQ=255)
                                                                               aCCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGt                                                                 >  2:277128/1‑149 (MQ=255)
                                                                                ccTTGGGGTCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTc                                                                <  5:284971/149‑1 (MQ=255)
                                                                                 cTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCg                                                               >  7:85946/1‑149 (MQ=255)
                                                                                    ggggCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTAc                                                            >  4:155958/1‑149 (MQ=255)
                                                                                      ggCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCt                                                          >  6:284993/1‑149 (MQ=255)
                                                                                           ccagccagGATGCTATGGTTATCGGTCTGGTTATGCAGCTTGGCGATGCCCTCGATTTGTTCGAGGCCGACGTCGGCAAGGTGGCTGATACCCTGTCGCAGCAGGGCTTGAATCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTg                                                     <  8:228672/149‑1 (MQ=255)
                                                                                            cagccagGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGc                                                    <  4:121138/149‑1 (MQ=255)
                                                                                                cagGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGTTTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAaca                                                <  3:270474/149‑1 (MQ=255)
                                                                                                 agGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAgg                                                                                                                                        <  6:142952/60‑1 (MQ=255)
                                                                                                 agGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAgg                                                                                                                                        >  5:142952/1‑60 (MQ=255)
                                                                                                        gATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCAccc                                        >  5:128779/1‑149 (MQ=255)
                                                                                                            gATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGAGCGATGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTTGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGt                                    <  6:261587/149‑1 (MQ=255)
                                                                                                                ccGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGAc                                 <  1:63150/148‑1 (MQ=255)
                                                                                                                   gTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACCCGCCACCCCGGTGCCCCTg                             >  5:248600/1‑149 (MQ=255)
                                                                                                                                 cTGCGCGATGCCCTCGATTTGATCGATGCCGAGCTCGGCAAGCGGTCAGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTTGGGCGTACCTGGTGGCAACACGCCACCCCTGTGTCCCTGGGCATGAAACTGGc               <  6:173423/149‑1 (MQ=255)
                                                                                                                                 cTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTgg                            >  2:267646/1‑136 (MQ=255)
                                                                                                                                 cTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTgg                            <  1:267646/136‑1 (MQ=255)
                                                                                                                                       gtcggccTCGATCAAATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTgg                            >  6:270394/6‑130 (MQ=255)
                                                                                                                                       gtcggccTCGATCAAATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTgg                            <  5:270394/125‑1 (MQ=255)
                                                                                                                                         tGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTAc       <  6:128779/149‑1 (MQ=255)
                                                                                                                                          ggccTCGATCAAATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACt      >  5:288861/3‑149 (MQ=255)
                                                                                                                                            ccTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGAGGGTCGTACCAGGTTGCAAAACGCCACCCCCGTCACCCAGGGCATCAAACTGGCCGGTGTACTgg    >  4:8759/1‑149 (MQ=255)
                                                                                                                                             cTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGGCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCACCACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTggg   >  6:266516/1‑149 (MQ=255)
                                                                                                                                             cTCGATTTGATCGAGGCCGAACTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTCGTTGCAACAAGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTggg   >  3:223870/1‑149 (MQ=255)
                                                                                                                                              tCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGt  >  5:159457/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GGCAACTCGGCCATTCCGCTGGTGAAAGCGTTGGGCAAGGTGATCGCCACCGGCGTGCCAGAGGCTGAGCGCTATGTGCACCTTGGGGCCACCAGCCAGGATGCGATGGATACCGGTCTGGTTCTGCAGCTGCGCGATGCCCTCGATTTGATCGAGGCCGACCTCGGCAAGCTGGCCGATACCCTGTCGCAGCAGGCCTTGAAGCACGCCGATACGCCCTTGGTGGGTCGTACCTGGTTGCAACACGCCACCCCGGTGACCCTGGGCATGAAACTGGCCGGTGTACTGGGT  >  NC_002947/1572101‑1572391

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