Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 452,478 C→T 100% L356L (CTG→CTA)  ampG ← muropeptide transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913452,4780CT88.2% 267.6 / 26.0 112L356L (CTG→CTA) ampGmuropeptide transporter
Reads supporting (aligned to +/- strand):  ref base C (9/4);  new base T (48/50);  total (58/54)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.70e-01

Only 100 of 113 total aligned reads displayed.

CGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAG  >  NC_000913/452381‑452577
                                                                                                 |                                                                                                   
 ggAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAgg                                                                                                  >  1:1141547/1‑100 (MQ=255)
 ggAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAgg                                                                                                  <  2:512709/100‑1 (MQ=255)
    gaATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCa                                                                                               >  1:2054638/1‑100 (MQ=255)
      atagatagAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCaaa                                                                                             >  2:2563029/1‑100 (MQ=255)
         gatagaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTg                                                                                          >  2:1950864/1‑100 (MQ=255)
         cttagaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTg                                                                                          <  1:1186829/98‑1 (MQ=255)
         cttagaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTg                                                                                          <  1:431366/98‑1 (MQ=255)
          aaagaATGTCGACCAGCCCTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCTCTCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTg                                                                                          <  1:2088601/97‑1 (MQ=255)
           tagaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAg                                                                                        >  2:1186834/1‑100 (MQ=255)
           tagaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAg                                                                                        >  2:431368/1‑100 (MQ=255)
            agaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGt                                                                                       >  1:2304396/1‑100 (MQ=255)
            agaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGt                                                                                       <  2:624691/100‑1 (MQ=255)
             gaATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGt                                                                                       >  1:2182601/1‑99 (MQ=255)
                 gTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGTTGAGAGTAGGGCAAATTGAGTAGCgg                                                                                  <  2:1772566/100‑1 (MQ=255)
                 gTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCg                                                                                   <  1:1975942/99‑1 (MQ=255)
                   cGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGaa                                                                                <  2:2737518/100‑1 (MQ=255)
                   cGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGaa                                                                                <  2:2744876/100‑1 (MQ=255)
                    gACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGaaa                                                                               >  2:2207600/1‑100 (MQ=255)
                     aCCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGaaaa                                                                              <  2:1405445/100‑1 (MQ=255)
                      ccAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAATCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAAt                                                                             >  2:132174/1‑100 (MQ=255)
                         gCCGTGTGCTTCAACAAACCTACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGAc                                                                          <  2:1878936/100‑1 (MQ=255)
                         gCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGAc                                                                          >  2:1169238/1‑100 (MQ=255)
                         gCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGAc                                                                          >  1:51597/1‑100 (MQ=255)
                          ccGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACt                                                                         >  2:448301/1‑100 (MQ=255)
                            gtgtGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACtta                                                                       >  1:151861/1‑100 (MQ=255)
                               tgCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACttatta                                                                    <  1:600091/100‑1 (MQ=255)
                               tgCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACttatta                                                                    <  1:364517/100‑1 (MQ=255)
                               tgCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACttatta                                                                    <  2:1929821/100‑1 (MQ=255)
                               tgCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACttatta                                                                    <  1:1950855/100‑1 (MQ=255)
                               tgCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGAAGGGCAAAATGAGAAGCGGAAAATAAAttataa                                                                    >  1:228823/1‑98 (MQ=255)
                               tgCTTCAACAAACCAACCCGCCACCGCGCCAACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACttatt                                                                     <  1:2057190/99‑1 (MQ=255)
                                gCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACttatta                                                                    <  1:81517/99‑1 (MQ=255)
                                gCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTAc                                                                   >  2:2017751/1‑100 (MQ=255)
                                 cTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTAc                                                                   <  1:1864579/99‑1 (MQ=255)
                                 cTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACa                                                                  >  1:2575588/1‑100 (MQ=255)
                                 cTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACa                                                                  >  2:765284/1‑100 (MQ=255)
                                 cTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACa                                                                  >  2:1095136/1‑100 (MQ=255)
                                 cTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACa                                                                  >  2:499992/1‑100 (MQ=255)
                                 cTTCAACAAACCAACCCGCCACGGCGCCGACATAAACTCGCCCTACAGCAGAAAGCGCCGAGAGCACGGCAAATTCAGTAGCGGAAAACGACTTATTACa                                                                  >  1:38765/1‑100 (MQ=255)
                                  ttCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACAt                                                                 >  1:486484/1‑100 (MQ=255)
                                   tCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATa                                                                >  2:1154243/1‑100 (MQ=255)
                                    caacaaACCAACCCGCCACGGGGCCGACATAAAGTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACCTATTACATAg                                                               >  1:1915716/1‑100 (MQ=255)
                                    caacaaACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAg                                                               >  1:1648948/1‑100 (MQ=255)
                                    caacaaACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAg                                                               >  2:2585188/1‑100 (MQ=255)
                                      acaaACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCg                                                             <  2:399690/100‑1 (MQ=255)
                                      acaaACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCg                                                             >  2:20671/1‑100 (MQ=255)
                                        aaaCCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTc                                                           <  2:485319/100‑1 (MQ=255)
                                         aaccaaccCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCa                                                          >  1:1661283/1‑100 (MQ=255)
                                           ccaaccCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCa                                                          <  1:2773624/98‑1 (MQ=255)
                                           ccaaccCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCAt                                                         <  2:568082/99‑1 (MQ=255)
                                              accCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAAc                                                     <  2:283900/100‑1 (MQ=255)
                                               cccGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACa                                                    <  1:1845555/100‑1 (MQ=255)
                                               cccGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACa                                                    <  2:274815/100‑1 (MQ=255)
                                                ccGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAg                                                   <  2:2544039/100‑1 (MQ=255)
                                                  gCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAgcg                                                 <  1:2105217/100‑1 (MQ=255)
                                                  gCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAgcg                                                 <  2:2066921/100‑1 (MQ=255)
                                                  gCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAgcg                                                 <  2:745665/100‑1 (MQ=255)
                                                   ccACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAgcgc                                                <  2:65109/100‑1 (MQ=255)
                                                   ccACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAgcgc                                                <  2:1988423/100‑1 (MQ=255)
                                                    cACGTGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGcgcg                                               <  2:2131508/100‑1 (MQ=255)
                                                          gCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAg                                         <  2:702112/100‑1 (MQ=255)
                                                            cGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGc                                       <  1:1495520/100‑1 (MQ=255)
                                                             gACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCt                                      <  1:1319347/100‑1 (MQ=255)
                                                              aCATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTg                                     >  2:175835/1‑100 (MQ=255)
                                                              aCATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTg                                     >  2:357507/1‑100 (MQ=255)
                                                                  aaaCTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGAtgt                                 >  1:2858341/1‑100 (MQ=255)
                                                                   aaCTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGAtgtg                                <  2:1757247/100‑1 (MQ=255)
                                                                   aaCTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGAtgtg                                >  2:123787/1‑100 (MQ=255)
                                                                    gcTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGc                               >  1:2658551/2‑100 (MQ=255)
                                                                    aCTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGc                               >  2:771986/1‑100 (MQ=255)
                                                                    aCTCGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGc                               >  2:1359201/1‑100 (MQ=255)
                                                                       cGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGcc                              <  1:356718/98‑1 (MQ=255)
                                                                       cGCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCa                            <  2:2868578/100‑1 (MQ=255)
                                                                       cGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCa                            <  2:2861254/100‑1 (MQ=255)
                                                                        gCCCTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCAt                           <  1:2136223/100‑1 (MQ=255)
                                                                        gCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCAt                           >  2:2380412/1‑100 (MQ=255)
                                                                          ccTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTAAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATcc                         >  1:1385287/1‑100 (MQ=255)
                                                                          ccTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATc                          >  1:2533200/1‑99 (MQ=255)
                                                                           cTACAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATccc                        >  2:1345808/1‑100 (MQ=255)
                                                                             aCAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGc                      >  2:1511538/1‑100 (MQ=255)
                                                                             aCAGCAGAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGc                      <  1:872738/100‑1 (MQ=255)
                                                                                gcagAAAGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCcac                   <  2:1801539/100‑1 (MQ=255)
                                                                                gcagAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCcac                   >  1:1291355/1‑100 (MQ=255)
                                                                                    aaaGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACagag               <  2:1062662/100‑1 (MQ=255)
                                                                                     aaGCGCTGAGAGTATGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAgg              <  1:709316/100‑1 (MQ=255)
                                                                                     aaGCGCTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAgg              <  2:2485874/100‑1 (MQ=255)
                                                                                     aaGCGCGGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAgg              <  2:1058988/100‑1 (MQ=255)
                                                                                     aaGCGCAGAGAGTAAGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACACCGCGACATAGGCTGATGGGCCCTCCCCTCCACagacg              <  1:2345323/100‑3 (MQ=255)
                                                                                        cgcTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGttt           >  1:604352/1‑100 (MQ=255)
                                                                                          cTGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTc         >  2:271365/1‑100 (MQ=255)
                                                                                           tGAGAGTAGGGCAAATTGAGTAGCGGAACATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCg        <  1:1228611/100‑1 (MQ=255)
                                                                                           tGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCg        <  2:1574352/100‑1 (MQ=255)
                                                                                           tGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCg        >  1:1940441/1‑100 (MQ=255)
                                                                                           tGAGAGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCg        <  1:2555026/100‑1 (MQ=255)
                                                                                            gagaGTAGGGCAAATTGAGTAGCGGAACATGACTTATTACATAGCGTCAGTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTc         <  1:1908531/98‑1 (MQ=255)
                                                                                            gagaGTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGa       >  2:922155/1‑100 (MQ=255)
                                                                                              gagTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGaaa     >  1:1907565/1‑100 (MQ=255)
                                                                                              gagTAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGaaa     >  2:1789525/1‑100 (MQ=255)
                                                                                              gagCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGaaa     <  2:1946140/100‑1 (MQ=255)
                                                                                                 cAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAg  >  2:2503854/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
CGGAGAATAGATAGAATGTCGACCAGCCGTGTGCTTCAACAAACCAACCCGCCACGGGGCCGACATAAACTCGCCCTACAGCAGAAAGCGCTGAGAGCAGGGCAAATTGAGTAGCGGAAAATGACTTATTACATAGCGTCATTAACAGCGCGACAAAGGCTGATGTGCCCATCCCGCCACAGAGGTTTTCGAAAAAG  >  NC_000913/452381‑452577

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