Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,669,936 A→G 6.0% 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 6.0% 193.7 / 2.4 66D77G (GAT→GGT) fnrBCrp/Fnr family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (34/28);  new base G (2/2);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.57e-01

ACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGA  >  NC_002947/3669801‑3670083
                                                                                                                                       |                                                                                                                                                   
acgaAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCAcc                                                                                                                                        <  1:153988/149‑1 (MQ=255)
acgaAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCAcc                                                                                                                                        >  6:141516/1‑149 (MQ=255)
  gaAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGa                                                                                                                                      <  3:179795/149‑1 (MQ=255)
   aaGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAg                                                                                                                                     >  5:373566/1‑149 (MQ=255)
    aGTGGACAAACTGGTGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGGTGCCCTCGGCATCGAGc                                                                                                                                    <  8:43416/149‑1 (MQ=255)
        gACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGt                                                                                                                                >  8:239991/1‑149 (MQ=255)
              cTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCgg                                                                                                                          >  3:154561/1‑149 (MQ=255)
               tGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGc                                                                                                                         >  3:363474/1‑149 (MQ=255)
                 gAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGc                                                                                                                         <  4:154561/147‑1 (MQ=255)
                   gCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGGTGCCCTCGGCATCGAGc                                                                                                                                    >  1:172092/1‑134 (MQ=255)
                   gCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGGTGCCCTCGGCATCGAGc                                                                                                                                    <  2:172092/134‑1 (MQ=255)
                       aTTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTcctcct                                                                                                                 >  3:233509/1‑149 (MQ=255)
                             cgccgTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCgg                                                                                                            <  5:141516/148‑1 (MQ=255)
                                cgTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGAACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCTAGCGGGTGACCGGCTTCCTCCTGCGCGGTg                                                                                                         >  3:221404/1‑148 (MQ=255)
                                 gTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTtc                                                                                                                      <  7:305383/134‑1 (MQ=255)
                                      cGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCt                                                                                                  >  8:321280/1‑149 (MQ=255)
                                        gCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCg                                                                                                <  6:373566/149‑1 (MQ=255)
                                        gCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCg                                                                                                <  4:233509/149‑1 (MQ=255)
                                         ccATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATTGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCgg                                                                                               <  4:363474/149‑1 (MQ=255)
                                            ttGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGc                                                                                                               <  6:39601/130‑1 (MQ=255)
                                            ttGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGc                                                                                                               >  5:39601/1‑130 (MQ=255)
                                                       tGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTg                                                                                 >  2:129833/1‑149 (MQ=255)
                                                        gAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTgg                                                                                <  2:113092/149‑1 (MQ=255)
                                                           tACCTGTTCAAGGCCAATGAACCCATGTAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGc                                                                             >  2:240055/1‑149 (MQ=255)
                                                                gTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGt                                                                         >  8:227844/1‑148 (MQ=255)
                                                                  tCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCgg                                                                                                                                  <  8:147667/89‑1 (MQ=255)
                                                                  tCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCgg                                                                                                                                  >  7:147667/1‑89 (MQ=255)
                                                                     aGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTAc                                                                   >  4:192392/1‑149 (MQ=255)
                                                                       gCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCt                                                                                            >  2:36634/1‑122 (MQ=255)
                                                                       gCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCt                                                                                            <  1:36634/122‑1 (MQ=255)
                                                                       gCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTAcc                                                                  >  4:324663/1‑148 (MQ=255)
                                                                             gAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCtt                                                           >  2:321465/1‑149 (MQ=255)
                                                                              aaCCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCtt                                                           <  3:192392/148‑1 (MQ=255)
                                                                               aCCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTg                                                         <  7:239991/149‑1 (MQ=255)
                                                                                     ggAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGaaa                                                                                                      <  3:380083/98‑1 (MQ=255)
                                                                                     ggAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGaaa                                                                                                      >  4:380083/1‑98 (MQ=255)
                                                                                          aCGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTg                                              >  7:5485/1‑149 (MQ=255)
                                                                                            gtgtTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTa                                                                    <  4:186973/125‑1 (MQ=255)
                                                                                            gtgtTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTa                                                                    >  3:186973/1‑125 (MQ=255)
                                                                                            gtgtTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTgg                                             >  1:29036/1‑148 (MQ=255)
                                                                                            gtgtTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTgg                                             >  8:249315/1‑148 (MQ=255)
                                                                                             tgtTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGaa                                           >  6:405992/1‑149 (MQ=255)
                                                                                               tttGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTc                                                                        >  5:369471/1‑118 (MQ=255)
                                                                                               tttGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTc                                                                        <  6:369471/118‑1 (MQ=255)
                                                                                                   ccTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACtctc                                      <  3:312649/148‑1 (MQ=255)
                                                                                                       gcgcTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTg                                  <  7:227844/148‑1 (MQ=255)
                                                                                                       gcgcTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTg                                  <  1:321465/148‑1 (MQ=255)
                                                                                                        cgcTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGt                                >  6:74716/1‑149 (MQ=255)
                                                                                                             tGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCGTTACCGAAGCTGTGCCGTTGCCCTGGAACTCTCGCTGGTCTGc                            >  8:297832/1‑148 (MQ=255)
                                                                                                             tGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCt                           >  1:23437/1‑149 (MQ=255)
                                                                                                               gggCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTcc                         <  5:405992/149‑1 (MQ=255)
                                                                                                               gggCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTc                          >  8:26812/1‑148 (MQ=255)
                                                                                                                     tCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCAGGAACTCTCTCTGGTCTGCTCCATCAgg                   >  6:142207/1‑149 (MQ=255)
                                                                                                                       aaGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCt                 >  4:221551/1‑149 (MQ=255)
                                                                                                                        aGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAAc                                          >  1:241588/1‑123 (MQ=255)
                                                                                                                        aGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAAc                                          <  2:241588/123‑1 (MQ=255)
                                                                                                                           aTTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGAt             <  3:324663/149‑1 (MQ=255)
                                                                                                                                 cgCTCGGTGCCCGCGGCATCGAGCGGGTGGCCGGCTTCCTCCTGCCGGGTGAAATGCTGGGGCTGGCCGCGATTGGTGCCAGGCCTTACCGAAGCTTTTCCGTTGCCCGGGAAAGAGAGCTGGTCTGCTCCCTAACGCTGGATCAGcag       >  4:276289/2‑147 (MQ=255)
                                                                                                                                  gcgcGATGCCGCGGGCCCCGAGCGGGTGACCGGCTTCCTCCTGCCGGGGGAAATGCGCGGGCTGGCCGCTATTGGTGCCAGGCATTACCGACGCTTTGCCGTTGCCCTGGAACTCTCGCGGGTCTGCTCCATCAGGCTGGATCAGCtgg      <  4:87658/146‑1 (MQ=255)
                                                                                                                                   cTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCtg       <  3:221551/147‑1 (MQ=255)
                                                                                                                                     cGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCtggtgg   <  7:321280/149‑1 (MQ=255)
                                                                                                                                     cGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCtg       <  1:266257/145‑1 (MQ=255)
                                                                                                                                     cGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCtg       >  2:266257/1‑145 (MQ=255)
                                                                                                                                      gATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTg                                                         <  4:242764/94‑1 (MQ=255)
                                                                                                                                      gATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGa  >  2:211119/1‑149 (MQ=255)
                                                                                                                                      gATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGCTGCCAGTCATTACCGAAGCTTTg                                                         >  3:242764/1‑94 (MQ=255)
                                                                                                                                       |                                                                                                                                                   
ACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGA  >  NC_002947/3669801‑3670083

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