Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,890,626 G→A L425F (CTC→TTC)  cysJ ← sulfite reductase, flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,890,6260GA98.1% 158.6 / ‑4.4 52L425F (CTC→TTC) cysJsulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (27/24);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

TAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCG  >  NC_000913/2890385‑2890861
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
tAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTggg                                                                                                                                                                                                                                     >  1:164290/1‑250 (MQ=255)
 aGCGCCGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGa                                                                                                                                                                                                                                     <  2:164290/250‑1 (MQ=255)
             ccACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCa                                                                                                                                                                                                                        <  1:407865/251‑1 (MQ=255)
                       cAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGc                                                                                                                                                                                                              >  1:376446/1‑251 (MQ=255)
                        aGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGcc                                                                                                                                                                                                             >  1:426733/1‑251 (MQ=255)
                            aGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTcc                                                                                                                                                                                                         <  1:363124/251‑1 (MQ=255)
                             gTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGGTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTccc                                                                                                                                                                                                        >  2:457138/1‑250 (MQ=255)
                              tCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAc                                                                                                                                                                                                       >  1:440788/1‑251 (MQ=255)
                                     tAAAGTGCGGATTACCAAAGAACAGCCAGATTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACAGCGCAATACCGGTGCCTTGGCCAATCATAATCACCGGGGTTTCTGGATTGTCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGCACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGCGGCACCACCGGCACGGGCGCGGCCTTCCAcgacttac                                                                                                                                                                                                >  1:139246/1‑245 (MQ=255)
                                          ggctgATTCCCAAAGAACAGCCAGTTTTTACCGGGCGCTTCGGCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcac                                                                                                                                                                                           <  2:171579/247‑1 (MQ=255)
                                            cGGATGACCAAAGAACAGTCAGTGGTGTCCTGGCGCTTCGGCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTGCTTCCACGCGGGCAGCGAAGAAGCTGGAGGTACCCCCGGCACTGGCGCGGCCTTCCACGTCGTAACGcacca                                                                                                                                                                                         <  2:472014/251‑1 (MQ=255)
                                               aTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGGCGTAACGcaccacac                                                                                                                                                                                      >  1:233278/1‑251 (MQ=255)
                                                     aaaGAACAGCCAGTGTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTGCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACgg                                                                                                                                                                               <  1:215038/251‑1 (MQ=255)
                                                       aGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACgg                                                                                                                                                                               >  1:234531/1‑249 (MQ=255)
                                                         aaCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGAc                                                                                                                                                                           <  2:118672/251‑1 (MQ=255)
                                                          aCAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTCCATAAAGGCGCGGAACGGCGAAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACg                                                                                                                                                                          <  1:174979/251‑1 (MQ=255)
                                                               cAGTTTTTCCCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATCCCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTAc                                                                                                                                                                     <  2:376764/251‑1 (MQ=255)
                                                               cAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTAc                                                                                                                                                                     >  1:282411/1‑251 (MQ=255)
                                                                 gTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACtt                                                                                                                                                                   >  1:365474/1‑251 (MQ=255)
                                                                        ccTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtctc                                                                                                                                                            <  1:147851/251‑1 (MQ=255)
                                                                                  cgtcgGCGGCGCGTTGCTGCATAAAGGCGCGGCACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCCACGGTGACGTGTCCTTCGTTCTCGACTTCCGcc                                                                                                                                                  >  2:457794/1‑251 (MQ=255)
                                                                                           cgcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGACGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGccc                                                                                                                                                                                                                       >  2:187123/1‑172 (MQ=255)
                                                                                           cgcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCGCGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCAc                                                                                                                                                                                                                       <  1:187123/174‑1 (MQ=255)
                                                                                           cgcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGGGCCTGGGCCAATCAGAATCACCGGGGTTTCTGGAGTGTCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGGCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgaggag                                                                                                                                         <  2:479272/251‑1 (MQ=255)
                                                                                             cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCAACACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGc                                                                                                                                       <  1:420947/251‑1 (MQ=255)
                                                                                                                       cgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCg                                                                                                             <  1:56300/251‑1 (MQ=255)
                                                                                                                         cAATACCGGTGCCTGGGGCAATCATTATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCggg                                                                                                           <  1:431595/251‑1 (MQ=255)
                                                                                                                         cAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCggg                                                                                                           >  1:149967/1‑251 (MQ=255)
                                                                                                                                 gTGCCTGGGCCACTCATAATCTCCGGGGTTTCGGGATTGGCGGGCAGGCGGAAGTTATCGTTATGTTCGAGAAATACGCGGACTTCGCCCTCTTCGTCCGCGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTGCTACGTCGTAACGCACCACACCAACGGTGACGTGTGCTGCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgcagc                                                                                                     <  2:34850/249‑1 (MQ=255)
                                                                                                                                       gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGCCTTCCGCCTGCGCGGAGGCGCTGGAATACAGACGCGGCGTCAGCGGGCGCAGCCGattaat                                                                                             >  2:209706/1‑251 (MQ=255)
                                                                                                                                           cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAgc                                                                                         <  1:224171/251‑1 (MQ=255)
                                                                                                                                                  aaTCCCCGGGGTTGCTGGATGGGCTTGCAGGCGGAAGTTATCGTTAGGTGCGATAAATACGCGGACGTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGCTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCg                                                                                   <  2:291706/250‑1 (MQ=255)
                                                                                                                                                               tCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGc                                                                     >  2:406102/1‑251 (MQ=255)
                                                                                                                                                                 tGGATTGGCTGGCAGGCGGCAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCCCGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGCCTTCCGCCTGCGAGGAGGCGATGGCATACCGACGCGGCGTCCTCGGGCGCAGCAGCTTAATTAGCGCTTCGGCCTCAAGCTGTGCCg                                                                   >  2:249239/1‑251 (MQ=255)
                                                                                                                                                                 tGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCg                                                                   >  1:271559/1‑251 (MQ=255)
                                                                                                                                                                      tggctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCa                                                                                                                <  1:33157/201‑1 (MQ=255)
                                                                                                                                                                        gctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCgggggga                                                             >  2:189272/1‑247 (MQ=255)
                                                                                                                                                                              aGGCGGAATGTATCGTTATGTTCGATAAACACGCGGACTTCGCCCTCTCCGTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCAc                                                      <  2:424456/251‑1 (MQ=255)
                                                                                                                                                                                      gTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTca                                               >  1:349641/1‑250 (MQ=255)
                                                                                                                                                                                           cGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATc                                         >  1:242970/1‑251 (MQ=255)
                                                                                                                                                                                            gTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCg                                        <  1:406102/251‑1 (MQ=255)
                                                                                                                                                                                                tGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCAACACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACAATGTAAACAATCGGcct                                    >  2:374618/1‑249 (MQ=255)
                                                                                                                                                                                                     gATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGtcgtcg                               >  1:472879/1‑251 (MQ=255)
                                                                                                                                                                                                      aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCCTCAAGCTGTGcc                                                                    >  2:315570/1‑213 (MQ=255)
                                                                                                                                                                                                      aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGcc                                                                    <  1:315570/213‑1 (MQ=255)
                                                                                                                                                                                                      aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCCCCATGTCAACAATCGGCGTCGTcgc                              >  2:92992/1‑251 (MQ=255)
                                                                                                                                                                                                          aTACGAGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGACTACAGACGAGGCGTCAGCGGGCGCAGCAGATTACTTAGCGCTTCGGCATCAAGCTGTGCAGGGGCGAAACGACCCATGTCAACAAACGGCGTCGTCtcggcg                          >  2:462283/1‑251 (MQ=255)
                                                                                                                                                                                                                         cctcttcccccccgcgGTCACCTAAGAAGCTGGAGGCCCCCCCGGCCCGGGCGCGGCCGTCCACGTCGTAACGCACCACACCATCGGTGACGTTTACTTCGTTCTCGACGTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAGCAATCGGCGGCGTCGCGGCGTAATGCTGTAACTTc           <  2:242970/239‑1 (MQ=255)
                                                                                                                                                                                                                          ctctTCTTCCACGCGGTCAGCGAAGCAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCg          >  2:442513/1‑251 (MQ=255)
                                                                                                                                                                                                                              cctccctcgcgCTCATCGAAGAAGCGGGATGCACCACCGGCACTGGCGCGGCCTGCCACGTCGTAACGCACCAACTCAACGGTGACGTGTACTTCGTTCTCGTCTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGTGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGGAATGCTGTAACTTCGCtt       <  2:82741/243‑1 (MQ=255)
                                                                                                                                                                                                                                 tcCACGCGGTCAGCGAAGAAGCTGGAGGCCCCCCCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTAt    <  2:28656/250‑1 (MQ=255)
                                                                                                                                                                                                                                  ccACGCGGTCAGCGAAGAAGCTGGCGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCg  >  2:73975/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                            
TAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCG  >  NC_000913/2890385‑2890861

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