Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,813,854 G→C 5.7% R70P (CGT→CCT) ‡ frmC → S‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,813,8540GC 5.7% 227.6 / 3.8 70R70P (CGT→CCT) ‡frmCS‑formylglutathione hydrolase/S‑lactoylglutathione hydrolase
Reads supporting (aligned to +/- strand):  ref base G (36/30);  new base C (2/2);  total (38/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

ATCATTCCAAGGTGCTGGGTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCACCCAGCAACCCTGGGCCCAGCAC  >  NC_002947/1813707‑1814002
                                                                                                                                                   |                                                                                                                                                    
atcatTCCAAGGTGCTGGGTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGt                                                                                                                                                     <  8:166797/149‑1 (MQ=255)
    ttCCAAGGTGCTGGGTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTgg                                                                                                                                                 <  7:196000/149‑1 (MQ=255)
       cAAGGTGCTGGGTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGgccg                                                                                                                                              >  7:45708/1‑149 (MQ=255)
       cAAGGTGCTGGGTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGgcc                                                                                                                                               <  5:16202/148‑1 (MQ=255)
          ggTGCTGGGTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGgccgcc                                                                                                                                            >  1:284266/1‑148 (MQ=255)
           gTGCTGGGTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGGCCAGACGCTGGgcgccgg                                                                                                                                          <  4:258993/149‑2 (MQ=255)
                gggTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTgg                                                                                                                                     >  2:306380/1‑149 (MQ=255)
                 ggTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTggg                                                                                                                                    <  3:88672/149‑1 (MQ=255)
                  gTTGCGATATGGTGATCGCTGTGAACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTgggg                                                                                                                                   <  3:240161/149‑1 (MQ=255)
                      cGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGa                                                                                                                               <  8:210209/149‑1 (MQ=255)
                       gATATGGTGTTCGCTGTGTACCAGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGa                                                                                                                               >  5:148157/1‑148 (MQ=255)
                        atatGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGatc                                                                                                                             <  5:25134/149‑1 (MQ=255)
                              tgtTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCAtcgtc                                                                                                                       <  8:268316/149‑1 (MQ=255)
                               gtTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCAtcgtcg                                                                                                                      >  3:111762/1‑149 (MQ=255)
                                        gtACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATa                                                                                                             >  5:66100/1‑149 (MQ=255)
                                             tGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGc                                                                                                        >  8:116631/1‑149 (MQ=255)
                                             tGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGc                                                                                                        >  4:172052/1‑149 (MQ=255)
                                                      aGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGc                                                                                               <  4:111762/149‑1 (MQ=255)
                                                       ggCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATa                                                                                                             >  7:98350/1‑134 (MQ=255)
                                                       ggCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATa                                                                                                             <  8:98350/134‑1 (MQ=255)
                                                        gCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGa                                                                                             >  8:191061/1‑149 (MQ=255)
                                                             gCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAgg                                                                                        >  1:75832/1‑149 (MQ=255)
                                                                gggCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGc                                                                                     <  1:306380/149‑1 (MQ=255)
                                                                   cGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCgg                                                                                  >  4:87744/1‑149 (MQ=255)
                                                                    gAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCggg                                                                                 <  6:148157/149‑1 (MQ=255)
                                                                    gAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCggg                                                                                 >  5:107506/1‑149 (MQ=255)
                                                                      aaaGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCa                                                                                          <  7:152707/138‑1 (MQ=255)
                                                                      aaaGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCa                                                                                          >  8:152707/1‑138 (MQ=255)
                                                                        aGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCg                                                                                   <  6:4129/143‑1 (MQ=255)
                                                                        aGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCg                                                                                   >  5:4129/1‑143 (MQ=255)
                                                                         gCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATc                                                                            <  2:284266/149‑1 (MQ=255)
                                                                            gCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCg                                                                          >  7:276463/1‑148 (MQ=255)
                                                                                     gTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCACCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCt                                                                >  5:62955/1‑149 (MQ=255)
                                                                                      tACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTg                                                               <  2:75832/149‑1 (MQ=255)
                                                                                         tGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGa                                                            >  5:245253/1‑149 (MQ=255)
                                                                                           gCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGa                                                            <  7:29634/147‑1 (MQ=255)
                                                                                                gCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGc                                                     >  4:70960/1‑149 (MQ=255)
                                                                                                gCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGc                                                     <  6:66100/149‑1 (MQ=255)
                                                                                                cgtgCCTCACCTGCACCGACGACAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTTATCATCGTCGCCGCCGATACCAGCCCGCGTGGCGAGCTGGGGCCGTGCGATCCGGACGGCGCCTGGGACTTCg                                                       <  7:85229/144‑1 (MQ=255)
                                                                                                ccGGGCTCCCCTGCACCGACGAGCAGTTCAGGCAGAAGGAAGGCGCCCCGCGTCTGTCCGGCGCGCTGGGGCTGATCATCGTCGCGCCCGATACCAGACGGCCGGGGGAGCACGTGCCGCGCGATCAGGACGGCGCCTGGGAGTTct                                                       >  8:85229/2‑146 (MQ=255)
                                                                                                  ggCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGa                                                                        >  8:167571/1‑128 (MQ=255)
                                                                                                  ggCCTCACCTGCACCGACGAGAACTTCATCCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGa                                                                        <  7:167571/128‑1 (MQ=255)
                                                                                                     cTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTggg                                                >  5:97882/1‑149 (MQ=255)
                                                                                                      tCACCTCCACAGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCCATACCAGCCCGCGTGGCGAGCAGGTCCCGGGCGATCCGGACGGCGCATGGGACTTCGGCCTGGgc                                               >  8:36516/1‑149 (MQ=255)
                                                                                                        aCCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACg                                                                      >  2:248179/1‑124 (MQ=255)
                                                                                                        aCCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACg                                                                      <  1:248179/124‑1 (MQ=255)
                                                                                                         ccTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCt                                            <  8:276463/149‑1 (MQ=255)
                                                                                                             cACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGc                                         <  3:70960/148‑1 (MQ=255)
                                                                                                             cACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCt                                        <  8:184024/149‑1 (MQ=255)
                                                                                                              aCCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCt                                        >  4:54979/1‑148 (MQ=255)
                                                                                                                cgacgaGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGc                                                                                           >  3:272812/1‑95 (MQ=255)
                                                                                                                 gacgaGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGTTGGCTTCTa                                    <  3:87744/149‑1 (MQ=255)
                                                                                                                  acgaGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTAc                                   >  4:173484/1‑149 (MQ=255)
                                                                                                                   cgaGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTAcc                                  <  8:45708/149‑1 (MQ=255)
                                                                                                                   cgaGAACATCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGCGGCTGATCATCGTCCCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTGCTccc                                  >  3:270704/1‑149 (MQ=255)
                                                                                                                       aaCTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTgg                                                 >  4:201871/1‑130 (MQ=255)
                                                                                                                        aCTTCATGCAGAACGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAAc                             <  7:191061/149‑1 (MQ=255)
                                                                                                                         cTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACg                            >  7:58948/1‑149 (MQ=255)
                                                                                                                          ttCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGc                           >  3:275281/1‑149 (MQ=255)
                                                                                                                             aTGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCAc                        <  7:116631/149‑1 (MQ=255)
                                                                                                                              tGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCAcc                       >  7:184136/1‑149 (MQ=255)
                                                                                                                                cAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCAcc                       >  5:302140/1‑147 (MQ=255)
                                                                                                                                       tcggcggCCAGACGCTGGGCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGc                                         <  5:179613/121‑1 (MQ=255)
                                                                                                                                       tcggcggCCAGACGCTGGGCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGc                                         >  6:179613/2‑122 (MQ=255)
                                                                                                                                       tcggcggCCAGACGCTGGGCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCACCCAGCAAccc              >  4:20759/2‑149 (MQ=255)
                                                                                                                                          gcgcCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCACCCAGCAACCCTgg           <  3:172052/149‑1 (MQ=255)
                                                                                                                                             cccAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGTGCCTGGGACTTCGGCCTGGGCGCTGGc                                         <  8:272287/116‑1 (MQ=255)
                                                                                                                                             cccAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGTGCCTGGGACTTCGGCCTGGGCGCTGGc                                         >  7:272287/1‑116 (MQ=255)
                                                                                                                                               cAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCAACCAGCAACCCTGGGCCCa      >  6:266166/1‑149 (MQ=255)
                                                                                                                                               cAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCAGCCTGGGCGCTGGCTTCTACCTCAACGCCACCCAGCAACCCTGGGCCCa      >  6:84644/1‑149 (MQ=255)
                                                                                                                                                   gTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCACCCAGCAACCCTGGGCCCAGCAc  <  4:275281/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
ATCATTCCAAGGTGCTGGGTTGCGATATGGTGTTCGCTGTGTACCTGCCGCCGCAGGCCGAGCAGGGCGAAAAGCTGCCGGTGCTGTACTGGCTCAGCGGCCTCACCTGCACCGACGAGAACTTCATGCAGAAGGCCGGCGCCCAGCGTCTGGCCGCCGAGCTGGGGCTGATCATCGTCGCCCCCGATACCAGCCCGCGTGGCGAGCAGGTGCCGGGCGATCCGGACGGCGCCTGGGACTTCGGCCTGGGCGCTGGCTTCTACCTCAACGCCACCCAGCAACCCTGGGCCCAGCAC  >  NC_002947/1813707‑1814002

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