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.

ACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGC  >  NC_000913/2888673‑2888967
                                                                                                                                                    |                                                                                                                                                   
  acGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcg                                                                                                                                                  >  1:1880856/1‑149 (MQ=255)
    gCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcggc                                                                                                                                                >  2:946615/1‑149 (MQ=255)
      aTCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGt                                                                                                                                              >  1:1367736/1‑149 (MQ=255)
       tCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTg                                                                                                                                             >  2:643716/1‑149 (MQ=255)
       tCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTg                                                                                                                                             >  1:1692483/1‑149 (MQ=255)
        cACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGa                                                                                                                                            <  4:833989/149‑1 (MQ=255)
         aCGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTg                                                                                                                                                    >  4:1567567/1‑140 (MQ=255)
            atatGCTCATCGTTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑TCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGc                                                                                                                                        >  4:1189200/1‑149 (MQ=255)
              atGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcat                                                                                                                                      <  1:1437091/149‑1 (MQ=255)
                  tCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcattcat                                                                                                                                  >  2:615171/1‑149 (MQ=255)
                    aTCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTa                                                                                                                                >  1:1148372/1‑149 (MQ=255)
                     tCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTaa                                                                                                                               >  3:1320335/1‑149 (MQ=255)
                           acacCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGct                                                                                                                         <  2:1817263/149‑1 (MQ=255)
                               cATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTaa                                                                                                                               >  2:647302/1‑139 (MQ=255)
                                aTGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCtt                                                                                                                    <  3:911284/149‑1 (MQ=255)
                                 tGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCttt                                                                                                                   <  3:1667797/149‑1 (MQ=255)
                                   tttCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTgg                                                                                                                 <  1:643716/149‑1 (MQ=255)
                                     tCGCCATTAAATTATCGATGTTGTCGATAAAA‑ACGGCCGGGAACCCGCCGCTTTCGCCCTTGGCCAGCGGCAAAGTCGGGATTGCAACGCAGGCCACCAAGTTTCCCCCCTCGGGCGGCACGGCATCGATAAACCCGCCCCCTTTgggg                                                                                                               >  1:1770516/1‑147 (MQ=255)
                                          aTTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtctt                                                                                                          >  2:1883607/1‑149 (MQ=255)
                                              aaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcggc                                                                                                                                                >  2:1370383/1‑107 (MQ=255)
                                              aaTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTgcggc                                                                                                                                                <  1:1370383/107‑1 (MQ=255)
                                               aTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtcttc                                                                                                         >  4:14940/1‑145 (MQ=255)
                                               aTTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtcttc                                                                                                         <  3:14940/145‑1 (MQ=255)
                                                ttATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGAt                                                                                                    >  3:1396426/1‑149 (MQ=255)
                                                    cGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtc                                                                                                <  2:1692483/149‑1 (MQ=255)
                                                    cGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtc                                                                                                <  3:1106680/149‑1 (MQ=255)
                                                    cGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtc                                                                                                <  1:1034457/149‑1 (MQ=255)
                                                          tgtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttc                                                                                          <  1:615171/149‑1 (MQ=255)
                                                           gtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         >  2:581575/1‑149 (MQ=255)
                                                           gtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         <  2:347239/149‑1 (MQ=255)
                                                           gtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         <  1:20422/149‑1 (MQ=255)
                                                           gtCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGCGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         <  3:1189200/149‑1 (MQ=255)
                                                            tCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGAtc                                                                                                            <  3:147667/129‑1 (MQ=255)
                                                            tCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctc                                                                                        >  1:517335/1‑149 (MQ=255)
                                                            tCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctc                                                                                        <  2:850937/149‑1 (MQ=255)
                                                            tCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctc                                                                                        <  4:1396426/149‑1 (MQ=255)
                                                            gtgaataaaGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCACGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctc                                                                                        <  2:1770516/143‑1 (MQ=255)
                                                              gATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctcgc                                                                                      >  3:965406/1‑149 (MQ=255)
                                                              gATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctcgc                                                                                      >  4:1160665/1‑149 (MQ=255)
                                                              gATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCCGCAGCGTGaaaa                                                                                                                                         >  1:1697866/1‑95 (MQ=255)
                                                                tAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCACTCTTTGGCGATCTTCTCGATCTtcgctttctcgctt                                                                                    <  1:824835/149‑1 (MQ=255)
                                                                  aaaGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGATTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttctcgctttc                                                                                  >  3:1599892/1‑149 (MQ=255)
                                                                    aGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCCAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCg                                                                                <  4:965406/149‑1 (MQ=255)
                                                                    aGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCg                                                                                <  3:1160665/149‑1 (MQ=255)
                                                                        ggCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         >  1:438419/1‑136 (MQ=255)
                                                                        ggCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTtcgctttct                                                                                         <  2:438419/136‑1 (MQ=255)
                                                                        ggCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTac                                                                            >  2:1526742/1‑149 (MQ=255)
                                                                             gAAACGCTCTGCTTCTG‑CCATCGCCAGCGGGCAAGTCGGGAACGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCTCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCgg                                                                       <  2:548153/149‑1 (MQ=255)
                                                                             gAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTATTTGGCGATCTTCTCGATCTTCGCTTTATCGCTTTCCGGTACACCgg                                                                       >  2:1351621/1‑149 (MQ=255)
                                                                                aCGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCga                                                                    <  2:1843041/149‑1 (MQ=255)
                                                                                  gCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcatt                                                                                                                                     >  1:21639/1‑82 (MQ=255)
                                                                                  gCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGcatt                                                                                                                                     <  2:21639/82‑1 (MQ=255)
                                                                                   ctctGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatga                                                                 >  4:198344/1‑149 (MQ=255)
                                                                                    tctGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCgatgat                                                                <  4:376801/149‑1 (MQ=255)
                                                                                     ctGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATc                                                               <  4:973574/149‑1 (MQ=255)
                                                                                     ctGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATc                                                               <  2:315419/149‑1 (MQ=255)
                                                                                        cTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGa                                                            >  2:1108869/1‑149 (MQ=255)
                                                                                         ttCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGAt                                                           >  2:612445/1‑149 (MQ=255)
                                                                                           ccG‑CCATCGCCAGCGGTCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTc                                                         >  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)
                                                                                                                                                    |                                                                                                                                                   
ACACGCATCACGATATGCTCATCGCTGACACCATGTTTCGCCATTAAATTATCGATGTTGTCGATAAAAGACGGCAGGAAACGCTCTGCTTCCG‑CCATCGCCAGCGGGCAAGTCGGGAATGACACGCAAGCCATCGAGTTTTCACGCTGCGGCGTGACGGCATTCATTAACCCGCTCTCTTTGGCGATCTTCTCGATCTTCGCTTTCTCGCTTTCCGGTACACCGGCGATGATCAGATTCTGGTTCGCCGTAATGCGGAAATCGCCTTTGTGGATCTTCGCGATCTCCAGCAGGC  >  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: ‑