Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,669,941 G→C 5.9% E79Q (GAG→CAG)  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,9410GC 5.9% 220.5 / 3.8 68E79Q (GAG→CAG) fnrBCrp/Fnr family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (35/29);  new base C (2/2);  total (37/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.25e-01

ACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGACCTG  >  NC_002947/3669801‑3670087
                                                                                                                                            |                                                                                                                                                  
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                                                                                                                   >  5:39601/1‑130 (MQ=255)
                                            ttGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGc                                                                                                                   <  6:39601/130‑1 (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                                                                                                <  1:36634/122‑1 (MQ=255)
                                                                       gCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCt                                                                                                >  2:36634/1‑122 (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                                                                        >  3:186973/1‑125 (MQ=255)
                                                                                            gtgtTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTa                                                                        <  4:186973/125‑1 (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           >  2:266257/1‑145 (MQ=255)
                                                                                                                                     cGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCtg           <  1:266257/145‑1 (MQ=255)
                                                                                                                                      gATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTg                                                             <  4:242764/94‑1 (MQ=255)
                                                                                                                                      gATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGa      >  2:211119/1‑149 (MQ=255)
                                                                                                                                      gATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGCTGCCAGTCATTACCGAAGCTTTg                                                             >  3:242764/1‑94 (MQ=255)
                                                                                                                                          ccGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGAccgg  >  5:364057/1‑147 (MQ=255)
                                                                                                                                          ccGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGACCt   <  1:240055/148‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                  
ACGAAGTGGACAAACTGGAGCGAATTGTCCGCCGTAACCGGCCATTGAAGCGGGGTGAGTACCTGTTCAAGGCCAATGAACCCATGGAGCACGTGTTTGCCTTGCGCTCTGGGGCGATCAAGAATTTTCTGCTCGATGCCGAGGGCACCGAGCGGGTGACCGGCTTCCTCCTGCCGGGTGAAATGCTCGGGCTGGACGCTATTGGTGCCAGTCATTACCGAAGCTTTGCCGTTGCCCTGGAACTCTCGCTGGTCTGCTCCATCAGGCTGGATCAGCTGGTGGACCTG  >  NC_002947/3669801‑3670087

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: