Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,888,820 T→C 6.1% Q427R (CAG→CGG)  cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,888,8200TC 6.1% 501.6 / 15.6 148Q427R (CAG→CGG) cysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base T (65/74);  new base C (5/4);  total (70/78)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

Only 100 of 148 total aligned reads displayed.

ACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGC  >  NC_000913/2888673‑2888967
                                                                                                                                                    |                                                                                                                                                   
  acGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcg                                                                                                                                                  >  1:1880856/1‑149 (MQ=255)
    gCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcggc                                                                                                                                                >  2:946615/1‑149 (MQ=255)
      aTCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt                                                                                                                                              >  1:1367736/1‑149 (MQ=255)
       tCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTg                                                                                                                                             >  2:643716/1‑149 (MQ=255)
       tCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTg                                                                                                                                             >  1:1692483/1‑149 (MQ=255)
        cACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGa                                                                                                                                            <  4:833989/149‑1 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTg                                                                                                                                                    >  4:1567567/1‑140 (MQ=255)
            atatGCTCATCGTTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGTCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGc                                                                                                                                        >  4:1189200/1‑149 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcat                                                                                                                                      <  1:1437091/149‑1 (MQ=255)
                  tCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattcat                                                                                                                                  >  2:615171/1‑149 (MQ=255)
                    aTCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa                                                                                                                                >  1:1148372/1‑149 (MQ=255)
                     tCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTaa                                                                                                                               >  3:1320335/1‑149 (MQ=255)
                           acacCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGct                                                                                                                         <  2:1817263/149‑1 (MQ=255)
                               cATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTaa                                                                                                                               >  2:647302/1‑139 (MQ=255)
                                aTGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCtt                                                                                                                    <  3:911284/149‑1 (MQ=255)
                                 tGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCttt                                                                                                                   <  3:1667797/149‑1 (MQ=255)
                                   tttCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTgg                                                                                                                 <  1:643716/149‑1 (MQ=255)
                                     tCGCCATTAAATTATCGATGTTGTCGATAAAAACGGCCGGGAACCCGCCGCTTTCGCCCTTGGCCAGCGGCAAAGTCGGGATTGCAACGCAGGCCACCAAGTTTCCCCCCTCGGGCGGCACGGCATCGATAAACCCGCCCCCTTTgggg                                                                                                               >  1:1770516/1‑147 (MQ=255)
                                          aTTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtctt                                                                                                          >  2:1883607/1‑149 (MQ=255)
                                              aaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcggc                                                                                                                                                >  2:1370383/1‑107 (MQ=255)
                                              aaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcggc                                                                                                                                                <  1:1370383/107‑1 (MQ=255)
                                               aTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtcttc                                                                                                         >  4:14940/1‑145 (MQ=255)
                                               aTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtcttc                                                                                                         <  3:14940/145‑1 (MQ=255)
                                                ttATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGAt                                                                                                    >  3:1396426/1‑149 (MQ=255)
                                                    cGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtc                                                                                                <  2:1692483/149‑1 (MQ=255)
                                                    cGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtc                                                                                                <  3:1106680/149‑1 (MQ=255)
                                                    cGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtc                                                                                                <  1:1034457/149‑1 (MQ=255)
                                                          tgtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttc                                                                                          <  1:615171/149‑1 (MQ=255)
                                                           gtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         >  2:581575/1‑149 (MQ=255)
                                                           gtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         <  2:347239/149‑1 (MQ=255)
                                                           gtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         <  1:20422/149‑1 (MQ=255)
                                                           gtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGCGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         <  3:1189200/149‑1 (MQ=255)
                                                            tCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtc                                                                                                            <  3:147667/129‑1 (MQ=255)
                                                            tCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctc                                                                                        >  1:517335/1‑149 (MQ=255)
                                                            tCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctc                                                                                        <  2:850937/149‑1 (MQ=255)
                                                            tCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctc                                                                                        <  4:1396426/149‑1 (MQ=255)
                                                            gtgaataaaGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCACGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctc                                                                                        <  2:1770516/143‑1 (MQ=255)
                                                              gATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctcgc                                                                                      >  3:965406/1‑149 (MQ=255)
                                                              gATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctcgc                                                                                      >  4:1160665/1‑149 (MQ=255)
                                                              gATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCCGCAGCGTGaaaa                                                                                                                                         >  1:1697866/1‑95 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCACTCTTTGGCGATCTTCTCGATCTtcgctttctcgctt                                                                                    <  1:824835/149‑1 (MQ=255)
                                                                  aaaGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGATTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctcgctttc                                                                                  >  3:1599892/1‑149 (MQ=255)
                                                                    aGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCCAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCg                                                                                <  4:965406/149‑1 (MQ=255)
                                                                    aGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCg                                                                                <  3:1160665/149‑1 (MQ=255)
                                                                        ggCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         >  1:438419/1‑136 (MQ=255)
                                                                        ggCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         <  2:438419/136‑1 (MQ=255)
                                                                        ggCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTac                                                                            >  2:1526742/1‑149 (MQ=255)
                                                                             gAAACGCTCTGCTTCTGCCATCGCCAGCGGGCAAGTCGGGAACGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCTCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCgg                                                                       <  2:548153/149‑1 (MQ=255)
                                                                             gAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTATTTGGCGATCTTCTCGATCTTCGCTTTATCGCTTTCCGGTACACCgg                                                                       >  2:1351621/1‑149 (MQ=255)
                                                                                aCGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCga                                                                    <  2:1843041/149‑1 (MQ=255)
                                                                                  gCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcatt                                                                                                                                     >  1:21639/1‑82 (MQ=255)
                                                                                  gCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcatt                                                                                                                                     <  2:21639/82‑1 (MQ=255)
                                                                                   ctctGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatga                                                                 >  4:198344/1‑149 (MQ=255)
                                                                                    tctGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                                                <  4:376801/149‑1 (MQ=255)
                                                                                     ctGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATc                                                               <  4:973574/149‑1 (MQ=255)
                                                                                     ctGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATc                                                               <  2:315419/149‑1 (MQ=255)
                                                                                        cTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                                            >  2:1108869/1‑149 (MQ=255)
                                                                                         ttCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGAt                                                           >  2:612445/1‑149 (MQ=255)
                                                                                           ccGCCATCGCCAGCGGTCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTc                                                         >  2:1648128/1‑149 (MQ=255)
                                                                                                    gCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGc                                                 <  1:1883607/149‑1 (MQ=255)
                                                                                                     ccAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCg                                                                                <  2:509068/117‑1 (MQ=255)
                                                                                                     ccAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCg                                                                                >  1:509068/1‑117 (MQ=255)
                                                                                                     ccAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGcc                                                <  2:375501/149‑1 (MQ=255)
                                                                                                          gggCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTg                                                       >  1:1635211/1‑137 (MQ=255)
                                                                                                            gCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGc                                         <  1:1108869/149‑1 (MQ=255)
                                                                                                             cAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCg                                        <  2:1831104/149‑1 (MQ=255)
                                                                                                              aaGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctcgc                                                                                      >  4:174415/1‑102 (MQ=255)
                                                                                                              aaGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCt                                                        <  3:1840948/132‑1 (MQ=255)
                                                                                                                 tCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGaaa                                    >  3:536266/1‑149 (MQ=255)
                                                                                                                   gggAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATc                                  <  2:1391773/149‑1 (MQ=255)
                                                                                                                     gAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTc                                                         >  2:860097/1‑124 (MQ=255)
                                                                                                                     gAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTc                                                         <  1:860097/124‑1 (MQ=255)
                                                                                                                     gAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGc                                >  2:922568/1‑149 (MQ=255)
                                                                                                                       aTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACAcc                                                                         >  2:919668/1‑106 (MQ=255)
                                                                                                                       aTGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACAcc                                                                         <  1:919668/106‑1 (MQ=255)
                                                                                                                          acacGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtg                           <  3:198344/149‑1 (MQ=255)
                                                                                                                            acGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgtg                         >  1:338806/1‑149 (MQ=255)
                                                                                                                              gCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                                            <  3:1759958/112‑1 (MQ=255)
                                                                                                                              gCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                                            >  4:551141/1‑112 (MQ=255)
                                                                                                                              gCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                                            >  4:1759958/1‑112 (MQ=255)
                                                                                                                                aaGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgt                          >  2:464198/1‑144 (MQ=255)
                                                                                                                                aaGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgt                          <  1:1560061/144‑1 (MQ=255)
                                                                                                                                aaGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgt                          <  1:464198/144‑1 (MQ=255)
                                                                                                                                aaGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgt                          >  2:1560061/1‑144 (MQ=255)
                                                                                                                                     aTCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCgg                                                                       <  1:1825455/94‑1 (MQ=255)
                                                                                                                                     aTCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTcgc                >  2:521947/1‑149 (MQ=255)
                                                                                                                                     aTCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTcgc                >  3:855448/1‑149 (MQ=255)
                                                                                                                                      tCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTcgcg               >  1:373092/1‑149 (MQ=255)
                                                                                                                                        gAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGAt             <  1:612445/149‑1 (MQ=255)
                                                                                                                                         aGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgt                          >  4:476858/1‑135 (MQ=255)
                                                                                                                                         aGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTtgt                          <  3:476858/135‑1 (MQ=255)
                                                                                                                                           ttttCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGAtctc          >  4:1849908/1‑149 (MQ=255)
                                                                                                                                             ttCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGTCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCca        <  4:1320335/149‑1 (MQ=255)
                                                                                                                                              tCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCcag       <  1:1526742/149‑1 (MQ=255)
                                                                                                                                              tCACGCCGCAGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGa                                      <  2:1012206/118‑1 (MQ=255)
                                                                                                                                              tCACGCCGCAGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGa                                      >  1:1012206/1‑118 (MQ=255)
                                                                                                                                              tCACGCCGCAGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCtt                   <  2:119910/137‑1 (MQ=255)
                                                                                                                                              tCACGCCGCAGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCtt                   >  1:119910/1‑137 (MQ=255)
                                                                                                                                                aCGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCcagca     >  3:985535/1‑149 (MQ=255)
                                                                                                                                                 cGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTgg                        >  1:935476/1‑129 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCGCCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGC  >  NC_000913/2888673‑2888967

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: