Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,983,761 C→T 100% E437K (GAA→AAA)  araG ← L‑arabinose ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,983,7610CT99.2% 407.1 / ‑4.6 120E437K (GAA→AAA) araGL‑arabinose ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (61/58);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

Only 100 of 125 total aligned reads displayed.

TCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCTG  >  NC_000913/1983665‑1983857
                                                                                                |                                                                                                
tCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTtgtg                                                                                               <  1:2762239/100‑1 (MQ=255)
 ccGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGc                                                                                              >  1:2852158/1‑100 (MQ=255)
  cGGTCGTCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCTCCACGCCCCGCGCCGCCAGCGCATAAATTACGTTATGTCTTTTGTGCt                                                                                             >  1:115541/1‑100 (MQ=255)
      cGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgc                                                                                         >  1:1493139/1‑100 (MQ=255)
       ggCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcg                                                                                        <  2:1945187/100‑1 (MQ=255)
        gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCAGCGCCACGCCCTGCGCCGCCAGCGCATATATTACGTTCTATATTTTGTGCTGAgcgc                                                                                       >  1:2657720/1‑100 (MQ=255)
        gCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAgcgc                                                                                       >  2:2215220/1‑100 (MQ=255)
         cAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGcc                                                                                      <  2:868480/100‑1 (MQ=255)
         cAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGcc                                                                                      >  1:1958827/1‑100 (MQ=255)
            cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                   <  1:3576889/100‑1 (MQ=255)
            cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                   <  1:151332/100‑1 (MQ=255)
            cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                   <  1:1188563/100‑1 (MQ=255)
            cGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAAc                                                                                   >  2:135113/1‑100 (MQ=255)
             gCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACa                                                                                  <  2:2754334/100‑1 (MQ=255)
              ccGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACAt                                                                                 <  2:1127050/100‑1 (MQ=255)
              ccGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACAt                                                                                 >  2:3500018/1‑100 (MQ=255)
               cGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATc                                                                                <  2:2454717/100‑1 (MQ=255)
               cGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATc                                                                                <  1:2773027/100‑1 (MQ=255)
                gAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGTCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCa                                                                               <  1:2373599/100‑1 (MQ=255)
                 aGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCaa                                                                              <  2:4211928/100‑1 (MQ=255)
                 aGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGACCTGTGCCGCCAGCGCATAAATTACGTTATACATTTTGTGCTTAGCGGCAAAATc                                                                                <  1:4265976/98‑1 (MQ=255)
                  ggACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAAt                                                                             <  1:2862269/100‑1 (MQ=255)
                   gACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGTGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATg                                                                            <  1:707013/100‑1 (MQ=255)
                      ttCAGGTAAGTCGCTGGAGGCAAACAGCAACGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcg                                                                         <  1:2485646/100‑1 (MQ=255)
                       tCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgc                                                                        >  2:1614907/1‑100 (MQ=255)
                       tCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgc                                                                        >  1:681532/1‑100 (MQ=255)
                       tCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgc                                                                        >  2:735254/1‑100 (MQ=255)
                        cAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcg                                                                       <  1:1329157/100‑1 (MQ=255)
                         aGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgc                                                                      >  1:1382131/1‑100 (MQ=255)
                          ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCGGCGCATAAATTGCGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                     <  2:1240350/100‑1 (MQ=255)
                          ggTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCcgcgcg                                                                     <  1:1409014/100‑1 (MQ=255)
                           gTAAGTCGCTGGAGGCAAACAGCACTGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGt                                                                    >  2:495240/1‑100 (MQ=255)
                            tAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTa                                                                   >  2:3976451/1‑100 (MQ=255)
                            tAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTa                                                                   <  2:1445031/100‑1 (MQ=255)
                             aaGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAg                                                                  >  2:2225955/1‑100 (MQ=255)
                              aGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATGTATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAgg                                                                 >  1:110588/1‑100 (MQ=255)
                               gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                                >  1:494585/1‑100 (MQ=255)
                               gTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGt                                                                <  1:3278818/100‑1 (MQ=255)
                                tCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGtt                                                               <  2:308136/100‑1 (MQ=255)
                                tCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGtt                                                               <  1:4213527/100‑1 (MQ=255)
                                  gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                             >  1:3970656/1‑100 (MQ=255)
                                  gCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCa                                                             >  2:3425079/1‑100 (MQ=255)
                                    tGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATc                                                           <  1:1977331/100‑1 (MQ=255)
                                    tGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATc                                                           <  1:503926/100‑1 (MQ=255)
                                    tGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATc                                                           <  2:750955/100‑1 (MQ=255)
                                    tGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATc                                                           >  2:1173997/1‑100 (MQ=255)
                                     ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                          <  1:2940386/100‑1 (MQ=255)
                                     ggAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATcc                                                          >  2:2671678/1‑100 (MQ=255)
                                      gAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCa                                                         <  2:3469428/100‑1 (MQ=255)
                                        ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGc                                                       >  2:257255/1‑100 (MQ=255)
                                        ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGc                                                       <  1:3257404/100‑1 (MQ=255)
                                        ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGc                                                       <  2:669754/100‑1 (MQ=255)
                                        ggCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGAGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGc                                                       >  1:796148/1‑100 (MQ=255)
                                         gCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCa                                                      >  2:1116143/1‑100 (MQ=255)
                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                     >  1:3488353/1‑100 (MQ=255)
                                          cAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaa                                                     <  1:3582826/100‑1 (MQ=255)
                                            aaCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCaaaa                                                   <  2:3899161/100‑1 (MQ=255)
                                             aCAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAAt                                                  <  1:1779429/100‑1 (MQ=255)
                                              cAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATg                                                 >  2:3041620/1‑100 (MQ=255)
                                              cAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATg                                                 >  2:3118511/1‑100 (MQ=255)
                                               aGCACCGCCACGCCCTGCTCCGCCAGCTCATAAATTGCGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGa                                                <  2:3997542/100‑1 (MQ=255)
                                               aGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCTGCGCGTAGGTTCATCCAGCAAAATGa                                                <  2:2603874/100‑1 (MQ=255)
                                                 cACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGAcc                                              >  2:1363960/1‑100 (MQ=255)
                                                   ccGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCtt                                            >  1:32744/1‑100 (MQ=255)
                                                       cACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                        <  2:3852222/100‑1 (MQ=255)
                                                       cACGCCCTGCGCCGCCAGCGCATAAAGTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCAtc                                        <  1:1407471/100‑1 (MQ=255)
                                                          gCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATctct                                     >  1:1502619/1‑100 (MQ=255)
                                                            ccTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTc                                   >  1:2835794/1‑100 (MQ=255)
                                                             cTGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTcc                                  <  1:3950790/100‑1 (MQ=255)
                                                              tGCGCCGCCAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCg                                 <  1:4144608/100‑1 (MQ=255)
                                                                   cgccAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAAc                            >  1:549924/1‑100 (MQ=255)
                                                                      cAGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAg                         <  2:1635576/100‑1 (MQ=255)
                                                                       aGCGCATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGc                        >  2:432785/1‑100 (MQ=255)
                                                                        gcgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCg                       >  2:2726668/1‑100 (MQ=255)
                                                                        gcgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATAAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGc                        >  1:1779817/1‑99 (MQ=255)
                                                                         cgcATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCgg                      >  1:863703/1‑100 (MQ=255)
                                                                           cATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGcc                    <  1:2156423/100‑1 (MQ=255)
                                                                           cATAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGcc                    <  1:2669337/100‑1 (MQ=255)
                                                                             tAAATTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCa                  >  1:1707816/1‑100 (MQ=255)
                                                                             tAAATCACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAACATGACCTTCATCCCTTCCGATAACCATCGGCCCa                  >  1:1363686/1‑100 (MQ=255)
                                                                               aaTTACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGa                >  1:4067850/1‑100 (MQ=255)
                                                                                 ttACGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAAt              <  1:2298986/100‑1 (MQ=255)
                                                                                   aCGTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATgg            >  2:3045938/1‑100 (MQ=255)
                                                                                     gTTATATATTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCt          <  2:2569702/100‑1 (MQ=255)
                                                                                       tatataTTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttt        >  2:2217283/1‑100 (MQ=255)
                                                                                        atatatTTTGTGCTTAGCGCCAACATCAATGCCACGCGTAGGTTCATACAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCtttt       >  1:2078998/1‑100 (MQ=255)
                                                                                         tatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttttt      >  1:652885/1‑100 (MQ=255)
                                                                                         tatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttttt      >  1:3414381/1‑100 (MQ=255)
                                                                                         tatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttttt      >  2:1100110/1‑100 (MQ=255)
                                                                                         tatatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCttttt      >  1:3406950/1‑100 (MQ=255)
                                                                                          atatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTg     <  1:2090/100‑1 (MQ=255)
                                                                                           tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc    <  2:2830686/100‑1 (MQ=255)
                                                                                           tatTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGc    <  2:1126684/100‑1 (MQ=255)
                                                                                            atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt   >  2:1754061/1‑100 (MQ=255)
                                                                                            atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt   >  2:734554/1‑100 (MQ=255)
                                                                                            atTTTGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCt   >  1:1332702/1‑100 (MQ=255)
                                                                                             ttttGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATggagatcggaag  >  1:868773/1‑90 (MQ=255)
                                                                                             ttttGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCTg  >  1:3639120/1‑100 (MQ=255)
                                                                                             ttttGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCTg  >  1:1010729/1‑100 (MQ=255)
                                                                                             tttCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCTg  <  2:3554130/100‑1 (MQ=255)
                                                                                                |                                                                                                
TCCGGTCGGCAACGCCGAGGACTTCAGGTAAGTCGCTGGAGGCAAACAGCACCGCCACGCCCTGCGCCGCCAGCGCATAAATTACGTTATATATTTCGTGCTTAGCGCCAACATCAATGCCGCGCGTAGGTTCATCCAGCAAAATGACCTTCATCTCTTCCGATAACCAGCGGCCCAGAATGGCTTTTTGCTG  >  NC_000913/1983665‑1983857

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