Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,585,149 Δ1 bp 7.0% coding (1338/2280 nt) ydeP ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,585,1490G. 7.0% 756.9 / 49.6 214coding (1338/2280 nt)ydePputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (102/97);  new base . (9/6);  total (111/103)
Fisher's exact test for biased strand distribution p-value = 5.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-02

Only 100 of 216 total aligned reads displayed.

AGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCC  >  NC_000913/1585001‑1585296
                                                                                                                                                    |                                                                                                                                                    
    gTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCaa                                                                                                                                                  >  4:292508/1‑149 (MQ=255)
     ttAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAAt                                                                                                                                                 <  3:777916/149‑1 (MQ=255)
     ttAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAAt                                                                                                                                                 >  2:702657/1‑149 (MQ=255)
     ttAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAAt                                                                                                                                                 <  4:1304551/149‑1 (MQ=255)
         gCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgca                                                                                                                                            <  3:148967/149‑1 (MQ=255)
         gCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgca                                                                                                                                            >  3:216075/1‑149 (MQ=255)
         gCTTAGTGGCTACGTGTACCGCCAAACCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgca                                                                                                                                            <  1:379520/149‑1 (MQ=255)
           ttAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgcagc                                                                                                                                          >  3:1639350/1‑149 (MQ=255)
             aGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCAt                                                                                                                                        <  1:1724934/149‑1 (MQ=255)
                  cTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCa                                                                                                                                   >  2:42699/1‑149 (MQ=255)
                   tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGACCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAg                                                                                                                                  <  3:969963/149‑1 (MQ=255)
                       tgtACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGc                                                                                                                              <  1:702657/149‑1 (MQ=255)
                         tACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCAt                                                                                                                            <  3:828069/149‑1 (MQ=255)
                              ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGt                                                                                                                       <  4:1647597/149‑1 (MQ=255)
                                aaaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTgg                                                                                                                     >  2:19515/1‑149 (MQ=255)
                                 aaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTggg                                                                                                                    <  2:1312045/149‑1 (MQ=255)
                                 aaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTggg                                                                                                                    <  1:1044164/149‑1 (MQ=255)
                                   tCTAATTGCGTTAACGGTACAGCGCTCGCTTCGCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCAAGAGGTGGGGt                                                                                                                  >  4:1208364/1‑149 (MQ=255)
                                     tAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGa                                                                                                                >  2:613525/1‑149 (MQ=255)
                                      aaTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGaa                                                                                                               <  1:1210451/149‑1 (MQ=255)
                                       aTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTcc                                                                                                                                    >  3:295192/1‑127 (MQ=255)
                                       aTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                                                              >  1:1417441/1‑149 (MQ=255)
                                        ttGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                                                                             <  2:1264472/149‑1 (MQ=255)
                                         tGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGt                                                                                                                       <  3:1307690/138‑1 (MQ=255)
                                         tGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGcc                                                                                                            >  1:117008/1‑149 (MQ=255)
                                            gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                         >  1:730738/1‑149 (MQ=255)
                                            gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                         >  1:1136163/1‑149 (MQ=255)
                                             ttAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTggg                                                                                                                    <  2:1398200/137‑1 (MQ=255)
                                              tAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgc                                                                                                       >  4:891190/1‑149 (MQ=255)
                                               aaCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcg                                                                                                      <  4:1639350/149‑1 (MQ=255)
                                                        gcgcTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgcagc                                                                                                                                          <  1:1278422/104‑1 (MQ=255)
                                                         cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGaa                                                                                                               <  3:40074/130‑1 (MQ=255)
                                                         cgctcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGa                                                                                            <  1:195927/149‑1 (MQ=255)
                                                          gctcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGt                                                                                                                       <  3:1601001/121‑1 (MQ=255)
                                                            tcgctTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGa                                                                                         >  1:404732/1‑149 (MQ=255)
                                                              gctTCCCGATCTGGCATTGCCAGCGCAAAATTGACCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTgca                                                                                                                                            >  3:1609983/1‑96 (MQ=255)
                                                                ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                     >  1:870657/1‑149 (MQ=255)
                                                                  cccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGa                                                                                   >  1:1754731/1‑149 (MQ=255)
                                                                    cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                 >  2:1106476/1‑149 (MQ=255)
                                                                    cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                 <  1:941652/149‑1 (MQ=255)
                                                                    cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                 <  1:1106476/149‑1 (MQ=255)
                                                                     gATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAc                                                                                <  2:1417441/149‑1 (MQ=255)
                                                                           gCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCa                                                                          >  1:1777091/1‑149 (MQ=255)
                                                                                   aGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGttt                                                                  <  4:825401/149‑1 (MQ=255)
                                                                                     cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtct                                                                <  3:891190/149‑1 (MQ=255)
                                                                                      gcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGt                                                                                                                       >  2:905448/1‑93 (MQ=255)
                                                                                      gcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGt                                                                                                                       <  1:905448/93‑1 (MQ=255)
                                                                                      gcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtctc                                                               <  4:1886685/149‑1 (MQ=255)
                                                                                       cAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCg                                                              <  1:1387808/149‑1 (MQ=255)
                                                                                       cAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTCCTCg                                                              <  2:1136163/149‑1 (MQ=255)
                                                                                        aaaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGa                                                                                            >  2:1739793/1‑118 (MQ=255)
                                                                                          aaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg                                                           >  3:1049901/1‑149 (MQ=255)
                                                                                          aaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg                                                           >  4:1394295/1‑149 (MQ=255)
                                                                                          aaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg                                                           >  1:1788201/1‑149 (MQ=255)
                                                                                          aaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTg                                                           >  1:1863610/1‑149 (MQ=255)
                                                                                           aTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGa                                                          <  2:117008/149‑1 (MQ=255)
                                                                                                 ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGa                                                                                            >  2:1641899/1‑109 (MQ=255)
                                                                                                 ccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGa                                                    >  2:1767649/1‑149 (MQ=255)
                                                                                                      tGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGa                                                                                            >  4:962655/1‑104 (MQ=255)
                                                                                                      tGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATAc                                                       <  1:821652/141‑1 (MQ=255)
                                                                                                        cAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCg                                             <  2:340774/149‑1 (MQ=255)
                                                                                                         aGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACgg                                                                     >  2:1782743/1‑124 (MQ=255)
                                                                                                           aTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGa                                                                                   <  4:1517710/108‑1 (MQ=255)
                                                                                                           aTCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAActct                                                                             >  2:990806/1‑114 (MQ=255)
                                                                                                              aaTGCTCGAGCCTGCCCCGTACGTATTGCTTGCATGCTGCAATTgcca                                                                                                                                           <  2:1845336/48‑3 (MQ=255)
                                                                                                              aaTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGAc                                                                                           >  1:1574410/1‑97 (MQ=255)
                                                                                                              aaTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGcc                                       >  4:1787363/1‑149 (MQ=255)
                                                                                                              aaTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGcc                                       >  4:398542/1‑149 (MQ=255)
                                                                                                              aaTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGcc                                       >  3:1263864/1‑149 (MQ=255)
                                                                                                                tGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCt                                     <  4:160404/149‑1 (MQ=255)
                                                                                                                tGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCt                                     <  2:870657/149‑1 (MQ=255)
                                                                                                                   tCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCt                                     <  3:413822/146‑1 (MQ=255)
                                                                                                                   tCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTa                                  >  1:470605/1‑149 (MQ=255)
                                                                                                                      aGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACAt                               <  2:730738/149‑1 (MQ=255)
                                                                                                                       gCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACAtt                              <  2:1790830/149‑1 (MQ=255)
                                                                                                                        ccTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGTTACCGACGGTTCGGTCGCCCTGTACATTa                             <  3:1208364/149‑1 (MQ=255)
                                                                                                                          tGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                           <  4:1313970/149‑1 (MQ=255)
                                                                                                                          tGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                           <  4:1049901/149‑1 (MQ=255)
                                                                                                                               cGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTc                      <  3:1349374/149‑1 (MQ=255)
                                                                                                                                gTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACTCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTcc                     >  1:848779/1‑149 (MQ=255)
                                                                                                                                gTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTcc                     <  4:1263864/149‑1 (MQ=255)
                                                                                                                                gTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTcc                     >  3:663090/1‑149 (MQ=255)
                                                                                                                                  aCATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCag                   <  3:1394295/149‑2 (MQ=255)
                                                                                                                                    atatTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTagag                          >  4:410367/1‑140 (MQ=255)
                                                                                                                                    atatTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGt                 >  2:1554834/1‑149 (MQ=255)
                                                                                                                                       ttgcttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGt              >  4:720456/1‑149 (MQ=255)
                                                                                                                                        tgcttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTg             >  3:623383/1‑149 (MQ=255)
                                                                                                                                         gcttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTgg            >  3:1782262/1‑149 (MQ=255)
                                                                                                                                           ttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTggg                                                                                                                    >  4:852895/1‑43 (MQ=255)
                                                                                                                                           ttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTggg                                                                                                                    <  3:852895/43‑1 (MQ=255)
                                                                                                                                           ttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGa                                                          >  4:1136653/1‑101 (MQ=255)
                                                                                                                                           ttgcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGa                                                          <  3:1136653/101‑1 (MQ=255)
                                                                                                                                            gacATGCTGCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGccc                                      >  3:1873535/3‑120 (MQ=255)
                                                                                                                                             gcATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAg        <  1:58870/149‑1 (MQ=255)
                                                                                                                                             acATGCTGCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAg        >  3:1184367/2‑149 (MQ=255)
                                                                                                                                              cATGCTGCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGATGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGa       >  2:790452/1‑149 (MQ=255)
                                                                                                                                              cATGCTGCAATTGCCAGCATGTCCAGGTACATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTtga                           <  4:627579/129‑1 (MQ=255)
                                                                                                                                               aTGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGAt      >  4:1718895/1‑149 (MQ=255)
                                                                                                                                                   tGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATAccc  >  2:1202005/1‑149 (MQ=255)
                                                                                                                                                   tGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTGCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATAccc  >  4:424468/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGCGGTTAAGCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTGCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCC  >  NC_000913/1585001‑1585296

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