Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 557,288 G→A 7.6% R152C (CGT→TGT)  folD ← bifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913557,2880GA 7.6% 162.9 / 3.7 65R152C (CGT→TGT) folDbifunctional 5,10‑methylene‑tetrahydrofolate dehydrogenase/ 5,10‑methylene‑tetrahydrofolate cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base G (30/30);  new base A (3/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

TCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACC  >  NC_000913/557140‑557432
                                                                                                                                                    |                                                                                                                                                
tCTACGTGATGACGCAGATTTTTCGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCGCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACg                                                                                                                                                  <  2:932024/149‑1 (MQ=255)
  tACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCt                                                                                                                                                >  1:70671/1‑149 (MQ=255)
   aCGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCCCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTc                                                                                                                                               >  1:365847/1‑149 (MQ=255)
     gTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGACAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCGCaa                                                                                                                                             >  2:504242/1‑149 (MQ=255)
      agaggaCGCAGATTTTTCGGGAACGGGTGAGGCACAGTCGAGGGGCAACCGGCCAGCATAAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAag                                                                                                                                            <  1:441834/145‑1 (MQ=255)
       gatgaCGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAagc                                                                                                                                           >  2:8189/1‑149 (MQ=255)
             gCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCCTCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGCAGGTATCAATGTTGTAACGCTCAAGCAGCGTg                                                                                                                                     >  2:712505/1‑149 (MQ=255)
               aGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGAc                                                                                                                                   >  2:106903/1‑149 (MQ=255)
                 atatttAGGGCAGCGGTGAGGCACTGTAGTGGTGCACCCTGCCAGCAGCAGTTCCAAGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGa                                                                                                                                 <  1:869637/146‑1 (MQ=255)
                        gTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACACTCAAGCAGCGTGACGATACcgcg                                                                                                                          >  2:304858/1‑149 (MQ=255)
                        gTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACACTCAAGCAGCGTGACGATACcgcg                                                                                                                          >  2:993938/1‑149 (MQ=255)
                        gTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAACGTATCTATGTTGTAACGCTCAAGCAGCGTGACGATGCcgcg                                                                                                                          >  2:809359/1‑149 (MQ=255)
                                 tGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCCGCGTGACGATACCGCGCGGGGTGCa                                                                                                                 >  1:451676/1‑149 (MQ=255)
                                       aCTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACg                                                                                                           >  2:856796/1‑149 (MQ=255)
                                        cTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGc                                                                                                          <  1:545555/149‑1 (MQ=255)
                                            aGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGAc                                                                                                      <  1:106903/149‑1 (MQ=255)
                                              tggtgCACCCTGCCTGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCTCCAATCACCACGGCGTTGAGGCCGTAGGTATCAATGTTGTAACGCTCCAGCGGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGCcgc                                                                                                    >  2:565535/1‑149 (MQ=255)
                                              tggtgCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGAcgc                                                                                                    >  2:88004/1‑149 (MQ=255)
                                                 tgCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACgcggc                                                                                                 >  1:1017163/1‑149 (MQ=255)
                                                  gCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCATCGATCTTCGATGCGCCAATGACCACGGCGTGGAGGCGGAAGGTATCAATGTTGTAACGCTCCAGCAGCGTGCCGATACCGCGCGGGGTGCGGGGGCGCAGCCgcggcg                                                                                                >  2:758156/1‑149 (MQ=255)
                                                   cAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGCCGCGgcgc                                                                                               >  2:326050/1‑149 (MQ=255)
                                                         gCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTg                                                                                         <  2:296680/149‑1 (MQ=255)
                                                          cctgcagcagTTCCATGCGCTTCGGGCGGCGAACGATATTCGATGCGCCAATCACCACGGCGTTGCGGCCGAAGGTATCATTGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCTGGGACGCAGACGCGGCGCGCGCTgg                                                                                        <  1:121651/146‑1 (MQ=255)
                                                            agcagcagGTCCATGCTCATCGGGCGGCCACCGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGGGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGca                                                                                      <  2:274644/149‑1 (MQ=255)
                                                             gcagaagTTCCATGCTCATCGGGCGGCCAACGATATTCGCTGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGcac                                                                                     <  2:770537/149‑1 (MQ=255)
                                                                 ccgTTCCATGCGCATCCGGCGGCCAACGAAAGTCGATGCGCCAATCACCACGGCGTTGAGGCCGATGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGGGCAGGGACGCAGACGCGGCGCGCGCTGGCACAgac                                                                                 <  2:751634/147‑1 (MQ=255)
                                                                   gTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAgacga                                                                               <  1:86694/149‑1 (MQ=255)
                                                                   gTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACACTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAgacga                                                                               <  1:304858/149‑1 (MQ=255)
                                                                   gTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACACTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAgacga                                                                               <  1:993938/149‑1 (MQ=255)
                                                                    ttCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCGCAAGCAGCGTGACGATACCGCGCTGGGTGAAGGGACGCAGACGCTGCGCGCGCTGGCACCgacgac                                                                              >  2:670316/1‑149 (MQ=255)
                                                                       cacgCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGAcgaccga                                                                           <  1:271107/146‑1 (MQ=255)
                                                                          gCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGt                                                                        >  1:41947/1‑149 (MQ=255)
                                                                              aTCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTTAGGCAGAAGTTATCAATTTTTTAACGCTAAAGCAGCGTGACGATACCGAGCTGGGTGCATGTACGCAGACGCGGAGCGCTCTGGCACAGACGACCGACGTTGTa                                                                    >  1:500576/1‑149 (MQ=255)
                                                                               tCGGGCGGCCCACGAAATTCGATGCGCCAAGCACCACGGCGTGGAGGCCGAAGGTATCAATGTTGGAACGCACACGCAGCGTGACGATACCGCGCGGGGGGCAGGGAGGCAGACGCGGCGCGCGCTGGCACAGACGACCGAAGTAGTaa                                                                   <  2:316806/149‑1 (MQ=255)
                                                                                   gcggcCAACGATATTCGAGGCGCCAATCACCCCGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCCAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGAt                                                               <  1:517811/149‑1 (MQ=255)
                                                                                    cggcCAACGATATTCGAGGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATg                                                              <  2:85526/149‑1 (MQ=255)
                                                                                       ccAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGaa                                                           >  1:775594/1‑149 (MQ=255)
                                                                                            gATATTCGATGCGCCAATCACCAAGGCGTTGAGGCCGAAGTTATCAAAGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTCTAAGGATGGAAACCgg                                                      >  2:685484/1‑148 (MQ=255)
                                                                                             atatTCGATGCGCCTATCCCCACCGCGTTGCGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTc                                                     <  1:941947/149‑1 (MQ=255)
                                                                                              tatTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTcc                                                    <  1:216626/149‑1 (MQ=255)
                                                                                                   gATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGCCGACCGCCGTTGTAAGGATGGAAACCGTCCACGTc                                               >  1:45554/1‑149 (MQ=255)
                                                                                                      gcgcCAATCCCCACGGCGTGGAGGCCGAAGGTAGCAATGCTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCttt                                            <  2:512816/149‑1 (MQ=255)
                                                                                                        gcCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACACTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGt                                          >  1:745329/1‑149 (MQ=255)
                                                                                                           aaTCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCg                                       >  1:968242/1‑149 (MQ=255)
                                                                                                             tCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGa                                     >  1:1011313/1‑149 (MQ=255)
                                                                                                                  acGTCGTTGAGTCCGAAGGTATCAATGGTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAAt                                <  2:84570/149‑1 (MQ=255)
                                                                                                                  acGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGGCTTTGTCCGGATGAAt                                <  1:88004/149‑1 (MQ=255)
                                                                                                                   cGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACCGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATa                               >  1:954027/1‑149 (MQ=255)
                                                                                                                    ggCGTTGAGGCCGAAGGTTTCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATAc                              <  2:405603/149‑1 (MQ=255)
                                                                                                                        ttGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTGTTTGTCCGGATGAATACGTTc                          >  1:245880/1‑149 (MQ=255)
                                                                                                                        ttGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTc                          <  2:881623/149‑1 (MQ=255)
                                                                                                                         tGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTcc                         <  2:451676/149‑1 (MQ=255)
                                                                                                                          gacgCCGCAGGTATCAATGGTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGTCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCa                        <  2:534787/146‑1 (MQ=255)
                                                                                                                          gAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCa                        <  1:388953/149‑1 (MQ=255)
                                                                                                                           aGGCCGAAGGTATCATTGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATTGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAg                       <  1:859305/149‑1 (MQ=255)
                                                                                                                             gCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCa                     >  2:484392/1‑149 (MQ=255)
                                                                                                                               cGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACt                   >  1:502277/1‑149 (MQ=255)
                                                                                                                                       taaaTGTTGTAACGCTCACGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGtt           <  2:112740/147‑1 (MQ=255)
                                                                                                                                       tCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGtt           <  1:856796/149‑1 (MQ=255)
                                                                                                                                        cAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTa          <  1:712505/149‑1 (MQ=255)
                                                                                                                                          aTGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGTGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATc        >  1:93465/1‑149 (MQ=255)
                                                                                                                                          aTGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAATGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATc        <  1:8189/149‑1 (MQ=255)
                                                                                                                                           tgttgtAACGCTCAAGCAGCGTGACGATACCGCGCGGGCTGCAGCGACGCAGACGCGGCGCGCGCGCGCACAGACGACCGCCGTTGAAAGGAAGGAAACCGTCCACGTCAATGACCGCATGAATACGTTCCAGCACTTTGCCGTAATCa       >  2:275425/1‑149 (MQ=255)
                                                                                                                                             ttgtAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGCCGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTGCAGCACTTTGCCGTTATCAAt     >  1:842257/1‑149 (MQ=255)
                                                                                                                                                tAACGCTCAAGCAGCGTGACGATACCTCGCGGGGTGCATGGACGCAGACGCGGCTCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATAcc  <  2:70671/149‑1 (MQ=255)
                                                                                                                                                tAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTGGTCCGGATGAATAAGTTCCAGCACTTTGACGTTATAAATAcc  >  2:838587/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TCTACGTGATGACGCAGATTTTTAGTGAAGCGGTGAGTCACTGTAGTGGTGCAACCTGCCAGCAGCAGTTCCATGCTCATCGGGCGGCCAACGATATTCGATGCGCCAATCACCACGGCGTTGAGGCCGAAGGTATCAATGTTGTAACGCTCAAGCAGCGTGACGATACCGCGCGGGGTGCAGGGACGCAGACGCGGCGCGCGCTGGCACAGACGACCGACGTTGTAAGGATGGAAACCGTCCACGTCTTTGTCCGGATGAATACGTTCCAGCACTTTGACGTTATCAATACC  >  NC_000913/557140‑557432

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑