Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,669,936 A→G 7.9% D77G (GAT→GGT)  fnrB → Crp/Fnr family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,669,9360AG 7.9% 179.6 / 4.9 64D77G (GAT→GGT) fnrBCrp/Fnr family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base G (2/3);  total (35/29)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

TCAGCGCCGACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGAC  >  NC_002947/3669792‑3670084
                                                                                                                                                |                                                                                                                                                      
tCAGCGCCGACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGGTGcc                                                                                                                                                    >  2:11402/1‑149 (MQ=255)
 cAGCGCCGACGAAGTGGACAAACTGGAGCGAATTGTCAGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCACTGGGGCGATCAAGAATTTTCTGCTCGATGcc                                                                                                                                                    >  2:4460/1‑148 (MQ=255)
  aGCGCCGACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGcc                                                                                                                                                    >  4:85139/1‑147 (MQ=255)
   gcgcCGACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAg                                                                                                                                                 >  2:114666/1‑149 (MQ=255)
    cgcCGACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAgg                                                                                                                                                <  4:369377/149‑1 (MQ=255)
       cgacgaAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGGTGccctcggca                                                                                                                                             >  3:321433/1‑142 (MQ=255)
        gacgaAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCAc                                                                                                                                            >  7:362929/1‑149 (MQ=255)
            aaGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGGTGCCCTCGGCATCGAg                                                                                                                                        <  4:321433/149‑1 (MQ=255)
            aaGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGGTGCCCTCGGCATCGAg                                                                                                                                        <  1:11402/149‑1 (MQ=255)
                  aCAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGt                                                                                                                                   >  1:170399/1‑148 (MQ=255)
                   cAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGGTGCCCTCGGCATCGAGcagaaaa                                                                                                                                 <  4:292699/149‑7 (MQ=255)
                    aaaCTGGAGCGAATTCTCCGCCGTAACCGGCCATTGCAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGAc                                                                                                                                <  5:369373/149‑1 (MQ=255)
                        tGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGc                                                                                                                            >  7:196740/1‑149 (MQ=255)
                          gAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGAc                                                                                                                                <  1:209761/143‑1 (MQ=255)
                          gAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGAc                                                                                                                                >  2:209761/1‑143 (MQ=255)
                          gAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCtt                                                                                                                          >  8:53759/1‑149 (MQ=255)
                               aaTTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTtcctcc                                                                                                                     >  2:327632/1‑149 (MQ=255)
                               aaTTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTtcctc                                                                                                                      <  4:200935/148‑1 (MQ=255)
                               aaTTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTtcctc                                                                                                                      >  3:200935/1‑148 (MQ=255)
                                       gccgTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCAc                                                                                                                                            <  8:310757/118‑1 (MQ=255)
                                       gccgTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCAc                                                                                                                                            >  7:310757/1‑118 (MQ=255)
                                        ccgTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCAGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑Gt                                                                                                            >  4:225485/1‑148 (MQ=255)
                                             aCCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAAt                                                                                                       <  7:53759/148‑1 (MQ=255)
                                               cGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCt                                                                                                    >  2:23047/1‑149 (MQ=255)
                                               cGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCt                                                                                                    >  4:362222/1‑149 (MQ=255)
                                                  ccATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCgg                                                                                                 >  5:222449/1‑149 (MQ=255)
                                                  ccATTGAAGCGGGGTGAGTACCTGTTCAAGGACAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCgg                                                                                                 <  2:357379/149‑1 (MQ=255)
                                                     ttGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCt                                                                                              >  1:353624/1‑149 (MQ=255)
                                                          gctgggTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACg                                                                                         <  8:196740/146‑1 (MQ=255)
                                                                tGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGATGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGCGTGAAATGCTCGGGCTGGACGCTAtt                                                                                    <  1:4460/149‑1 (MQ=255)
                                                                    tACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCGTGCCGG‑GTGATATGCTCGGGCTGGACGCTATTGGTgg                                                                               >  8:254368/1‑148 (MQ=255)
                                                                    tACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGc                                                                               >  2:266783/1‑149 (MQ=255)
                                                                          ttCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCa                                                                         <  1:172116/149‑1 (MQ=255)
                                                                             aaGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGcc                                                                              <  2:91364/141‑1 (MQ=255)
                                                                             aaGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGcc                                                                              >  1:91364/1‑141 (MQ=255)
                                                                             aaGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTa                                                                      <  2:170399/149‑1 (MQ=255)
                                                                                gCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTAcc                                                                    >  4:218414/1‑148 (MQ=255)
                                                                                gCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTAcc                                                                    >  1:224674/1‑148 (MQ=255)
                                                                                  cAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑aa                                                                <  3:362222/149‑1 (MQ=255)
                                                                                   aaTGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAg                                                               >  3:291903/1‑149 (MQ=255)
                                                                                          ccATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGcc                                                        >  5:386539/1‑149 (MQ=255)
                                                                                            aTGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGt                                                      >  4:253203/1‑149 (MQ=255)
                                                                                            aTGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGt                                                      >  7:347479/1‑149 (MQ=255)
                                                                                            aTGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGCGCTGGACGCTATTGGTGGCAGTCATTCCCG‑AAGCTTTGCCGt                                                      >  6:81870/1‑149 (MQ=255)
                                                                                                      tgtTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGaa                                            <  3:253203/149‑1 (MQ=255)
                                                                                                             cTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGCAAGCTTTGCCGTTGCCCTGGAACTCTCg                                      <  1:23047/149‑1 (MQ=255)
                                                                                                             cTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGc                                     <  7:254774/149‑1 (MQ=255)
                                                                                                              ttGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTAc                                                                     >  4:266443/1‑117 (MQ=255)
                                                                                                              ttGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCt                                    >  5:233085/1‑149 (MQ=255)
                                                                                                                gcgcTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTg                                   <  6:233085/148‑1 (MQ=255)
                                                                                                                gcgcTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTg                                   <  1:327632/148‑1 (MQ=255)
                                                                                                                gcgcTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTg                                   <  8:347479/148‑1 (MQ=255)
                                                                                                                 cgcTCTGTGGCGATGAAGAATTTTCTGCTCGTTGCCGAGGGCACCGTGCGGGTGAGCGGCTTGCTCCTGCCGG‑GTGAAATGCTGGGGCTGGACGCTATTGGTGCGAGTGATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGt                                 <  2:94898/149‑1 (MQ=255)
                                                                                                                 cgcTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTgg                                  <  6:386539/148‑1 (MQ=255)
                                                                                                                        gggCGATCAAGAATTTTCTGCTCGATGCCGATGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCt                            <  3:85139/147‑1 (MQ=255)
                                                                                                                        gggCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGc                             <  2:353624/146‑1 (MQ=255)
                                                                                                                          gCGATCAAGAATTTTCTGCTCGATGCCGAGGGAACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑CTGAAATGCTCGGACAGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTAGCCGAAGCCCTGCACATCTCGCTCGTCAGCTACAt                        >  6:354877/1‑149 (MQ=255)
                                                                                                                             aTCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAg                     >  6:330488/1‑149 (MQ=255)
                                                                                                                                aaGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGc                   <  3:218414/148‑1 (MQ=255)
                                                                                                                                aaGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCt                  <  6:222449/149‑1 (MQ=255)
                                                                                                                                aaGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCt                  >  4:170210/1‑149 (MQ=255)
                                                                                                                                aaGAAATTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCt                  >  8:141476/1‑149 (MQ=255)
                                                                                                                                    aTTTTCTGCTCGATGCCGATGTCACCGAGCGGGTGACGTGCTTCCTCCTGCGGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGGCGTGGAACTCTCGCTGGTCTGCTCCATCAGGCTTGAt              <  5:354877/149‑1 (MQ=255)
                                                                                                                                           gCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCtgg       >  4:320327/1‑149 (MQ=255)
                                                                                                                                                aTGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGAc  <  1:266783/149‑1 (MQ=255)
                                                                                                                                                aTGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGAc  <  1:114666/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
TCAGCGCCGACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGG‑GTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCG‑AAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGAC  >  NC_002947/3669792‑3670084

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