Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,602,017 (C)7→8 coding (200/2019 nt) hyfB → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,602,0101.C86.7% 550.2 / 85.6 180coding (193/2019 nt)hyfBhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base . (13/11);  new base C (67/89);  total (80/100)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 202 total aligned reads displayed.

ATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTC  >  NC_000913/2601911‑2602104
                                                                                                    |                                                                                                
aTTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTg                                                                                                   >  2:1969269/1‑100 (MQ=255)
aTTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTg                                                                                                   >  1:2025059/1‑100 (MQ=255)
aTTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTg                                                                                                   >  1:2523127/1‑100 (MQ=255)
 ttAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑c                                                                                                <  1:2309784/100‑1 (MQ=255)
  tAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑c‑c                                                                                              <  1:63657/100‑1 (MQ=255)
  tAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTA‑‑c‑c                                                                                              <  1:5389774/100‑1 (MQ=255)
   aagttttCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑c‑cc                                                                                             <  2:2339430/96‑1 (MQ=255)
   aaGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑c‑cc                                                                                             >  2:2582244/1‑100 (MQ=255)
     gCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑c‑cccc                                                                                           >  2:4955702/1‑100 (MQ=255)
     gCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑c‑cccc                                                                                           >  2:385116/1‑100 (MQ=255)
      cTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑c‑ccccc                                                                                          <  2:2000980/100‑1 (MQ=255)
       tttCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑c‑cccccc                                                                                         >  1:641517/1‑100 (MQ=255)
          ccGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCctt                                                                                      <  1:5511794/100‑1 (MQ=255)
           cGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCcttt                                                                                     >  1:3370032/1‑100 (MQ=255)
            ggCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCctttt                                                                                    >  2:699280/1‑99 (MQ=255)
            ggCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCctttt                                                                                    <  1:786927/100‑2 (MQ=255)
            ggCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCctttt                                                                                    <  2:275731/100‑2 (MQ=255)
             gCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTg                                                                                    >  1:2987685/1‑100 (MQ=255)
             gCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCcttttg                                                                                   <  2:4743679/100‑3 (MQ=255)
             gCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTGTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTg                                                                                    <  1:4896519/100‑1 (MQ=255)
             gCATCACATCACTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTg                                                                                    <  1:2893183/100‑1 (MQ=255)
               aTCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGaa                                                                                  >  1:350057/1‑100 (MQ=255)
                 cacaTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAAtt                                                                                >  1:2932117/1‑100 (MQ=255)
                  acaTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAAttt                                                                               >  1:4158249/1‑100 (MQ=255)
                   caTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTg                                                                              <  1:5280515/100‑1 (MQ=255)
                   caTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTg                                                                              <  1:4368885/100‑1 (MQ=255)
                   caTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTg                                                                              <  2:3611676/100‑1 (MQ=255)
                   caTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTg                                                                              >  1:2669031/1‑100 (MQ=255)
                   caTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCTTTTGAATTTGc                                                                             >  1:1600646/1‑100 (MQ=255)
                    aTCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGc                                                                             <  2:4433343/100‑1 (MQ=255)
                     ccGCTGGTGGGCGGCGTGTTTGGCATCCTCAGCGGAATTGCGGCATTACATGCTCGTGTATCTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGcc                                                                            <  2:5057650/99‑1 (MQ=255)
                      ngCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCg                                                                           >  1:907609/2‑100 (MQ=255)
                      cGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCg                                                                           >  1:3301887/1‑100 (MQ=255)
                        cTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGAt                                                                         <  2:2003281/100‑1 (MQ=255)
                         tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGAtt                                                                        >  1:364558/1‑100 (MQ=255)
                         tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGAtt                                                                        <  2:1724030/100‑1 (MQ=255)
                         tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGAtt                                                                        <  1:707989/100‑1 (MQ=255)
                         tggtggGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGAtt                                                                        <  2:4146261/100‑1 (MQ=255)
                             gggCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAAc                                                                    >  2:5394850/1‑100 (MQ=255)
                             gggCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAAc                                                                    <  2:903029/100‑1 (MQ=255)
                             gggCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAAc                                                                    <  1:2039681/100‑1 (MQ=255)
                                   cgTGTTTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                                              <  1:4990768/100‑1 (MQ=255)
                                   cgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                                              >  2:5025/1‑100 (MQ=255)
                                   cgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                                              <  2:2959693/100‑1 (MQ=255)
                                   cgTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCg                                                              <  1:873626/100‑1 (MQ=255)
                                      gATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGAGAAt                                                           <  1:2233116/100‑1 (MQ=255)
                                      gATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACCTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAAt                                                           >  2:3161373/1‑100 (MQ=255)
                                       aTTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATg                                                          <  2:3806400/100‑1 (MQ=255)
                                         tGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGAt                                                       <  2:4884410/100‑1 (MQ=255)
                                         tGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACCTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGa                                                        <  1:3190884/100‑1 (MQ=255)
                                          ggCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGAt                                                       >  1:4755935/1‑100 (MQ=255)
                                          ggCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGAt                                                       <  1:4226181/100‑1 (MQ=255)
                                          ggCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGAt                                                       >  1:998125/1‑100 (MQ=255)
                                           gCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATa                                                      >  2:4552878/1‑100 (MQ=255)
                                           gCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATa                                                      <  2:3882082/100‑1 (MQ=255)
                                               catcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCct                                                  >  1:1852245/1‑100 (MQ=255)
                                                atcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCctc                                                 >  1:1183802/1‑100 (MQ=255)
                                                 tcaGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCtctc                                               >  2:3619203/1‑100 (MQ=255)
                                                   aGCGGAATTACGCCATGACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTATCCCTGCGAATGGATAGCCTCTCg                                              <  2:4741341/100‑1 (MQ=255)
                                                   aGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCg                                              <  1:1053103/100‑1 (MQ=255)
                                                   aGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCg                                              <  1:1656835/100‑1 (MQ=255)
                                                   aGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCg                                              >  2:1312632/1‑100 (MQ=255)
                                                    gCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCgg                                             <  1:3128928/100‑1 (MQ=255)
                                                    gCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCgg                                             <  1:5059/100‑1 (MQ=255)
                                                       gAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCAt                                          >  1:5206345/1‑100 (MQ=255)
                                                        aaTTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCAtt                                         <  2:2796646/100‑1 (MQ=255)
                                                         aTTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCAttt                                        >  1:1687297/1‑100 (MQ=255)
                                                           tACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTAt                                      >  2:2244493/1‑100 (MQ=255)
                                                           tACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCTTTTGGATTTGCCGATTTAACCCTGCGAATGGATAGCCCCTCGGCATTTATg                                     <  1:4063900/100‑1 (MQ=255)
                                                            aCGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATg                                     >  2:3061259/1‑100 (MQ=255)
                                                            aCGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATg                                     <  2:1231685/100‑1 (MQ=255)
                                                             cGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATgg                                    <  2:4670312/100‑1 (MQ=255)
                                                             cGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATgg                                    <  1:3659980/100‑1 (MQ=255)
                                                               cAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTg                                  >  1:5073764/1‑100 (MQ=255)
                                                                 aTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCt                                >  2:4044803/1‑100 (MQ=255)
                                                                 aTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCt                                <  2:1494462/100‑1 (MQ=255)
                                                                 aTTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑CACCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTgat                                >  2:3578625/1‑98 (MQ=255)
                                                                    aCATGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGt                             >  2:4812664/1‑100 (MQ=255)
                                                                      aTGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTa                           <  1:686955/100‑1 (MQ=255)
                                                                      aTGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTa                           <  1:37736/100‑1 (MQ=255)
                                                                      aTGCTGGTGTAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTa                           <  2:4401949/100‑1 (MQ=255)
                                                                       tGCTGGTGTAACTTTAGTCGCCCGTTTTGCCC‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTa                           <  2:157629/100‑1 (MQ=255)
                                                                            gtgtAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCt                     >  2:1864938/1‑100 (MQ=255)
                                                                              gtAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCCTg                   >  1:1992739/1‑100 (MQ=255)
                                                                               tAACTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGc                  <  1:2265348/100‑1 (MQ=255)
                                                                                aaCTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCt                 >  1:1407872/1‑100 (MQ=255)
                                                                                 aCTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGAATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtg                >  2:5239970/1‑100 (MQ=255)
                                                                                  cTTTAGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtgg               >  2:4007510/1‑100 (MQ=255)
                                                                                      aGTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtggt           <  1:1131979/100‑1 (MQ=255)
                                                                                       gTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtggtg          >  2:3284203/1‑100 (MQ=255)
                                                                                       gTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggtggtg          >  1:5371100/1‑100 (MQ=255)
                                                                                       gTCGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCtggcggcg          >  1:2193662/1‑98 (MQ=255)
                                                                                         cGCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGTCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTggtggtggt        <  1:4538208/100‑1 (MQ=255)
                                                                                          gCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTTCCGATTTAACCCTTCGGATGGATAGCCTCTCGGCCTTTATTGTGCTGGTTATTTCCTTGCTGGTGGTGGtt       >  1:1238169/1‑100 (MQ=255)
                                                                                          gCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGtt       >  1:1875856/1‑100 (MQ=255)
                                                                                          gCCCGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGtt       <  2:3625638/100‑1 (MQ=255)
                                                                                           cccGTTTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGttt      >  1:2897287/1‑100 (MQ=255)
                                                                                           cccGATTTGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCCGCATTTATGGTGCTGGTTATCTCCCTGCTGGTGGTGGttt      >  2:778633/1‑100 (MQ=255)
                                                                                            cctttttGC‑C‑CCCCCCTTTTTAACTTGCCGATTTAGCCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGCTTTg     <  2:4646042/97‑1 (MQ=255)
                                                                                               ttttGC‑C‑CCCCCCTTTTGAATTTGCCGATTTAGCCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTtc  <  2:5496352/100‑1 (MQ=255)
                                                                                                    |                                                                                                
ATTAAGCTTTCCGGCATCACATCGCTGGTGGGCGGCGTGATTGGCATCATCAGCGGAATTACGCAATTACATGCTGGTGTAACTTTAGTCGCCCGTTTTG‑‑C‑CCCCCCTTTTGAATTTGCCGATTTAACCCTGCGAATGGATAGCCTCTCGGCATTTATGGTGCTGGTTATCTCCTTGCTGGTGGTGGTTTGTTC  >  NC_000913/2601911‑2602104

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