Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,889,365 C→T K245K (AAG→AAA cysI ← sulfite reductase, hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,3650CT97.4% 132.4 / ‑2.6 38coding (735/1713 nt)cysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (14/23);  minor base . (0/1);  total (14/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.34e-01

GACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCA  >  NC_000913/2889124‑2889603
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
gACGCTCCACTTCCGCTTTAAACGTCTCCACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGttt                                                                                                                                                                                                                                         >  2:224134/1‑251 (MQ=255)
        aCTTCCGCTTTAAACGTCTCCACCCCCACGCGCTCCAGCGTGTATTGGGTTTTGGAATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGGGTTGTACCGACGGCGTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGAt                                                                                                                                                                                                                                 <  2:275828/251‑1 (MQ=255)
              gCTTTAACCGTCTCACCCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTAAGCTCGGTTCGGTTACCCCAGTCACGCTGGGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCAc                                                                                                                                                                                                                           <  2:156723/251‑1 (MQ=255)
                        ccgcatccaccACTCGCTCCAGCTTGTATTGGGTGTTGGCATTTTTACGATCGGTGCGGTTCCCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTGTCGGCGATCGCCACGAAGTTCATg                                                                                                                                                                                                                 <  2:355850/242‑1 (MQ=255)
                               ccccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAACGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTgg                                                                                                                                                                                                          >  2:130893/1‑251 (MQ=255)
                                 ccACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCg                                                                                                                                                                                                        >  1:295481/1‑251 (MQ=255)
                                            aGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGat                                                                                                                                                                                             >  2:5143/1‑251 (MQ=255)
                                             gCGTGTATTTGGTTTTGGCATTTTTCCGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTGGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                                                                                                                                                                                            <  1:210602/251‑1 (MQ=255)
                                                   aTTTGGTTTTGGCATTTTTACGCTCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTc                                                                                                                                                                                      >  1:270037/1‑251 (MQ=255)
                                                   aTTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGtt                                                                                                                                                                                       <  1:5143/250‑1 (MQ=255)
                                                                 tttttACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCggggtacc                                                                                                                                                                        >  1:285765/1‑246 (MQ=255)
                                                                         aTCGGTTCGGATACCCCAGTCAAGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCAACGCAGTTCATGTCGTTGGCGTGCAGATCGATCTCGTTCTGTGGCggggggacgatcgggg                                                                                                                                                                >  2:161250/1‑238 (MQ=255)
                                                                            gtagcgtttccccgaGTCGCGCTGAGTTGGCACGCCGCCTTGGTCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACGCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGGGTTCGATGGAAAGCCCACCGCCCCCCAACAGGTTAAAGCCCACCAGTTTGCCTTTTTCGGCGATCGCCTCGAAGGGCATGTCGTTGGCGTGCAGATGGATATCTTTCTGTGGCGGGATCACTACCGTGGttg                                                                                                                                                            <  2:353793/237‑2 (MQ=255)
                                                                               tcggtgACCCCAGTCTCGCGGAGTGGTCACGTCGGCTTCGGCCGCCGGCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCGCCGCCCATGAGCAGGTTATAGCCCACCTGTTGTCCCTTTTCGGCGATCGCCACGAAGTTCATGTCGTGGGCGTGCGTATCGATATCTTTCTGTGGCGGGATCACTACCGTGGTTTTg                                                                                                                                                          <  1:346585/250‑1 (MQ=255)
                                                                                     cccccAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTa                                                                                                                                                    <  2:113744/249‑1 (MQ=255)
                                                                                           gTCACGCTGAGTTGTCCCGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCa                                                                                                                                             <  2:79187/251‑1 (MQ=255)
                                                                                                        gTCACGACGGCTTCGGCCACCGCCAGCGTATGCGCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTGTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGAGATCGGTCTGTGGCGGGATCACTACCGTTGTTTTGAATTTACGCGGCAGGTAGGTCTGGc                                                                                                                                <  1:221734/250‑1 (MQ=255)
                                                                                                                ggCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTTCCGTGGGTTTGAATTTACGCGGCGGGTGGGTCTgggcggggggcg                                                                                                                       >  2:159570/1‑241 (MQ=255)
                                                                                                                ggCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCg                                                                                                                       <  1:345947/251‑1 (MQ=255)
                                                                                                                 gCTTCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATTTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGAt                                                                                                                                                                           <  2:241218/200‑1 (MQ=255)
                                                                                                                            cGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCag                                                                                                           <  1:194515/251‑1 (MQ=255)
                                                                                                                            cGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCag                                                                                                           <  1:232613/251‑1 (MQ=255)
                                                                                                                              ccAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACACTCGCGGGTTTGAATTTATGCGGCGGGTGGGGCTGGcggcgggtagggtttttatcact                                                                                                          >  2:203855/1‑228 (MQ=255)
                                                                                                                                       gCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGttc                                                                                                                   <  1:96133/232‑1 (MQ=255)
                                                                                                                                                  aGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTGTATTTATGCGGCAGGTCGGGCTGGCCGAGGAGCGGGTCTTCTTCTGTAGTGGCgccctttttcttgtc                                                                                     >  2:354939/1‑237 (MQ=255)
                                                                                                                                                   gATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCg                                                                                    <  2:231466/251‑1 (MQ=255)
                                                                                                                                                   gATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACGACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGAACGGTTCTTCATCAGGAGTGGCGACTTTTTCCTGGTCg                                                                                    <  2:160028/251‑1 (MQ=255)
                                                                                                                                                             ctcactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTGAATTTACGCGGCAGGGAGGGCTGGCCGAGGGTCGGTTGTTCATCAGTAGTGGCGGCttttttcttgttggggcaggttt                                                                          >  1:99448/1‑233 (MQ=255)
                                                                                                                                                                     cgGTGCGGGCGGAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCCCCGCCCACCAACAGGTTAAAGCCCACCAGTTTGACGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGTCTTTTTCCTGGTCGAGCCAGATCTCCGCATAc                                                                  <  1:219369/251‑1 (MQ=255)
                                                                                                                                                                               gTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTAc                                                        <  1:130893/251‑1 (MQ=255)
                                                                                                                                                                                     tttCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTTGTTTTTTATTTTCGCGGCCGGGAGGGCTGGCCGAGGGTCGGTTCTTCATCCGTAGGTGCGACTTTTTCCTTGTCGAGCCAGATCTCCGCCTACGCGCGGGTACGAggg                                                   >  1:73633/1‑249 (MQ=255)
                                                                                                                                                                                      ttCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGACGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGGCACTACCGGGGTTTTTAATTTACGCGGCAGGGGAGGCTGGGCGCGGGACGGTTCTTCATTCGTAGTGGCGACTTTTTTCTGGGCGAGCCAGAACTCCGCATCCCCCCGGGGACGAGGCa                                                  >  2:141952/1‑250 (MQ=255)
                                                                                                                                                                                                cGTGTTCGATGGAAAGCCCACCGCCCCCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCa                                       <  2:105236/251‑1 (MQ=255)
                                                                                                                                                                                                   gtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGaa                                    <  1:349223/251‑1 (MQ=255)
                                                                                                                                                                                                        aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACt                                                                                                                                                                       <  2:120055/117‑1 (MQ=255)
                                                                                                                                                                                                        aTGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACt                                                                                                                                                                       >  1:120055/1‑117 (MQ=255)
                                                                                                                                                                                                              aGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGc                         <  1:297276/251‑1 (MQ=255)
                                                                                                                                                                                                                                     aaaGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCa  <  1:159570/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                 
GACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCA  >  NC_000913/2889124‑2889603

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