Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 574,998 C→A 7.8% A12D (GCC→GAC)  fdoG → formate dehydrogenase‑O major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947574,9980CA 7.8% 202.2 / 6.8 64A12D (GCC→GAC) fdoGformate dehydrogenase‑O major subunit
Reads supporting (aligned to +/- strand):  ref base C (31/28);  new base A (2/3);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.84e-01

ACTAGTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGC  >  NC_002947/574864‑575140
                                                                                                                                      |                                                                                                                                              
aCTAGTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGCCGACGGCGGCGACCtt                                                                                                                                   <  2:235995/148‑1 (MQ=255)
    gTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGa                                                                                                                                <  4:347983/147‑1 (MQ=255)
       ttAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTc                                                                                                                           >  4:199510/1‑149 (MQ=255)
       ttAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTc                                                                                                                           >  6:89545/1‑149 (MQ=255)
       ttAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTc                                                                                                                           >  1:311174/1‑149 (MQ=255)
         aaCCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTc                                                                                                                           <  2:60255/147‑1 (MQ=255)
          aCCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCATCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGa                                                                                                                         <  7:312662/148‑1 (MQ=255)
           ccTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGc                                                                                                                       <  7:275031/149‑1 (MQ=255)
            cTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGcc                                                                                                                      >  4:214938/1‑149 (MQ=255)
                ggcggcCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCTACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGACTCGAGCCTTgc                                                                                                                  >  1:304703/1‑149 (MQ=255)
                    gcCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTgcggcg                                                                                                              >  6:173892/1‑149 (MQ=255)
                          aaTTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTggg                                                                                                         <  2:311174/148‑1 (MQ=255)
                          aaTTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTggg                                                                                                         >  3:81804/1‑148 (MQ=255)
                           aTTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCa                                                                                                       >  8:59722/1‑149 (MQ=255)
                                      ccGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGc                                                                                                                       >  2:145993/1‑122 (MQ=255)
                                      ccGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGc                                                                                                                       <  1:145993/122‑1 (MQ=255)
                                       cGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCtt                                                                                                                                   >  7:316698/1‑109 (MQ=255)
                                       cGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCtt                                                                                                                                   <  8:316698/109‑1 (MQ=255)
                                       cGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGc                                                                                           >  4:192781/1‑149 (MQ=255)
                                               ccAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTc                                                                                   >  7:102025/1‑149 (MQ=255)
                                                      gcgcCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAg                                                                                                                        <  3:149281/105‑1 (MQ=255)
                                                      gcgcCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAg                                                                                                                        >  4:149281/1‑105 (MQ=255)
                                                        gcCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGTGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAg                                                                          <  5:234825/149‑1 (MQ=255)
                                                         tccctattGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAg                                                                          <  5:160964/141‑1 (MQ=255)
                                                          cGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGt                                                                        <  3:199510/149‑1 (MQ=255)
                                                              aaTGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGt                                                                    <  2:146960/149‑1 (MQ=255)
                                                                  ttCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACt                                                                <  8:259479/149‑1 (MQ=255)
                                                                  ttCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACt                                                                >  3:176308/1‑149 (MQ=255)
                                                                     cgcgcgTCCACCACCTGATCAGGAGGGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTc                                                              <  2:160056/148‑1 (MQ=255)
                                                                        gcgTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGc                                                          >  3:110404/1‑149 (MQ=255)
                                                                          gTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTg                                                        <  4:110404/149‑1 (MQ=255)
                                                                             caccacCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTggg                                                                                                         >  8:69753/1‑97 (MQ=255)
                                                                             caccacCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTggg                                                                                                         <  7:69753/97‑1 (MQ=255)
                                                                              accaccTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGgcg                                                     >  7:47101/1‑148 (MQ=255)
                                                                              accaccTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGgcg                                                     >  5:227721/1‑148 (MQ=255)
                                                                               ccaccTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCaa                                                            >  6:229158/1‑140 (MQ=255)
                                                                               ccaccTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCaa                                                            <  5:229158/140‑1 (MQ=255)
                                                                               ccaccTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCa                                                   >  2:347343/1‑149 (MQ=255)
                                                                                  ccTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATAc                                                >  5:336112/1‑149 (MQ=255)
                                                                                   cTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGt                                                                                                             >  5:313345/1‑87 (MQ=255)
                                                                                   cTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGt                                                                                                             <  6:313345/87‑1 (MQ=255)
                                                                                   cTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATAc                                                >  3:13986/1‑148 (MQ=255)
                                                                                      aTCAGGAGTGCACGATGGACCTGATCCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCaaa                                            <  2:55012/149‑1 (MQ=255)
                                                                                      aTCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCaaa                                            <  6:336112/149‑1 (MQ=255)
                                                                                         aGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGaa                                         <  5:173892/149‑1 (MQ=255)
                                                                                           gAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAAc                                       <  6:227721/149‑1 (MQ=255)
                                                                                              tGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCg                                    >  8:312802/1‑149 (MQ=255)
                                                                                               gCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGACGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACt                                                                <  5:209974/120‑1 (MQ=255)
                                                                                               gCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGACGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACt                                                                >  6:209974/1‑120 (MQ=255)
                                                                                                cACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGt                                   >  5:165370/1‑148 (MQ=255)
                                                                                                cACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGt                                   <  1:347343/148‑1 (MQ=255)
                                                                                                    aTGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAacac                              >  8:89785/1‑149 (MQ=255)
                                                                                                        aCCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATa                                                 <  6:97668/126‑1 (MQ=255)
                                                                                                        aCCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATa                                                 >  5:97668/1‑126 (MQ=255)
                                                                                                             aaCCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTa                     <  8:47101/149‑1 (MQ=255)
                                                                                                                cGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTg                  >  6:208144/1‑149 (MQ=255)
                                                                                                                  gcGGCAGTTCATCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGtt                <  2:304703/148‑1 (MQ=255)
                                                                                                                     gCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTtcgg             >  5:212807/1‑149 (MQ=255)
                                                                                                                       aGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTtcggtc           >  7:277766/1‑149 (MQ=255)
                                                                                                                            tccAAGGCCGACGGCGGCGACCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGcccc                                                                                                <  4:354638/57‑1 (MQ=21)
                                                                                                                            tccAAGGCCGACGGCGGCGACCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGcccc                                                                                                >  3:354638/3‑59 (MQ=21)
                                                                                                                             tcAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGc     >  4:212720/1‑149 (MQ=255)
                                                                                                                               aaGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCg    <  7:59722/148‑1 (MQ=255)
                                                                                                                                aGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGc  <  4:155522/149‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                              
ACTAGTCTTAACCTTTGGCGGCCATCAATTGGCCGCCACCGTGACTACCAAGAGGCGCCGGCAATGTTCCGCGCGTCCACCACCTGATCAGGAGTGCACGATGGACCTGAACCGTCGGCAGTTCTTCAAGGTCGCCGCCGTCGGCCTTGGAGGCTCGAGCCTTGCGGCGTTGGGCATGGCCCCAACGCCGGCATTCGCCGAGCAGGTGCGTCACTTCAAGCTGGCGCATACCAAAGAAACCCGTAACACCTGCCCCTACTGTTCGGTCGGCTGCGGC  >  NC_002947/574864‑575140

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: