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,3650CT89.5% 92.9 / 4.2 39K245K (AAG→AAAcysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (15/19);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.51e-01

GACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACG  >  NC_000913/2889124‑2889591
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    
gACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGttt                                                                                                                                                                                                                             >  1:70953/1‑251 (MQ=255)
   gctctccgtCCGCTTTAAACGTCCCACCCCCCACGCCCTCCAGCGGTTATTTGGTTTTGTCATTTTTCCGATCGGTTCGGTTCCCCCAGTCCCGCTGAGTTGTCTCGCCGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCg                                                                                                                                                                                                                          <  2:469461/243‑1 (MQ=255)
     tCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGc                                                                                                                                                                                                                        >  1:9880/1‑251 (MQ=255)
        ccTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGACCGGTTCGGTTACCCCAGTCACGCTGATTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGAt                                                                                                                                                                                                                     <  1:207641/250‑1 (MQ=255)
             cGCTTTAAACATCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGAGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGTTTACAGCCCCCCAGTTTGCCGTTTTCGTCGATCGccc                                                                                                                                                                                                                >  2:450015/1‑250 (MQ=255)
                         ctcAACCCCCACGCGCTCCAGCGTGTGTTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCTGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGt                                                                                                                                                                                                    >  1:121361/1‑251 (MQ=255)
                             ccccccACGCGCCCCAGCGTGTTTTGGGTTTTGTCATTTTTACGATCGGTTCGGTCACCCCAGTCACGCGGAGTTTTCACGCCGGCTCGGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGtt                                                                                                                                                                                                <  2:156783/250‑1 (MQ=255)
                                   aCGCGCCCCAGCGTTTATTGGTTTTTGGCATTTTTACGATCGGTCCGGTTACCCCAGTCACCCTGGGTTGTCACGACGGCACCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACGCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTg                                                                                                                                                                                          <  2:344993/251‑1 (MQ=255)
                                         tCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGACGTAGGATTTCTTGTTGCCGTGTTCGATGGAAAGACCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCa                                                                                                                                                                                                                >  1:171912/1‑223 (MQ=255)
                                            gggtttttttttgTTTTTGCATTTCTACGATCGGTCATGTTACCCCATTCGCGTTGCGTTGTCTCGGCGGCTGCGGCCGCCGCCAGCTTATGCTCCAGCGGCAGATAGCCAAACTCAATCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTGTCGGCGATCGCCACGAAGTTCATGTCGTGGGCGTGCAGATCGat                                                                                                                                                                                 <  2:740/239‑1 (MQ=255)
                                            aGCGTGTATTTGGTTGTGGCTTTTTTCCGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGat                                                                                                                                                                                 <  1:102833/251‑1 (MQ=255)
                                             gCGTGTATTTGGTTGTGGCATTTTTACGAGCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGata                                                                                                                                                                                <  1:442726/251‑1 (MQ=255)
                                              cGTGTATTTGGTTTTGGAATTTTTACGATCGGCTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGGAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGGAGTTACTGTCCTTGGCGTGCAGCTCGatat                                                                                                                                                                               >  2:431080/1‑251 (MQ=255)
                                                        gTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTgg                                                                                                                                                                     >  1:11444/1‑251 (MQ=255)
                                                          tttGGCATTTTTCCGATCGGTTGGGTCACCCCAGTCACGCTGATTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGATCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGc                                                                                                                                                                    <  1:124902/250‑1 (MQ=255)
                                                           ttGGCATTTTTACGACCGGTTCGGTTACCCCTGTCACGTGGAGTTTACACGCCGCCTTCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCCCTCTCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACAAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCg                                                                                                                                                                   <  1:450617/250‑1 (MQ=255)
                                                                 tttttACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCTATCGCCCCGAAGTTCATGTCGTTGGCGTGCTGATCGGTATCGTTCTGTGGCggtggcgc                                                                                                                                                            >  2:186814/1‑245 (MQ=255)
                                                                  ttttACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCggggtgcac                                                                                                                                                           >  2:236330/1‑245 (MQ=255)
                                                                             gttcggttACCCACGTCACGCTGGGTTGTCACGCCGGCTTCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCGCACCAGTTTGCCGTTTGCGGCGATCGCCACGAAGTTCATTTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGtttt                                                                                                                                              <  1:469759/251‑1 (MQ=255)
                                                                                  gttACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCTGGATCACTACCGTGGTTTTGAAtt                                                                                                                                         >  1:297259/1‑251 (MQ=255)
                                                                                                    aGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATCTCGTTCTGTGGCGGGGTCCCTACCCGTGGttttcatttttatggggggtgggggt                                                                                                                        >  2:220330/1‑229 (MQ=255)
                                                                                                        cccgatctgtcctCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGcc                                                                                                                   <  2:272591/239‑1 (MQ=255)
                                                                                                          ccccacggttccgcccACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCg                                                                                                                  <  2:287700/236‑1 (MQ=255)
                                                                                                             gacgGCTTAGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGACGTTCATGTCGTTGGCGTGCCGATCGATATCGTTCTGTGGCGGGGGCACCACACGTGGTTTGTATTTACGCGGCaggaagagggggccgcggg                                                                                                              >  2:466874/1‑235 (MQ=255)
                                                                                                             gacgGCTCCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGa                                                                                                              <  1:220330/251‑1 (MQ=255)
                                                                                                                ggCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATTACTACCGGGGTTTTGGATTTTAGCGGCAGGTGGGTGTGGcccggtcggg                                                                                                           >  1:272591/1‑243 (MQ=255)
                                                                                                                ggCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAACCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTTCCGTGGGTTTGGATTTTCGCGGCAGGTAGGGCTGGCCgggggtcg                                                                                                           >  1:197783/1‑244 (MQ=255)
                                                                                                                        ccACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGAGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGATAGGTCTGGCCGAGGATCGGTTCTtca                                                                                                   <  1:58784/251‑1 (MQ=255)
                                                                                                                          aCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGCAGTTCATGTCGTTGGc                                                                                                                                                                                             >  2:210462/1‑161 (MQ=255)
                                                                                                                          aCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGc                                                                                                                                                                                             <  1:210462/161‑1 (MQ=255)
                                                                                                                               cAGCGTATGCTCCAGCGGCAGATAGCCAAACGCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGGTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCagtag                                                                                            <  2:407372/251‑1 (MQ=255)
                                                                                                                                   gtaggCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTtcatca                                                                                                <  2:66703/239‑1 (MQ=255)
                                                                                                                                   gTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGGTCACTACCGTGGGTTTTAATTTTCGCGGCCGGGAGGGCTGGCCGGGGGTCGGTTCTtcattaagatcggg                                                                                        >  1:66703/1‑241 (MQ=255)
                                                                                                                                                 cAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCTGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGt                                                                          <  2:217540/251‑1 (MQ=255)
                                                                                                                                                  aGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTc                                                                         <  2:9880/251‑1 (MQ=255)
                                                                                                                                                         cAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTAACGGGGGTTTGAATTTACGCGGCAGGGAGGGCTGGCCGAGGGTCGGTTCTTCCTCCGTAGGGGCGACTTTTTCCTgggcgtgcccg                                                                  >  1:470736/1‑242 (MQ=255)
                                                                                                                                                           aacgcactcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACAGTGGTGTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGAt                                                                <  1:279225/247‑1 (MQ=255)
                                                                                                                                                                     cggttccatcgtcaggtttCTTTTTGCCGTGTTCTAAGGAAAGCCCACCCACCATCAACAGGTTAAAGCCCACCAGGTTGCCGTTGTCGGCGCTCGCCCCGAAGTTTATGTCGTGTGCGTGCAGATCGTGACCGTTCTGGGGCTGTATCACTACAGTGGTTTTGAATTTACGCGGCAGGTATTTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGATCGATCTTGATCTCCGTATAc                                                      <  2:315483/235‑1 (MQ=255)
                                                                                                                                                                                                                         ccccaACAGGTTAAAGCCCTCCATTTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTGTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGATGCAACAGATTCTCAGAAATCTTCTTCGCCCACTCGTACg  <  1:190433/249‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                    
GACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGACGGCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCACTCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGAAGTTCATGTCGTTGGCGTGCAGATCGATATCGTTCTGTGGCGGGATCACTACCGTGGTTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACG  >  NC_000913/2889124‑2889591

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