Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 965,227 T→C 6.1% M70T (ATG→ACG)  ptxC → phosphonate transport system permease PtxC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947965,2270TC 6.1% 193.8 / 2.2 65M70T (ATG→ACG) ptxCphosphonate transport system permease PtxC
Reads supporting (aligned to +/- strand):  ref base T (35/26);  new base C (2/2);  total (37/28)
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.49e-01

TCGCGTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTGGTGTGGGCAGCGC  >  NC_002947/965080‑965373
                                                                                                                                                   |                                                                                                                                                  
tCGCGTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATg                                                                                                                                                   >  5:212738/1‑149 (MQ=255)
  gcgTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTc                                                                                                                                                 >  4:364676/1‑149 (MQ=255)
   cgTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGTGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCt                                                                                                                                                <  1:88165/149‑1 (MQ=255)
   cgTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCt                                                                                                                                                <  4:108778/149‑1 (MQ=255)
    gTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTg                                                                                                                                               <  3:21245/149‑1 (MQ=255)
              cTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCAc                                                                                                                                     <  6:212738/149‑1 (MQ=255)
                      gcacgcacTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTgg                                                                                                                             >  2:333371/1‑149 (MQ=255)
                          gcacTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAAt                                                                                                                         <  8:317515/149‑1 (MQ=255)
                                   ggCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTgg                                                                                                                 <  1:333371/148‑1 (MQ=255)
                                    gCAATGGCGGGCTGGGGCAGATGGGAGAGTACGCAGGGCGCTTCCTGCAGCCGGATCTGTCTGCGCGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTCGAAACCCTGGCCACGTCTGCCATGGCCACCTTGTTGTCAATGGTATTggg                                                                                                                >  4:264633/1‑148 (MQ=255)
                                     cAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGa                                                                                                              >  1:104977/1‑149 (MQ=255)
                                       aTGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATg                                                                                                            >  1:35897/1‑149 (MQ=255)
                                        tGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGt                                                                                                           >  3:8541/1‑149 (MQ=255)
                                        tGGAGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGt                                                                                                           >  8:239402/1‑149 (MQ=255)
                                         ggCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCACGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGtt                                                                                                          >  5:32924/1‑149 (MQ=255)
                                          gCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTg                                                                                                         >  8:127472/1‑149 (MQ=255)
                                              gCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTgg                                                                                                     >  4:200156/1‑149 (MQ=255)
                                              gCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTgg                                                                                                     >  6:358530/1‑149 (MQ=255)
                                                   ggCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGc                                                                                                                                             <  4:371210/104‑1 (MQ=255)
                                                   ggCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGc                                                                                                                                             >  3:371210/1‑104 (MQ=255)
                                                      aGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCcg                                                                                             >  3:220632/1‑149 (MQ=255)
                                                           ggCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCcggcggc                                                                                        >  3:301194/1‑149 (MQ=255)
                                                            gCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCcggcggcg                                                                                       <  8:208013/149‑1 (MQ=255)
                                                              gAGTACGCCGGGCGCTTCCTGCAGTCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCggg                                                                                     <  2:28702/149‑1 (MQ=255)
                                                                    gCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCtt                                                                               <  4:8541/149‑1 (MQ=255)
                                                                       gggCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGGTGCTGGCACTGCCGGCGGCGGGGCGCTTCgg                                                                            <  2:2946/149‑1 (MQ=255)
                                                                         gcgcTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGtt                                                                          <  4:220632/149‑1 (MQ=255)
                                                                             ttCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGcc                                                                      >  4:148797/1‑149 (MQ=255)
                                                                               ccTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGcc                                                                      <  2:104977/147‑1 (MQ=255)
                                                                                cTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGtt                                                                   >  2:376880/1‑149 (MQ=255)
                                                                                  gCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTg                                                                  >  6:334635/1‑148 (MQ=255)
                                                                                   cAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTg                                                                  <  3:364676/147‑1 (MQ=255)
                                                                                   cAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGc                                                                 >  3:190236/1‑148 (MQ=255)
                                                                                    aGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCaa                                                               <  3:148797/149‑1 (MQ=255)
                                                                                     gCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCACGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAg                                                              <  6:32924/149‑1 (MQ=255)
                                                                                      ccGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAgg                                                             >  4:156766/1‑149 (MQ=255)
                                                                                      ccGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAgg                                                             >  1:374007/1‑149 (MQ=255)
                                                                                      ccGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTCGCCATGTCTGGCATCGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAgg                                                             >  5:265449/1‑149 (MQ=255)
                                                                                        ggATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAgggg                                                           >  6:88154/1‑149 (MQ=255)
                                                                                             tgtctgCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCcgcg                                                      >  4:4269/1‑149 (MQ=255)
                                                                                                ctgCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGt                                                   >  8:92509/1‑149 (MQ=255)
                                                                                                      ggCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTtgctg                                             <  7:127472/149‑1 (MQ=255)
                                                                                                       gCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTtgctgc                                            <  3:4269/149‑1 (MQ=255)
                                                                                                         aCCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCcggcggcgg                                                                                      >  8:55663/1‑105 (MQ=255)
                                                                                                         aCCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCcggcggcgg                                                                                      <  7:55663/105‑1 (MQ=255)
                                                                                                           cTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCaa                                        >  4:192584/1‑149 (MQ=255)
                                                                                                            tCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAAt                                       >  6:286068/1‑149 (MQ=255)
                                                                                                                   ggTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGc                                >  4:134984/1‑149 (MQ=255)
                                                                                                                   ggTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGc                                >  5:366999/1‑149 (MQ=255)
                                                                                                                       tgGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTCCTCAATGCCTTGCGTg                             <  7:92509/148‑1 (MQ=255)
                                                                                                                         gcgcGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTgg                                                                                                                 <  3:107056/62‑1 (MQ=255)
                                                                                                                         gcgcGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTgg                                                                                                                 >  4:107056/1‑62 (MQ=255)
                                                                                                                         gcgcGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCa                          >  2:308130/1‑149 (MQ=255)
                                                                                                                           gcgGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCAtt                        >  2:223545/1‑149 (MQ=255)
                                                                                                                              gggCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCt                     <  3:200156/149‑1 (MQ=255)
                                                                                                                                 cgcTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGaa                  <  2:35897/149‑1 (MQ=255)
                                                                                                                                 cgcTGGAAACCCTGGCCACGTCTGGCATTGGCACCTTGTTGGCAATGGTTTTTGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGTCCGTTGCAAGGGGCCGCGCGGTTGCTGCTCAATGCCTTGCGTGCCATTCCTGaa                  <  3:264633/149‑1 (MQ=255)
                                                                                                                                      gAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTg               <  3:22466/147‑1 (MQ=255)
                                                                                                                                      gAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTg               >  4:22466/1‑147 (MQ=255)
                                                                                                                                      gAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTGgt             >  5:161394/1‑149 (MQ=255)
                                                                                                                                      gAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTGgt             >  4:24070/1‑149 (MQ=255)
                                                                                                                                        aaCCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTGgtgt           <  5:334635/149‑1 (MQ=255)
                                                                                                                                           ccTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTGGTGTggg        <  5:286068/149‑1 (MQ=255)
                                                                                                                                            cTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTGGTGTggg        >  4:99138/1‑148 (MQ=255)
                                                                                                                                                 cATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTGGTGTGGGCAgcgc  <  3:99138/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TCGCGTATCTGCAGCTGGACCTGCACGCACTGGTCGGCAATGGCGGGCTGGGGCAGATGGGCGAGTACGCCGGGCGCTTCCTGCAGCCGGATCTGTCTGCGGGGCACCTCAAGGCGGTGTGGCGCGGGGCGCTGGAAACCCTGGCCATGTCTGGCATGGGCACCTTGTTGGCAATGGTATTGGGGATGTTGCTGGCACTGCCGGCGGCGGGGCGCTTCGGTTGGCCGTTGCAAGGGGCCGCGCGTTTGCTGCTCAATGCCTTGCGTGCCATTCCTGAACTGGTGTGGGCAGCGC  >  NC_002947/965080‑965373

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