Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,890,617 G→T 6.7% R428S (CGC→AGC)  cysJ ← sulfite reductase, alpha subunit, flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,890,6170GT 6.7% 497.4 / 18.1 146R428S (CGC→AGC) cysJsulfite reductase, alpha subunit, flavoprotein
Reads supporting (aligned to +/- strand):  ref base G (65/70);  new base T (2/8);  total (67/79)
Fisher's exact test for biased strand distribution p-value = 1.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 150 total aligned reads displayed.

TCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGA  >  NC_000913/2890469‑2890765
                                                                                                                                                    |                                                                                                                                                    
tcgGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCAcgcg                                                                                                                                                      <  3:64191/149‑1 (MQ=255)
     ggcgCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAg                                                                                                                                                 <  1:724068/149‑1 (MQ=255)
        gcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGa                                                                                                                                              <  2:1563392/149‑1 (MQ=255)
        gcgTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGa                                                                                                                                              <  3:64631/149‑1 (MQ=255)
            tgctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGaa                                                                                                                                          <  4:612426/149‑1 (MQ=255)
              ctgcATAAAGGCGCGGAACGGCGCAATAGCGGTGCCTGGGCCAATCATAATCACCGTGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGTCTTCGCCCTCTTCTTCCACGCGGGCAGCGAGGAAGc                                                                                                                                        <  2:1137130/149‑1 (MQ=255)
              ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAga                                                                                                                                        <  3:492464/149‑2 (MQ=255)
              ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGc                                                                                                                                        <  3:1665301/149‑1 (MQ=255)
              ctgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTCCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGc                                                                                                                                        >  2:309177/1‑149 (MQ=255)
               tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaa                                                                                                                                       <  2:438958/149‑3 (MQ=255)
               tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaa                                                                                                                                       <  3:745501/149‑3 (MQ=255)
               tgcATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCt                                                                                                                                       >  3:801864/1‑149 (MQ=255)
                 cATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaaga                                                                                                                                     <  2:486671/149‑5 (MQ=255)
                  aTAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGa                                                                                                                                    >  3:151239/1‑149 (MQ=255)
                   tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaagagg                                                                                                                                   <  3:225077/149‑7 (MQ=255)
                   tAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                   >  1:1334683/1‑149 (MQ=255)
                      aGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCTCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGca                                                                                                                                <  2:857706/149‑1 (MQ=255)
                      aGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAg                                                                                                                                   >  2:62833/1‑146 (MQ=255)
                                  gcgcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACgg                                                                                                                    <  4:1390479/149‑1 (MQ=255)
                                    gcAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGa                                                                                                                                    >  4:783777/1‑131 (MQ=255)
                                     cAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGgcg                                                                                                                 >  4:1277518/1‑149 (MQ=255)
                                      aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGgcgc                                                                                                                >  1:256408/1‑149 (MQ=255)
                                      aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGgcgc                                                                                                                >  2:691288/1‑149 (MQ=255)
                                      aaTACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGgcgc                                                                                                                >  1:1042481/1‑149 (MQ=255)
                                             gTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGAGGCCTTc                                                                                                         <  2:926108/149‑1 (MQ=255)
                                                 cTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcg                                                                                                     >  3:1343612/1‑149 (MQ=255)
                                                  tGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgt                                                                                                    <  2:1042481/149‑1 (MQ=255)
                                                   gggCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCAcgtc                                                                                                   <  3:361081/149‑1 (MQ=255)
                                                      ccAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTa                                                                                                >  1:1557922/1‑149 (MQ=255)
                                                       cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATATATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTaa                                                                                               <  2:201614/149‑1 (MQ=255)
                                                       cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTaa                                                                                               <  4:1632520/149‑1 (MQ=255)
                                                       cAATCATAATCACCGGGGTTTCTGGATTGGCTGGCACGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCAGTCAGCGAGAAAGCTGGAGGCATCACCGGCACGGGCGCGGCCTACCACGTCGTaa                                                                                               <  4:357418/149‑1 (MQ=255)
                                                        aaTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                              >  2:690675/1‑149 (MQ=255)
                                                        aaTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                              >  3:336765/1‑149 (MQ=255)
                                                        aaTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTACCACGTCGTAAc                                                                                              >  1:1359694/1‑149 (MQ=255)
                                                        aaTCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTCCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAAc                                                                                              >  2:800731/1‑149 (MQ=255)
                                                        aaTCATAATCACCGGGGCTTCTGGATTGGCCGGCATGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTGGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGTACGGGCGCGGCCTTCCACGTCGTAAc                                                                                              >  4:1412651/1‑149 (MQ=255)
                                                          tCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaa                                                                                                                                       >  1:827732/1‑104 (MQ=255)
                                                          tCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCTCGCTGACCGCGTGGAAgaa                                                                                                                                       <  2:827732/106‑3 (MQ=255)
                                                           cATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGca                                                                                           >  2:1222974/1‑149 (MQ=255)
                                                             tAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGa                                                                                                                                    >  3:1623383/1‑106 (MQ=255)
                                                             tAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGa                                                                                                                                    <  4:1623383/106‑1 (MQ=255)
                                                              aaTCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcacca                                                                                        >  3:787028/1‑149 (MQ=255)
                                                                tCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcaccaca                                                                                      <  1:691288/149‑1 (MQ=255)
                                                                    cGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCaa                                                                                  <  2:216071/149‑1 (MQ=255)
                                                                     ggggTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAAc                                                                                 <  2:1108613/149‑1 (MQ=255)
                                                                     ggggTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAAc                                                                                 >  1:421110/1‑149 (MQ=255)
                                                                      gggTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACg                                                                                <  3:953754/149‑1 (MQ=255)
                                                                        gTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGt                                                                              <  2:421110/149‑1 (MQ=255)
                                                                         tttCTGGATTGGTTGGAAGGCGGAGGTTATGGTTTGGTTGGATAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGa                                                                            <  2:522572/149‑1 (MQ=255)
                                                                           tCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGAc                                                                           <  2:1335768/149‑1 (MQ=255)
                                                                             tGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACgt                                                                         <  1:798972/149‑1 (MQ=255)
                                                                              ggATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACgtg                                                                        >  1:235132/1‑149 (MQ=255)
                                                                                 ttGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACg                                                                                >  3:558734/1‑138 (MQ=255)
                                                                                    gctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTc                                                                  >  1:525394/1‑149 (MQ=255)
                                                                                    gctggcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTc                                                                  >  1:1058750/1‑149 (MQ=255)
                                                                                        gcAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtc                                                              >  3:378977/1‑149 (MQ=255)
                                                                                          aGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtctc                                                            >  4:601011/1‑149 (MQ=255)
                                                                                               gAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcac                                                                                          >  3:554273/1‑114 (MQ=255)
                                                                                               gAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGcac                                                                                          <  4:554273/114‑1 (MQ=255)
                                                                                                aaGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGAc                                                                           <  3:519694/128‑1 (MQ=255)
                                                                                                aaGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGAc                                                                           >  4:519694/1‑128 (MQ=255)
                                                                                                  gTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCg                                                    >  2:573100/1‑149 (MQ=255)
                                                                                                       cGTTATGTTCGATAAATACGCGGACTTCGCCCTTTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGc                                               >  3:1067902/1‑149 (MQ=255)
                                                                                                         ttATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCga                                             >  1:851605/1‑149 (MQ=255)
                                                                                                         ttATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCga                                             <  1:1611377/149‑1 (MQ=255)
                                                                                                            tGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgagga                                          <  4:985006/149‑1 (MQ=255)
                                                                                                            tGTTCGATAAATACGCGGACTTCGCCCTCTGCGTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCgagga                                          <  1:1707573/149‑1 (MQ=255)
                                                                                                                 gATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                     <  1:690675/149‑1 (MQ=255)
                                                                                                                 gATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGa                                     >  3:968783/1‑149 (MQ=255)
                                                                                                                  aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAt                                    <  4:378977/149‑1 (MQ=255)
                                                                                                                  aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAt                                    <  2:1082803/149‑1 (MQ=255)
                                                                                                                  aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAt                                    >  1:1029757/1‑149 (MQ=255)
                                                                                                                  aTAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAt                                    <  2:1029757/149‑1 (MQ=255)
                                                                                                                      aTACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGaa                                >  2:1107586/1‑149 (MQ=255)
                                                                                                                        aCGCGGACTTCGCCTTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATa                              <  4:1247045/149‑1 (MQ=255)
                                                                                                                        aCGCGGACATCGCCGTCTTCTTCCACGCGGTGAGCGAGGAAGCTTGAGGCACCACCGGCACGGGCGCGGCCTTCGACGTCGTAACGCACCACCCCAACGGTGACGTGTACTTCGTTCTCGTCTTCCGCCTGCGAGGAGGCGATGGAATa                              <  3:1412651/149‑1 (MQ=255)
                                                                                                                          gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACg                                                                                <  3:1686745/97‑1 (MQ=255)
                                                                                                                          gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACg                                                                                >  4:1686745/1‑97 (MQ=255)
                                                                                                                          gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACgtgt                                                                       >  1:1228829/1‑106 (MQ=255)
                                                                                                                          gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                            >  2:981538/1‑149 (MQ=255)
                                                                                                                          gcgGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                            <  1:179594/149‑1 (MQ=255)
                                                                                                                            ggACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGa                          >  2:1593457/1‑149 (MQ=255)
                                                                                                                              aCTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGCGGCAACACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACTACACCAACGGTGACGTGTACTTCGTTCTCGACTGCCGCCTGCGAGGAGGCGATGGAATACAGAcg                        <  2:478479/149‑1 (MQ=255)
                                                                                                                               cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc                       <  1:373283/149‑1 (MQ=255)
                                                                                                                               cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc                       <  1:1222974/149‑1 (MQ=255)
                                                                                                                               cTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgc                       <  2:1024404/149‑1 (MQ=255)
                                                                                                                                ttCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAt                                    >  4:1325234/1‑135 (MQ=255)
                                                                                                                                ttCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGAt                                    <  3:1325234/135‑1 (MQ=255)
                                                                                                                                ttCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGAcgcg                      >  3:1058869/1‑149 (MQ=255)
                                                                                                                                  cGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggc                    <  1:1107586/149‑1 (MQ=255)
                                                                                                                                   gCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGAACAACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACgcggcg                   >  3:725221/1‑149 (MQ=255)
                                                                                                                                     ccTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTtctc                                                            >  2:971595/1‑106 (MQ=255)
                                                                                                                                     ccTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGt                  >  2:286504/1‑148 (MQ=255)
                                                                                                                                     ccTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTc                 >  4:754189/1‑149 (MQ=255)
                                                                                                                                       tcttctTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAg               >  4:851782/1‑149 (MQ=255)
                                                                                                                                        cttcttCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACa                            >  2:1512195/1‑135 (MQ=255)
                                                                                                                                                aCGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCaa                                                                                  >  1:621694/1‑73 (MQ=255)
                                                                                                                                                aCGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCaa                                                                                  <  2:621694/73‑1 (MQ=255)
                                                                                                                                                aCGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTCCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCgcag      >  3:156230/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCGGCGGCGCGTTGCTGCATAAAGGCGCGGAACGGCGCAATACCGGTGCCTGGGCCAATCATAATCACCGGGGTTTCTGGATTGGCTGGCAGGCGGAAGTTATCGTTATGTTCGATAAATACGCGGACTTCGCCCTCTTCTTCCACGCGGTCAGCGAGGAAGCTGGAGGCACCACCGGCACGGGCGCGGCCTTCCACGTCGTAACGCACCACACCAACGGTGACGTGTACTTCGTTCTCGACTTCCGCCTGCGAGGAGGCGATGGAATACAGACGCGGCGTCAGCGGGCGCAGCAGA  >  NC_000913/2890469‑2890765

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