Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,381,396 G→C 5.7% P308R (CCG→CGG)  frdA ← anaerobic fumarate reductase catalytic and NAD/flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,381,3960GC 5.7% 484.0 / 15.0 139P308R (CCG→CGG) frdAanaerobic fumarate reductase catalytic and NAD/flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (56/75);  new base C (3/5);  total (59/80)
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.69e-01

Only 100 of 140 total aligned reads displayed.

GTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTA  >  NC_000913/4381248‑4381538
                                                                                                                                                    |                                                                                                                                              
gtAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCAcgcg                                                                                                                                                >  2:954049/1‑149 (MQ=255)
       gCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTgg                                                                                                                                          <  4:131131/148‑1 (MQ=255)
       gCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTgg                                                                                                                                          <  4:232204/148‑1 (MQ=255)
         ggTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGgaga                                                                                                                                       <  1:483235/149‑1 (MQ=255)
          gTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCTCCACGCCGCGTGGCGAt                                                                                                                                      <  1:596238/149‑1 (MQ=255)
          gTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGgaga                                                                                                                                       <  3:551520/148‑1 (MQ=255)
          gTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGgaga                                                                                                                                       <  3:115099/148‑1 (MQ=255)
           tCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACTCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGAt                                                                                                                                      <  4:613585/148‑1 (MQ=255)
           tCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCTCCACGCCGCGTGGCGAt                                                                                                                                      <  4:466108/148‑1 (MQ=255)
            cGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCTCCACGCCGCGTGGCGATgt                                                                                                                                    <  2:495408/149‑2 (MQ=255)
             ggACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGgt                                                                                                                                   <  1:738377/149‑1 (MQ=255)
              gACGTACCGGAATCGGTTCTTTAACCGGACCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACACCATCTCCACGCCGCGTGGCGATgtg                                                                                                                                   <  4:789852/148‑3 (MQ=255)
                cGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGtt                                                                                                                                <  1:954049/149‑1 (MQ=255)
                              ttCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTg                                                                                                                               <  4:588424/136‑1 (MQ=255)
                                 tttAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGt                                                                                                               >  1:734226/1‑149 (MQ=255)
                                  ttAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGt                                                                                                               >  3:235292/1‑148 (MQ=255)
                                  ttAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGt                                                                                                               >  3:88027/1‑148 (MQ=255)
                                   tAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGc                                                                                                             <  2:561398/149‑1 (MQ=255)
                                    aaCCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGc                                                                                                             <  3:925719/148‑1 (MQ=255)
                                      ccGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAg                                                                                                          >  2:454801/1‑149 (MQ=255)
                                       cGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCa                                                                                                                    >  3:1004563/1‑138 (MQ=255)
                                       cGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCa                                                                                                                    <  4:1004563/138‑1 (MQ=255)
                                       cGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAg                                                                                                          >  3:21303/1‑148 (MQ=255)
                                         gATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGaa                                                                                                        <  4:857734/148‑1 (MQ=255)
                                           tCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAgg                                                                                                      >  3:938696/1‑148 (MQ=255)
                                           tCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGc                                                                                                     >  2:303354/1‑149 (MQ=255)
                                               cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGTGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTga                                                                                                 <  2:185479/149‑1 (MQ=255)
                                               cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTga                                                                                                 <  2:1049827/149‑1 (MQ=255)
                                               cGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTg                                                                                                  <  4:692685/148‑1 (MQ=255)
                                                 ccAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgag                                                                                                <  4:626555/148‑1 (MQ=255)
                                                    aCGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgagaga                                                                                             >  4:937437/1‑148 (MQ=255)
                                                    aCGTACGCTGTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGGCGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgagaga                                                                                             >  4:562986/1‑148 (MQ=255)
                                                        aCGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTgaga                                                                                               <  1:471248/142‑1 (MQ=255)
                                                             ttCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTcgc                                                                                    >  3:348779/1‑148 (MQ=255)
                                                              tCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGt                                                                                  <  2:827614/149‑1 (MQ=255)
                                                                gCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTgg                                                                                >  2:930885/1‑149 (MQ=255)
                                                                gCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTgg                                                                                >  1:995623/1‑149 (MQ=255)
                                                                  cAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCAtt                                                                                                                  <  2:294690/113‑1 (MQ=255)
                                                                  cAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCAtt                                                                                                                  >  1:294690/1‑113 (MQ=255)
                                                                    gTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGAcc                                                                             <  3:562986/148‑1 (MQ=255)
                                                                        gCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGtt                                                                        <  2:734226/149‑1 (MQ=255)
                                                                        ctgcaTGAACGTCTGCCGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGt                                                                         >  4:767232/5‑148 (MQ=255)
                                                                          aGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTcc                                                                      <  2:69080/149‑1 (MQ=255)
                                                                          aGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTc                                                                       <  4:938696/148‑1 (MQ=255)
                                                                             tGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCata                                                                   >  2:471296/1‑149 (MQ=255)
                                                                              gAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTga                                                                                                             >  1:756828/1‑105 (MQ=255)
                                                                              gAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTga                                                                                                             <  2:756828/106‑2 (MQ=255)
                                                                               aaCGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTcgcg                                                                                   >  3:597698/1‑131 (MQ=255)
                                                                                  ggCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATAtt                                                               <  3:507014/148‑1 (MQ=255)
                                                                                   gCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATAttt                                                              >  3:796989/1‑148 (MQ=255)
                                                                                   gCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTg                                                             >  1:787888/1‑149 (MQ=255)
                                                                                    cAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTg                                                             <  3:411343/148‑1 (MQ=255)
                                                                                     aGACGTTCATACAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACGTTGTCGCGTGGACCCAGTTCCATATATTTGt                                                            <  4:676122/148‑1 (MQ=255)
                                                                                       aCGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGttt                                                          <  4:753040/148‑1 (MQ=255)
                                                                                         gTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGttttt                                                        <  3:1023683/148‑1 (MQ=255)
                                                                                          ttCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTc                                                       >  3:685251/1‑148 (MQ=255)
                                                                                           tCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTcgcg                                                                                   <  3:843499/119‑1 (MQ=255)
                                                                                           tCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTcgcg                                                                                   >  4:843499/1‑119 (MQ=255)
                                                                                              tGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCt                                                                                                   <  1:116041/100‑1 (MQ=255)
                                                                                                  gTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGccc                                              >  2:342123/1‑149 (MQ=255)
                                                                                                  gTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGcc                                               >  3:651275/1‑148 (MQ=255)
                                                                                                       ttCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCgg                                         <  2:787888/149‑1 (MQ=255)
                                                                                                        tctcGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCgg                                         <  3:838162/148‑1 (MQ=255)
                                                                                                        tctcGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCgg                                         <  4:651275/148‑1 (MQ=255)
                                                                                                        tctcGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGa                                        <  2:331669/149‑1 (MQ=255)
                                                                                                           cGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGt                                      >  4:123437/1‑148 (MQ=255)
                                                                                                               gAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCg                                 <  1:271346/149‑1 (MQ=255)
                                                                                                                aGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCgg                                <  1:454801/149‑1 (MQ=255)
                                                                                                                aGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCg                                 <  4:318786/148‑1 (MQ=255)
                                                                                                                  gTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCggg                               <  4:348779/148‑1 (MQ=255)
                                                                                                                  gTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGc                              >  2:182074/1‑149 (MQ=255)
                                                                                                                     aCGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCa                           >  2:1026986/1‑149 (MQ=255)
                                                                                                                      cGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCCTATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCa                           >  4:838670/1‑148 (MQ=255)
                                                                                                                      cGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGCTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCa                           <  4:88027/148‑1 (MQ=255)
                                                                                                                       gCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTga                                                                                                 >  1:43095/1‑77 (MQ=255)
                                                                                                                       gCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATg                         >  1:521930/1‑149 (MQ=255)
                                                                                                                        cAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATg                         >  3:152217/1‑148 (MQ=255)
                                                                                                                        cAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTCTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATg                         <  4:21303/148‑1 (MQ=255)
                                                                                                                           gTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGa                                        <  1:299695/130‑1 (MQ=255)
                                                                                                                            tCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGt                     >  4:1003395/1‑148 (MQ=255)
                                                                                                                            tCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTa                    >  1:55131/1‑149 (MQ=255)
                                                                                                                             cGAGATAACACACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTaa                   <  2:84032/149‑1 (MQ=255)
                                                                                                                               agaTAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAAt                  <  3:24796/148‑1 (MQ=255)
                                                                                                                               agaTAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATc                 >  1:976335/1‑149 (MQ=255)
                                                                                                                                 aTAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGGCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCtt               <  1:187642/149‑1 (MQ=255)
                                                                                                                                 aTAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCt                <  3:937437/148‑1 (MQ=255)
                                                                                                                                  tAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTg                                                                                                  <  4:88143/65‑1 (MQ=255)
                                                                                                                                  tAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTg                                                                                                  >  3:88143/1‑65 (MQ=255)
                                                                                                                                    aacaccaTCTCCACGCCGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCa            >  1:139032/6‑149 (MQ=255)
                                                                                                                                     aCCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAg           >  1:198262/1‑149 (MQ=255)
                                                                                                                                      ccACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAg           >  3:855190/1‑148 (MQ=255)
                                                                                                                                       cacaTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGa          >  3:773998/1‑148 (MQ=255)
                                                                                                                                       cacaTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGAt         >  2:467838/1‑149 (MQ=255)
                                                                                                                                       accaTCTCCACGCCGCGTGGCGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGAt         >  2:234450/3‑149 (MQ=255)
                                                                                                                                        acaTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGTGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATa        >  2:703530/1‑149 (MQ=255)
                                                                                                                                        acaTCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGAt         >  3:195844/1‑148 (MQ=255)
                                                                                                                                            cGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACgg    <  1:471296/149‑1 (MQ=255)
                                                                                                                                            cGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACgg    <  1:1020605/149‑1 (MQ=255)
                                                                                                                                              ccACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGt   <  4:152217/148‑1 (MQ=255)
                                                                                                                                              ccACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTa  <  2:137295/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
GTAGTGTGCGGTCGGACGTACCGGAATCGGTTCTTTAACCGGATCGACGCCAACGTACGCTTTCGCCAGTTCGCAGATGAACGGCAGACGTTCATGCAGTTTTTTCTCGCCGAGGTGACGCAAGTCGAGATAAACCACATCGCCACGCGGCGTGGAGATGGTGTTGCCTTTACGCCATTCGTGCCAGAAGGCCTGAGAGACTTTGTCGCGTGGACCCAGTTCCATATATTTGTTTTTCGGCTCGCCCAGCGGAGTTTCCGGGCCCATGCCGTAATCTTGCAGATAACGGTA  >  NC_000913/4381248‑4381538

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: