Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,340,378 C→G 6.1% A76G (GCC→GGC) ‡ folK → 2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,340,3780CG 6.1% 213.0 / 4.0 66A76G (GCC→GGC) ‡folK2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base G (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.21e-01

GACCAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCT  >  NC_002947/5340233‑5340525
                                                                                                                                                 |                                                                                                                                                   
gACCAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGcc                                                                                                                                                    <  1:51367/147‑1 (MQ=255)
gACCAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCa                                                                                                                                                   <  4:138054/148‑1 (MQ=255)
 aCCAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATc                                                                                                                                                 <  8:135952/149‑1 (MQ=255)
  ccAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCCACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATc                                                                                                                                                 >  8:12215/1‑148 (MQ=255)
  ccAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGGCCTGCAGCGCATc                                                                                                                                                 >  5:328588/1‑148 (MQ=255)
   cAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCACCCGCGCTACACCAATGCTGTCGCCGCCCTGGACACAGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCAGCg                                                                                                                                                >  1:64190/1‑148 (MQ=255)
   cAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGGCCTGCAGCGCATCGa                                                                                                                                               <  6:328588/149‑1 (MQ=255)
   cAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGGCCTGCAGCGCATCGa                                                                                                                                               <  3:231537/149‑1 (MQ=255)
      gTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCCCCACTGGCATTGCTCGATGCGCTGCAGGCCATCGaaaa                                                                                                                                            >  7:219550/1‑149 (MQ=255)
              cACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAgg                                                                                                                                    <  6:249997/149‑1 (MQ=255)
                     cTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGt                                                                                                                             >  3:245236/1‑149 (MQ=255)
                     cTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCACCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGt                                                                                                                             >  3:292250/1‑149 (MQ=255)
                         tAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAg                                                                                                                                     >  6:105581/1‑137 (MQ=255)
                         tAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAg                                                                                                                                     <  5:105581/137‑1 (MQ=255)
                            ccgccTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCaaa                                                                                                                      <  7:164814/149‑1 (MQ=255)
                              gccTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAaga                                                                                                                    <  3:153680/149‑1 (MQ=255)
                              gccTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAaga                                                                                                                    <  4:347252/149‑1 (MQ=255)
                                           acacCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCACCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGccc                                                                                                        <  4:292250/148‑1 (MQ=255)
                                             acCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCcacg                                                                                                     >  3:75639/1‑149 (MQ=255)
                                               cAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCcacgca                                                                                                   >  5:98998/1‑149 (MQ=255)
                                               cAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCcacgca                                                                                                   <  5:242273/149‑1 (MQ=255)
                                                aGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCcacgcac                                                                                                  >  6:241411/1‑149 (MQ=255)
                                                aGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCcacgcac                                                                                                  >  4:58151/1‑149 (MQ=255)
                                                  cGATTCCCTGCGGCCCTGCCAACCGCGCTACACCAATGCTGTCGGCGCCCTGGACACCGCCGTGGCACCACTCGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAATGAGCGCTGGGGCCCacgcacgc                                                                                                <  8:81199/149‑1 (MQ=255)
                                                        gctgctGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGctggac                                                                                           <  3:58151/148‑1 (MQ=255)
                                                           gctgcCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGctggacctg                                                                                        >  5:113544/1‑148 (MQ=255)
                                                             tgcCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGcgggacctgga                                                                                      >  7:235030/1‑148 (MQ=255)
                                                                      aaCCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCctgctg                                                                            >  7:109940/1‑149 (MQ=255)
                                                                      aaCCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCctgctg                                                                            <  4:75639/149‑1 (MQ=255)
                                                                       aCCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCCCCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGt                                                                           >  8:62799/1‑149 (MQ=255)
                                                                         cgcgCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTc                                                                         >  2:277207/1‑149 (MQ=255)
                                                                            gcTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCTAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGt                                                                      <  1:277207/149‑1 (MQ=255)
                                                                                  ccAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAg                                                                <  5:241411/149‑1 (MQ=255)
                                                                                        cTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGctgg                                                                                             >  5:247657/1‑114 (MQ=255)
                                                                                        cTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGctgg                                                                                             <  6:247657/114‑1 (MQ=255)
                                                                                           tccccgccCTGGACACCGCCCTGGCACCACTGGCATTGCGCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCGACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGAc                                                       <  8:219550/146‑1 (MQ=255)
                                                                                             gccgccCTGGACACCGCCCTGGCTCCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCa                                                     <  7:62799/149‑1 (MQ=255)
                                                                                             gccgccCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGAcc                                                      <  7:12215/148‑1 (MQ=255)
                                                                                              ccgccCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAg                                                    >  8:44960/1‑149 (MQ=255)
                                                                                               cgccCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGc                                                   <  6:113544/149‑1 (MQ=255)
                                                                                               cgccCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGc                                                   <  3:219080/149‑1 (MQ=255)
                                                                                               cgccCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGc                                                   <  7:44960/149‑1 (MQ=255)
                                                                                                gccCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGcc                                                  >  6:208444/1‑149 (MQ=255)
                                                                                                   cTGGACGCCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCg                                                                        <  6:176858/124‑1 (MQ=255)
                                                                                                   cTGGACGCCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCg                                                                        >  5:176858/1‑124 (MQ=255)
                                                                                                   cTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACg                                               >  7:166910/1‑149 (MQ=255)
                                                                                                            gCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACCTCCTGCTGTTCGGTGACCAGGTGCTCGACCCGCCACGCCTGAAAGt                                      >  5:64816/1‑149 (MQ=255)
                                                                                                             cccTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGCCCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGAcc                                                      >  5:48046/1‑132 (MQ=255)
                                                                                                             cccTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGAcc                                                      <  6:48046/132‑1 (MQ=255)
                                                                                                             cccTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGa                                          >  7:120426/1‑144 (MQ=255)
                                                                                                             cccTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGa                                          <  8:120426/144‑1 (MQ=255)
                                                                                                                tGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGa                                          >  4:132836/1‑141 (MQ=255)
                                                                                                                tGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGa                                          <  3:132836/141‑1 (MQ=255)
                                                                                                                    accacTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCctgct                                                                             <  2:79757/102‑1 (MQ=255)
                                                                                                                    accacTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCctgct                                                                             >  1:79757/1‑102 (MQ=255)
                                                                                                                      cacTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTAc                            >  5:294185/1‑149 (MQ=255)
                                                                                                                           gCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACcac                         <  6:294185/147‑1 (MQ=255)
                                                                                                                           gCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCAcat                       <  1:24671/149‑1 (MQ=255)
                                                                                                                                tttCGATGGCCTGCAGCGCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTGCGGTGACCAGGTGCTCGACCCGCCACGCCGGAAAGTGCCGCATTACCAcgtgcatg                  >  8:348703/3‑149 (MQ=255)
                                                                                                                                      tGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGa                                                        <  7:247663/105‑1 (MQ=255)
                                                                                                                                      tGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGa                                                        >  8:247663/1‑105 (MQ=255)
                                                                                                                                           tGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTc       >  5:356164/1‑149 (MQ=255)
                                                                                                                                            gCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTc                                                          <  5:239962/97‑1 (MQ=255)
                                                                                                                                            gCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTc                                                          >  6:239962/1‑97 (MQ=255)
                                                                                                                                             cAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTc       <  6:98998/147‑1 (MQ=255)
                                                                                                                                              aGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTg    <  6:356164/149‑1 (MQ=255)
                                                                                                                                                gCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCt  >  7:246298/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GACCAGGTCGAAGCCACGCGCCTGGTAGCCGCCTCGGCCCTGTACACCAGCGATTCGCTGCTGCCCGGCCAACCGCGCTACACCAATGCTGTCGCCGCCCTGGACACCGCCCTGGCACCACTGGCATTGCTCGATGCGCTGCAGGCCATCGAAAACGACCAGGGCCGCGTGCGCAAAGAGCGCTGGGGCCCACGCACGCTGGACCTGGACATCCTGCTGTTCGGTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCT  >  NC_002947/5340233‑5340525

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: