Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,340,604 C→G 5.9% C151W (TGC→TGG) ‡ 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,6040CG 5.9% 216.7 / 3.0 68C151W (TGC→TGG) ‡folK2‑amino‑4‑hydroxy‑6‑ hydroxymethyldihydropteridine diphosphokinase
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base G (2/2);  total (31/37)
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.84e-01

GTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGACTATAGGCGTCC  >  NC_002947/5340456‑5340745
                                                                                                                                                    |                                                                                                                                              
gTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGc                                                                                                                                                <  6:240141/149‑1 (MQ=255)
 tGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGCATGcc                                                                                                                                               >  3:42408/1‑149 (MQ=255)
     cAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGtt                                                                                                                                           >  5:243433/1‑149 (MQ=255)
     cAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGtt                                                                                                                                           >  5:191298/1‑149 (MQ=255)
     cAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGtt                                                                                                                                           >  2:331802/1‑149 (MQ=255)
         tGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTtcgtc                                                                                                                                       <  2:165474/149‑1 (MQ=255)
          gCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTtcgtcg                                                                                                                                      >  4:297256/1‑149 (MQ=255)
           cTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATCCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTtcgtcg                                                                                                                                      <  6:24211/148‑1 (MQ=255)
            tCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTtcgtcg                                                                                                                                      >  3:132456/1‑147 (MQ=255)
            tCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTtcgtcg                                                                                                                                      <  4:132456/147‑1 (MQ=255)
              gACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCt                                                                                                                                  >  3:180170/1‑149 (MQ=255)
                ccAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTgg                                                                                                                                <  8:2300/149‑1 (MQ=255)
                  aGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGaa                                                                                                                              <  2:228158/149‑1 (MQ=255)
                   gCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAAc                                                                                                                             >  3:243512/1‑149 (MQ=255)
                    ccACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGAT‑GCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGCACGc                                                                                                                           >  7:285866/1‑149 (MQ=255)
                      aCGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGcc                                                                                                                          <  7:114407/149‑1 (MQ=255)
                      aCGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGcc                                                                                                                          <  4:243512/149‑1 (MQ=255)
                       cGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCt                                                                                                                         <  2:223127/149‑1 (MQ=255)
                       cGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCt                                                                                                                         <  1:324763/149‑1 (MQ=255)
                             aaaGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAg                                                                                                                     <  6:191298/147‑1 (MQ=255)
                                 tGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCggg                                                                                                               <  2:39545/149‑1 (MQ=255)
                                   ccGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTg                                                                                                                        >  7:119633/1‑138 (MQ=255)
                                   ccGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTg                                                                                                                        <  8:119633/138‑1 (MQ=255)
                                      cATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTc                                                                                                          <  4:180170/149‑1 (MQ=255)
                                      cATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTc                                                                                                          >  3:332110/1‑149 (MQ=255)
                                           ccACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGcc                                                                                                     <  2:252798/149‑1 (MQ=255)
                                           ccACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGc                                                                                                      >  6:235922/1‑148 (MQ=255)
                                               atgcatgcCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCa                                                                                                 >  2:37098/1‑149 (MQ=255)
                                                  catgcCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATc                                                                                               <  2:133720/148‑1 (MQ=255)
                                                      cccgcccgTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCAc                                                                                          >  2:98901/1‑149 (MQ=255)
                                                        cgcccgTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCt                                                                                        >  3:225114/1‑149 (MQ=255)
                                                        cgcccgTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAACGGGCATtcgtc                                                                                                                                       >  7:60552/1‑102 (MQ=255)
                                                        cgcccgTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAACGGGCATtcgtc                                                                                                                                       <  8:60552/102‑1 (MQ=255)
                                                          cccgTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCg                                                                                      <  7:117045/149‑1 (MQ=255)
                                                                   gCTGGATCCGCTGGCCGAACTGGGACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCa                                                                             <  8:159488/149‑1 (MQ=255)
                                                                     tGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAAc                                                                                 <  6:337639/143‑1 (MQ=255)
                                                                     tGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGATCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAAc                                                                                 >  5:337639/1‑143 (MQ=255)
                                                                       tATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGcc                                                                                                                                   >  6:5235/1‑91 (MQ=255)
                                                                       tATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGcc                                                                                                                                   <  5:5235/91‑1 (MQ=255)
                                                                        aTCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTaag                                                                        >  5:227721/1‑148 (MQ=255)
                                                                           cGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACAcc                                                                     <  2:180566/149‑1 (MQ=255)
                                                                                gCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑a                                                                <  1:331802/148‑1 (MQ=255)
                                                                                 ccGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATc                                                              <  1:37098/149‑1 (MQ=255)
                                                                                   gAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGt                                                            <  3:297256/149‑1 (MQ=255)
                                                                                    aaCTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTa                                                           <  4:161972/149‑1 (MQ=255)
                                                                                     aCTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTaa                                                          <  1:186935/149‑1 (MQ=255)
                                                                                         gTACCGGGGGATTTCCGCCTGGCCGAT‑GCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAAt                                                      <  8:285866/148‑1 (MQ=255)
                                                                                           aCCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGt                                                                                    <  7:121878/118‑1 (MQ=255)
                                                                                           aCCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGt                                                                                    >  8:121878/1‑118 (MQ=255)
                                                                                                     ttCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAggg                                          >  6:226020/1‑149 (MQ=255)
                                                                                                      tCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGt                                         >  4:333298/1‑149 (MQ=255)
                                                                                                          ccTGTCCGATTGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGAtt                                     <  8:128087/149‑1 (MQ=255)
                                                                                                          ccTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCAGTAACAGGGTGAtt                                     <  4:332110/149‑1 (MQ=255)
                                                                                                              gCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACt                                 <  4:225114/149‑1 (MQ=255)
                                                                                                                  aTGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTccc                             <  4:9547/149‑1 (MQ=255)
                                                                                                                   tGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCca                            >  4:224458/1‑149 (MQ=255)
                                                                                                                           cccTCGCACAACTGCTCGACGAATGCCCGTTCGTCGCCCTGGAACGCCAGTAGACCCGGCGGGACACCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑AACATAACAAAGCGGTAAAAGGCTGATTGACTTCCCAAACCtagg                    >  4:37413/1‑146 (MQ=255)
                                                                                                                             cTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGGTTGACTTCCCACACCTCGCt                   >  5:225864/1‑148 (MQ=255)
                                                                                                                              tCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGCCCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCa                 >  1:218171/1‑149 (MQ=255)
                                                                                                                                gCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACCCCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCtaccg               >  4:160854/1‑145 (MQ=255)
                                                                                                                                gCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCAc                <  5:235922/148‑1 (MQ=255)
                                                                                                                                     aCTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGACta           >  3:14784/1‑148 (MQ=255)
                                                                                                                                      cTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGGTTGACTTCCCACACCTCGCTCACGACtata         <  6:225864/149‑1 (MQ=255)
                                                                                                                                        gCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGACTATAgg       >  8:156117/1‑149 (MQ=255)
                                                                                                                                         cTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTGCTTCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGACTATAg        <  1:98901/148‑1 (MQ=255)
                                                                                                                                         cTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGACTATAg        <  3:224458/147‑1 (MQ=255)
                                                                                                                                           cgacgaACGGGCATTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGa              >  7:49388/1‑139 (MQ=255)
                                                                                                                                           cgacgaACGGGCATTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGa              <  8:49388/139‑1 (MQ=255)
                                                                                                                                            gacgaATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAAAACCCTG‑ATCGTAAAAATGAGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGACTATAGGCGTc   >  2:291454/1‑149 (MQ=255)
                                                                                                                                             acgaATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCCCTCACGACTATAGGCGTcc  >  7:89190/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GTGACCAGGTGCTCGACCAGCCACGCCTGAAAGTGCCGCATTACCACATGCATGCCCGCCCGTTCGTGCTGTATCCGTTGGCCGAACTGGTACCGGGGGATTTCCGCCTGGCCGATGGCCGCGCCCTCGCACAACTGCTCGACGAATGCCCGTTCGTCGGCCTGGAACGCCTGTAGACCGGGCGTTCCAGCCTGCATCGCCACCTCGGTAACGCCAGTAACACCCTG‑ATCGTAACAATGCGGTAACAGGGTGATTGACTTCCCACACCTCGCTCACGACTATAGGCGTCC  >  NC_002947/5340456‑5340745

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