Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 20,724 G→A intergenic (‑216/+91) insA‑1 ← / ← rpsT IS1 protein InsA/30S ribosomal subunit protein S20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091320,7240GA100.0% 1432.2 / NA 335intergenic (‑216/+91)insA‑1/rpsTIS1 protein InsA/30S ribosomal subunit protein S20
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (186/149);  total (186/149)

Only 100 of 341 total aligned reads displayed.

TTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCC  >  NC_000913/20626‑20821
                                                                                                  |                                                                                                 
ttCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGc                                                                                                 <  1:38606/101‑1 (MQ=255)
 tCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCa                                                                                                >  1:39018/1‑101 (MQ=255)
  cTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACcagcag                                                                                               <  1:37580/101‑1 (MQ=255)
  cTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACcagcag                                                                                               >  2:39020/1‑101 (MQ=255)
    tGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAgtg                                                                                             >  2:39025/1‑101 (MQ=255)
     gTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAgtgt                                                                                            >  1:39028/1‑101 (MQ=255)
     gTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAgtgt                                                                                            <  1:38269/101‑1 (MQ=255)
      ttCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAgtgtg                                                                                           >  2:39032/1‑101 (MQ=255)
       tCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAgtgtgt                                                                                          >  1:39034/1‑101 (MQ=255)
       tCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAgtgtgt                                                                                          <  2:38739/101‑1 (MQ=255)
        cAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGtt                                                                                         >  2:39040/1‑101 (MQ=255)
         aGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTtc                                                                                        <  1:38375/101‑1 (MQ=255)
          gAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTtct                                                                                       <  2:37941/101‑1 (MQ=255)
            aTGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTtctct                                                                                     <  1:38301/101‑1 (MQ=255)
            aTGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTtctct                                                                                     >  1:39044/1‑101 (MQ=255)
              gTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTtctctc                                                                                    >  1:39048/1‑100 (MQ=255)
              gTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCa                                                                                   <  2:38486/101‑1 (MQ=255)
                 tCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGc                                                                                >  2:39055/1‑101 (MQ=255)
                  cAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACACCGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCa                                                                               >  1:39056/1‑101 (MQ=255)
                       gATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATcaca                                                                          >  1:39064/1‑101 (MQ=255)
                       gATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATcaaa                                                                          <  1:38815/101‑3 (MQ=255)
                        aTGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACaa                                                                         >  1:39067/1‑101 (MQ=255)
                         tGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACaaa                                                                        <  1:38403/101‑1 (MQ=255)
                         tGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACaaa                                                                        >  2:39068/1‑101 (MQ=255)
                          ggTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACaaaa                                                                       <  2:38562/101‑1 (MQ=255)
                          ggTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACaaaa                                                                       <  1:38237/101‑1 (MQ=255)
                           gTTCTACGATTCTTAAGCCGCGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAag                                                                      <  2:38043/101‑1 (MQ=255)
                           gTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAag                                                                      <  1:38646/101‑1 (MQ=255)
                            ttCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAagc                                                                     >  1:39073/1‑101 (MQ=255)
                             tCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAagca                                                                    >  2:39076/1‑101 (MQ=255)
                               tACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAagcagc                                                                  <  1:37923/101‑1 (MQ=255)
                                aCGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAgtg                                                                                             >  2:39081/1‑73 (MQ=255)
                                  gATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAgg                                                               >  1:39084/1‑101 (MQ=255)
                                  gATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAgg                                                               <  2:38756/101‑1 (MQ=255)
                                   aTTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGc                                                              >  2:39086/1‑101 (MQ=255)
                                    ttCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCa                                                             <  2:38335/101‑1 (MQ=255)
                                    ttCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCa                                                             <  2:38278/101‑1 (MQ=255)
                                      cTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATa                                                           >  1:39089/1‑101 (MQ=255)
                                       ttAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaa                                                          >  2:39093/1‑101 (MQ=255)
                                        tAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaa                                                         <  1:38891/101‑1 (MQ=255)
                                        tAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaa                                                         <  1:38128/101‑1 (MQ=255)
                                         aaGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaaa                                                        >  1:39098/1‑101 (MQ=255)
                                         aaGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaaa                                                        >  2:39099/1‑101 (MQ=255)
                                           gCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaaaaa                                                      >  1:39103/1‑101 (MQ=255)
                                            ccACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAagca                                                                    <  2:38390/86‑1 (MQ=255)
                                            ccACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaaaaaa                                                     >  2:39105/1‑101 (MQ=255)
                                             cACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaaaa                                                       >  2:39109/1‑98 (MQ=255)
                                             cACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAAc                                                    <  2:38304/101‑1 (MQ=255)
                                                 aaGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaaaaa                                                      >  2:39116/1‑95 (MQ=255)
                                                 aaGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGc                                                <  2:38788/101‑1 (MQ=255)
                                                 aaGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGc                                                >  1:39115/1‑101 (MQ=255)
                                                   gagTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCtt                                              >  1:39119/1‑101 (MQ=255)
                                                      ttCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCtt                                              <  2:38647/98‑1 (MQ=255)
                                                      ttCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTgcg                                           >  1:39123/1‑101 (MQ=255)
                                                       tCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTgcgc                                          >  1:39127/1‑101 (MQ=255)
                                                        cAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGcgcg                                         >  1:39130/1‑101 (MQ=255)
                                                        cAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGcgcg                                         <  2:38697/101‑1 (MQ=255)
                                                         aGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCgg                                        <  2:38836/101‑1 (MQ=255)
                                                           aTAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGc                                      <  2:38522/101‑1 (MQ=255)
                                                            tAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCt                                     <  2:39047/101‑1 (MQ=255)
                                                            tAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCt                                     <  2:38240/101‑1 (MQ=255)
                                                            tAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCt                                     >  1:39135/1‑101 (MQ=255)
                                                             aGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCtt                                    <  1:38034/101‑1 (MQ=255)
                                                             aGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCtt                                    >  2:39137/1‑101 (MQ=255)
                                                               tACAACGGCATGTCTCTTTTGACTATCTGGCAACCAGc                                                                                                 >  2:39141/1‑38 (MQ=255)
                                                                aCAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCttttt                                 <  1:38550/101‑1 (MQ=255)
                                                                 cAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTc                                <  2:38191/101‑1 (MQ=255)
                                                                 cAACGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTc                                <  1:38215/101‑1 (MQ=255)
                                                                  aaCGGCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCttttt                                 >  2:39144/1‑99 (MQ=255)
                                                                     ggCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACaa                            <  1:38483/101‑1 (MQ=255)
                                                                      gCATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCGGCTTGCGCGGGCTTTTTCACaaa                           >  2:39149/1‑101 (MQ=255)
                                                                       cATGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAg                          >  1:39151/1‑101 (MQ=255)
                                                                        aTGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGc                         >  1:39156/1‑101 (MQ=255)
                                                                         tGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTgcgc                                          >  2:39160/1‑83 (MQ=255)
                                                                         tGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGcgcg                                         >  1:39159/1‑84 (MQ=255)
                                                                         tGTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCt                        <  1:38177/101‑1 (MQ=255)
                                                                          gTCTCTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCtt                       <  1:38744/101‑1 (MQ=255)
                                                                           tctcTTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTc                      <  2:38438/101‑1 (MQ=255)
                                                                            ctctTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCa                     <  1:38244/101‑1 (MQ=255)
                                                                             tctTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAg                    >  1:39164/1‑101 (MQ=255)
                                                                             tctTTTGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAg                    <  1:38696/101‑1 (MQ=255)
                                                                                tttGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATaaaaaa                                                      >  2:39169/1‑64 (MQ=255)
                                                                                 ttGACTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCaaa                >  1:39171/1‑101 (MQ=255)
                                                                                    aCTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTg             <  1:38455/101‑1 (MQ=255)
                                                                                     cTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTgg            <  2:38625/101‑1 (MQ=255)
                                                                                     cTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTgg            >  2:39175/1‑101 (MQ=255)
                                                                                     cTATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTgg            >  1:39177/1‑101 (MQ=255)
                                                                                      tATCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGc           >  2:39182/1‑101 (MQ=255)
                                                                                       aTCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCg          >  2:39187/1‑101 (MQ=255)
                                                                                       aTCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCg          >  1:39188/1‑101 (MQ=255)
                                                                                       aTCTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCg          <  1:38436/101‑1 (MQ=255)
                                                                                         cTGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAAtt              >  2:39193/1‑95 (MQ=255)
                                                                                          tGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGAtt       >  2:39198/1‑101 (MQ=255)
                                                                                          tGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGAtt       >  1:39197/1‑101 (MQ=255)
                                                                                          tGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGAtt       <  2:38156/101‑1 (MQ=255)
                                                                                          tGGCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGAtt       <  1:38511/101‑1 (MQ=255)
                                                                                            gCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTaa     >  2:39203/1‑101 (MQ=255)
                                                                                            gCAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCGTTTTCACAAAGCTTCAGCAAATTGGCGATTaa     <  2:38132/101‑1 (MQ=255)
                                                                                             cAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAg    >  1:39204/1‑101 (MQ=255)
                                                                                             cAACCAGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAg    <  1:38588/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
TTCTTGTTCAGAATGTATCAGCCGATGGTTCTACGATTCTTAAGCCACGAAGAGTTCAGATAGTACAACGGCATGTCTCTTTTGACTATCTGGCAACCGGCAGTGTGTTCTCTCACGCATCACAAAAGCAGCAGGCATAAAAAAACCCGCTTGCGCGGGCTTTTTCACAAAGCTTCAGCAAATTGGCGATTAAGCC  >  NC_000913/20626‑20821

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