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,6260GA95.5% 210.0 / NA 68L425F (CTC→TTC) cysJsulfite reductase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base A (36/28);  minor base C (3/0);  total (40/28)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACC  >  NC_000913/2890384‑2890867
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
gTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTgg                                                                                                                                                                                                                                           <  2:163294/251‑1 (MQ=255)
  aGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGGTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTgg                                                                                                                                                                                                                                           >  1:82335/1‑249 (MQ=255)
     gcTGCCACTCCCCCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGGTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGca                                                                                                                                                                                                                                      <  2:59626/251‑1 (MQ=255)
        gCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGcac                                                                                                                                                                                                                                     <  1:42502/249‑1 (MQ=255)
            cTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCgg                                                                                                                                                                                                                               >  1:28721/1‑251 (MQ=255)
                                       aaaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                                                                                                                                    <  1:232383/251‑1 (MQ=255)
                                       aaaGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                                                                                                                                    >  1:94062/1‑251 (MQ=255)
                                           tGCGGATTACCAACGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcac                                                                                                                                                                                                >  2:197426/1‑251 (MQ=255)
                                            gCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcacc                                                                                                                                                                                               <  1:316634/251‑1 (MQ=255)
                                              ggATTACCAAAGAACAGCCAGTGTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGAAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTCTCGTTATGTTCGATAAATACGCGGACTTATCCCTCTTCTTCCACGCGGTCAGCGAGGAAGATGGAGGCACCACCGGCACGGGCGCGGCCTTCCTCGTCGTAACGcacca                                                                                                                                                                                              >  2:83021/1‑250 (MQ=255)
                                                aTTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATACAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGCAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccacac                                                                                                                                                                                           >  2:56541/1‑251 (MQ=255)
                                                        aGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCCCCACACCAACggcg                                                                                                                                                                                   >  2:92583/1‑249 (MQ=255)
                                                             aGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACgtg                                                                                                                                                                              >  1:397646/1‑251 (MQ=255)
                                                                                       ggcggcGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCg                                                                                                                                                    >  1:136927/1‑251 (MQ=255)
                                                                                        gcggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGTTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACATGTACTTCGTTCTCGACTTCCGCCTGCga                                                                                                                                                   <  1:162116/251‑1 (MQ=255)
                                                                                              cgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGc                                                                                                                                             >  1:259950/1‑251 (MQ=255)
                                                                                                 tgctgcATAAAGGAGCGGAACGGCGCAATACCGGTGCCGGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGTGGAAGTTATCGTTATGTTCGATAAATACGCGGACTGCGCCCTCTTCCTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCTCCGGCACGGGCGCGTCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGTCTTCCGCCTGCGAGGAGGCGa                                                                                                                                           <  2:73762/250‑1 (MQ=255)
                                                                                                   ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATgg                                                                                                                                        >  2:66389/1‑251 (MQ=255)
                                                                                                              cgcgGAACGACGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc                                                                                                                             >  1:86288/1‑251 (MQ=255)
                                                                                                               gcgGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgcg                                                                                                                            <  1:234715/251‑1 (MQ=255)
                                                                                                                 ggAACGGCGCAATCCCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGAGGACGATATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGACGAACCTGGAGGCACCACCGGCACGGGCGCGGCCTTCACCGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGCATTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggg                                                                                                                          >  2:402376/1‑250 (MQ=255)
                                                                                                                  gAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggcg                                                                                                                         >  2:79345/1‑251 (MQ=255)
                                                                                                                     cGGCGCAATACCGGTGCCTGGACCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTCCTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTc                                                                                                                       >  2:319052/1‑250 (MQ=255)
                                                                                                                                  gTGCCTGGGCCAATCAGAATCGCCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAGGGAATACAGACGCGGCGTCAGCGGGCGcagcag                                                                                                         <  2:223052/251‑1 (MQ=255)
                                                                                                                                    gCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCCGCGGGCgcagc                                                                                                           >  1:398971/1‑247 (MQ=255)
                                                                                                                                        gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCTCCACACCAACGGTTACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGattaat                                                                                                   <  1:97593/251‑1 (MQ=255)
                                                                                                                                        gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGatta                                                                                                     <  1:245701/249‑1 (MQ=255)
                                                                                                                                         ggCCAATCATAATCACCGGGGTTTCTGGAGTGGCTGGCAGGCGGAAGTTATCTGTATGTTCGATAAATACCCGGCCTTCGCCCTCTTCGTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTGCCGCCTGCGAAGAGGCGATGGCATACAGACGCGGCGTCAGCGGGCGCAGCAGattaatt                                                                                                  <  2:310416/251‑1 (MQ=255)
                                                                                                                                           ccAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGattaatt                                                                                                  >  2:113145/1‑249 (MQ=255)
                                                                                                                                             aagcgtAATCACCGGGGTTTCGGGATTGGCTGGCAGGCGGAAGTTATCGGTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTGCCACGCGGTCAGCGAAGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTGCTCTACTTCCGCCTGCGAGGAGGGGATGGAATACAGACGCGGCGTCAGCGGGCGCAGAACAGTAATTAgcg                                                                                              <  2:57053/246‑1 (MQ=255)
                                                                                                                                             aaTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAgcg                                                                                              <  1:63908/251‑1 (MQ=255)
                                                                                                                                                        ccGGGGTTTCTGGAGTGGCTGGAATGGGGAAGTTATCGGTATGTTCGATAAATACGCGGACTTCGCCCTCGTCTCCCACGCGGTCAGCGAAGAAGCGGGAGGCACCACAGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGACTGCGAGGAGGCGATGGAAGACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCa                                                                                   <  2:417926/251‑1 (MQ=255)
                                                                                                                                                             gTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGTAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCtg                                                                              <  2:397646/251‑1 (MQ=255)
                                                                                                                                                                          ctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGaa                                                                  >  1:23935/1‑250 (MQ=255)
                                                                                                                                                                           tggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAAc                                                                >  1:413949/1‑251 (MQ=255)
                                                                                                                                                                                 gCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCAcc                                                           >  2:104693/1‑250 (MQ=255)
                                                                                                                                                                                   ggAAGTTATCGTTATGTTCGATAAATACGCGGACTGCGCCCTCTTCTTCCACGCGGTAAGCGAAGAAGCTGGCGGCACCACCGGCACGGGCGCGGCCTTCCACGACGTAACGCACCACACCAAAGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGACGATGGCATACAGACGCGGCGTCATCGGGCGCAGCAGATTAATTAGCGCTTCGGCCTCATGCTCTGCCGGGGAGCAACGCCCCATg                                                        >  2:322632/1‑251 (MQ=255)
                                                                                                                                                                                    gAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGt                                                       >  2:250615/1‑251 (MQ=255)
                                                                                                                                                                                     aaGTTATCGTTATGTTCGATAAATACGCGGACTTCGCGCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCGCCGGCACGGGCGCGGCCTTCCACGTCGTAACGCCCCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTTCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTc                                                      <  2:357996/251‑1 (MQ=255)
                                                                                                                                                                                       gTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGACGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGCGAAACGCCCCATGTcac                                                    >  2:338824/1‑250 (MQ=255)
                                                                                                                                                                                         tATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGCTGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTcaac                                                   >  2:3988/1‑250 (MQ=255)
                                                                                                                                                                                          atttgtATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGCCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACTGTGACGTGTACTTCGTTCTCGACTTCCGCCGGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATGAATTATCGCTTCGGCATCAATCTGTGCCGGTGAGAAACGATCCATGTcaaca                                                  <  2:198074/245‑1 (MQ=255)
                                                                                                                                                                                             gTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCa                                                                                                                                                                                                             >  1:332876/1‑92 (MQ=255)
                                                                                                                                                                                             gTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCa                                                                                                                                                                                                             <  2:332876/92‑1 (MQ=255)
                                                                                                                                                                                               tATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGc                                            >  1:167868/1‑251 (MQ=255)
                                                                                                                                                                                                aTGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGcg                                           <  1:44887/251‑1 (MQ=255)
                                                                                                                                                                                                aTGTTCGAGAAATACGCGGACGTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCAGGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCGGTTCCGGGGATAAACGCACCATGTCAACAATCGGc                                            <  2:314045/250‑1 (MQ=255)
                                                                                                                                                                                                     cggtAAATACGCGGGCTTCGCCCTCTTCTTCCACGCGGTCAGCGACGAAGCTGGAGGCACCACCGGCACGGACGCGGCCTTCCACGTCGTAACGCACCACACCAACGGCGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAAACGCGGCGTCAACGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCGCCATGTCAACAATCGGCgtcgtc                                      >  2:357622/4‑251 (MQ=255)
                                                                                                                                                                                                     cGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCgtcgtc                                      >  1:385856/1‑251 (MQ=255)
                                                                                                                                                                                                     cGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGAGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGCGAAACGCACCATATCACCACTCGGCgtcgtc                                      >  1:373235/1‑251 (MQ=255)
                                                                                                                                                                                                      gATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGtcgtcg                                     <  1:109082/251‑1 (MQ=255)
                                                                                                                                                                                                      gATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCAACACCGGCACGGGCGCGGAATTCCACGTCGTAACGCAACACAACAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGAAGATTAATTAGCGCATCGGCATCAAGCTGTGCCGGGGAGAAACGCACAATGTCAACAATCGGcgtcgt                                       >  1:249804/1‑249 (MQ=255)
                                                                                                                                                                                                       aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTcgc                                    <  1:113145/251‑1 (MQ=255)
                                                                                                                                                                                                       aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTAAGGCATCAAGCTGTGCCGGGGAGAAACGCAACAGGTCAACAATCGGCGTCGTcgc                                    >  1:17378/1‑251 (MQ=255)
                                                                                                                                                                                                          aaTACGCGGACTTCGCCCTCTTCTTCCACGCTGTCAGCGAAGACGCTGGAGGCACCACCGACACGTGAGCGGCCTGCCACGTCGTAACGCACCACACCAACGGTGCCGTGTACTTCGTCCTCGACTTCCTCCTCCGAGGCGGCGCTGGAATACAGACGCAGCCTCATCGCGCGCAGCGCATTAATTAGCGCTTCGCCATCAAGCTGTGCCGGGGAGAAACGCACCATCTCAACAATCGGCGTCGTCgcggc                                 >  1:385565/1‑251 (MQ=255)
                                                                                                                                                                                                           aTACGCGGACTTCGACCTATTCTTCCCCGCGGTCAGCGAAGAAGCTGGAGGCCCCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGcc                                                                                                                                                        >  1:104838/1‑131 (MQ=255)
                                                                                                                                                                                                             aCGCGGACTTCACCATCTTCTTCCACCCGGTCAGCGCAGACGCTCGAGGCAACACCGGCACGGGCGCGGCCCTCCACGCCGTAACGCACCACACCCACGGTGCCGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACCGACGCGGCGTCAGCGGGCGCAGCAGATGCATTAGCGCTTCGGCATCAAGCTATGCCGGGCAGAAACGCACCATGTCAACAATCGGCGACGTCGCGGCGt                               >  2:270908/1‑250 (MQ=255)
                                                                                                                                                                                                              cgcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGCCTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTa                              >  1:437747/1‑250 (MQ=255)
                                                                                                                                                                                                                cgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGCGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATg                           >  2:261864/1‑251 (MQ=255)
                                                                                                                                                                                                                  gtcTTCTCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCGTCCCCGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCCCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGTCGTAATGCt                         <  2:414666/249‑1 (MQ=255)
                                                                                                                                                                                                                       cGCCCTCTTCCTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTaa                     <  1:418575/250‑1 (MQ=255)
                                                                                                                                                                                                                        gCCCTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACt                   <  1:10626/251‑1 (MQ=255)
                                                                                                                                                                                                                         cccTCTTCTTCCACGCGGTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACtt                  <  1:75883/251‑1 (MQ=255)
                                                                                                                                                                                                                          ccTCTTCTTCCACGCGGTCAGCGAAGAAGATGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAAAAATCGGCGTCGTCGCGGCGTAATGCTGTAAcgtc                 >  2:438604/1‑248 (MQ=255)
                                                                                                                                                                                                                                  tcctcgcgGTCAGCGAAGAAGCTGGAGGCACCACCGGCGCGGGCGCGGCCTTCCACGGCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATATAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTAt          <  2:92209/246‑1 (MQ=255)
                                                                                                                                                                                                                                         ggTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCccacc  >  2:93416/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gTCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCccacc  >  2:337534/1‑250 (MQ=255)
                                                                                                                                                                                                                                           tCAGCGAAGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCccacc  >  2:426950/1‑249 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                 
GTAGCGCTGCCACTCCACCTGGTACAGGAAGTCTTCCGTAAAGTGCGGATTACCAAAGAACAGCCAGTTTTTACCTGGCGCTTCGTCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGATTAATTAGCGCTTCGGCATCAAGCTGTGCCGGGGAGAAACGCACCATGTCAACAATCGGCGTCGTCGCGGCGTAATGCTGTAACTTCGCTTTATCGCCCACC  >  NC_000913/2890384‑2890867

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