Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,655,922 C→T 16.2% intergenic (‑20/+44) yhiD ← / ← hdeB putative Mg(2+) transport ATPase, inner membrane protein/acid‑resistance protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,655,9220CT16.2% 449.3 / 61.7 178intergenic (‑20/+44)yhiD/hdeBputative Mg(2+) transport ATPase, inner membrane protein/acid‑resistance protein
Reads supporting (aligned to +/- strand):  ref base C (77/72);  new base T (15/14);  total (92/86)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.88e-01

Only 100 of 187 total aligned reads displayed.

CCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGC  >  NC_000913/3655823‑3656020
                                                                                                   |                                                                                                   
ccGCGCATTTGCCTTTCCATACCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGagc                                                                                                     >  2:2340344/1‑100 (MQ=255)
 cgcgCATTTGCCTTTCCATGCCAATGGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGagca                                                                                                    >  2:1462574/1‑100 (MQ=255)
 cgcgCATTTGCCTTTCCATACCAATAGCGCCACAGGCAATTGCCGCGAGTATAATACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGagc                                                                                                     <  2:765294/99‑1 (MQ=255)
  gcgcATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGagcag                                                                                                   >  2:884082/1‑100 (MQ=255)
   cgcATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGagcagc                                                                                                  <  1:2021313/100‑1 (MQ=255)
    gcATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAgcagca                                                                                                 >  2:3010881/1‑100 (MQ=255)
         tGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATg                                                                                            <  1:1993321/100‑1 (MQ=255)
          gCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATgg                                                                                          >  1:202838/1‑100 (MQ=255)
            cTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGTAGCAGCAAGATGGCt                                                                                        >  1:3757890/1‑100 (MQ=255)
            cTTTCCATGCAAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAAAAAAAAAAGGAGCAGCAAGATGGc                                                                                         <  2:198662/99‑1 (MQ=255)
             tttCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTc                                                                                       <  1:2802250/100‑1 (MQ=255)
                ccATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAAc                                                                                    >  1:1539154/1‑100 (MQ=255)
                 cATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACt                                                                                   >  1:2090167/1‑100 (MQ=255)
                 cATGCAAATAGCGCCACTGTCAATTGCCGCGAGTATAATACGTATAATAAATTCTGCTGTCAAAAAAAAAAAAAAAAAGGAGCAGCAAGATGGCTCAACt                                                                                   <  2:395200/100‑1 (MQ=255)
                  aTGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACtt                                                                                  <  1:597159/100‑1 (MQ=255)
                     ccAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCt                                                                               <  2:996310/100‑1 (MQ=255)
                     ccAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCt                                                                               <  1:1530793/100‑1 (MQ=255)
                      cAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCt                                                                               <  2:3333587/99‑1 (MQ=255)
                       aaTAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCAAAAAAAAAAAAAAAAAGGAGCAGCAAGATGGCTCAACTTGCt                                                                               <  2:2690147/98‑1 (MQ=255)
                             gcaacaGGCAATTGCCGTGAGTATAAGATGTATAATAAATTCTGCTGTTATAAAAAAAAAAAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCttt                                                                       <  2:1431764/97‑1 (MQ=255)
                             gcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCttt                                                                       >  1:781185/1‑100 (MQ=255)
                             gcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCttt                                                                       >  1:1131642/1‑100 (MQ=255)
                             gcCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAACATGGCTCAACTTGCTACTc                                                                           >  2:1889628/1‑96 (MQ=255)
                              ccACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCtttt                                                                      <  1:2492719/100‑1 (MQ=255)
                               cacaGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCttttt                                                                     <  1:332756/100‑1 (MQ=21)
                               cacaGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCttttt                                                                     <  1:50838/100‑1 (MQ=255)
                                acaGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTa                                                                    >  1:1288996/1‑100 (MQ=255)
                                acaGCCAAGTCCCGCAAGGAGAATACTTATAATAAATTCTTCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttg                                                                    <  1:1651830/100‑2 (MQ=2)
                                 caGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttgc                                                                   <  1:942026/100‑3 (MQ=18)
                                  aGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttgca                                                                  <  2:2923880/100‑4 (MQ=18)
                                  aGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttgca                                                                  <  2:195499/100‑4 (MQ=18)
                                  aGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCtttttgca                                                                  <  1:310429/100‑4 (MQ=21)
                                  aGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACt                                                                  >  1:1696796/1‑100 (MQ=255)
                                  aGGCAATTGCCAAAAGAATAAGACGTATAATAAAGGCTGCTGGCATATACAAAGGCAAAAATTAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTAc                                                                   <  2:3062446/99‑1 (MQ=255)
                                   ggCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACtt                                                                 >  2:2943627/1‑100 (MQ=255)
                                    gaaaTTGCCGCGAGTAAAAGACGTATAAGAAATTCTGCTGTCATATAAAAATGGCAAAAGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCcttttgca                                                                  <  2:1869235/96‑4 (MQ=12)
                                    gCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTg                                                                <  1:3298301/100‑1 (MQ=255)
                                    gCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTg                                                                <  1:30870/100‑1 (MQ=255)
                                       aTTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCAc                                                             <  2:2102910/100‑1 (MQ=255)
                                       aTTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCAc                                                             <  1:661125/100‑1 (MQ=255)
                                          gCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGc                                                               <  2:2378993/95‑1 (MQ=255)
                                          gCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCt                                                           <  1:1203842/99‑1 (MQ=255)
                                          gCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTc                                                          >  2:2188887/1‑100 (MQ=255)
                                          gCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTc                                                          >  2:3477860/1‑100 (MQ=255)
                                           ccGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCa                                                         >  1:2425776/1‑100 (MQ=255)
                                            cgcgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCAt                                                        >  1:2115540/1‑100 (MQ=255)
                                            cgcgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCAt                                                        >  1:674267/1‑100 (MQ=255)
                                            cgcgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCAt                                                        >  1:2273900/1‑100 (MQ=255)
                                            cgcgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCAt                                                        >  2:1375218/1‑100 (MQ=255)
                                              cgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTa                                                      <  1:450972/100‑1 (MQ=255)
                                              cgAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTa                                                      <  1:290038/100‑1 (MQ=255)
                                                aGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAAt                                                    >  1:321346/1‑100 (MQ=255)
                                                 gTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAAtt                                                   <  1:2650301/100‑1 (MQ=255)
                                                  tataAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTc                                                  <  2:545455/100‑1 (MQ=255)
                                                   ataAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAAGAGCAGCAAGAGGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCg                                                 >  2:2440681/1‑100 (MQ=255)
                                                   ataAGACGTATAATAAATTCGGCTGTCATACAAAAATGCAAAAAGGATCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTc                                                  <  2:2425850/99‑1 (MQ=255)
                                                     aaGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGc                                               <  1:1074838/100‑1 (MQ=255)
                                                     aaGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGc                                               <  1:3384188/100‑1 (MQ=255)
                                                      aGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGc                                               >  1:65572/1‑99 (MQ=255)
                                                      aGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCa                                              >  2:1214863/1‑100 (MQ=255)
                                                      aGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCa                                              <  1:1023807/100‑1 (MQ=255)
                                                       gACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCaa                                             >  1:2827008/1‑100 (MQ=255)
                                                          gTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTc                                          >  1:1930578/1‑100 (MQ=255)
                                                            ataataAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCAt                                        >  1:495142/1‑100 (MQ=255)
                                                              aataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTa                                      >  1:2360052/1‑100 (MQ=255)
                                                              aataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTa                                      >  2:2289750/1‑100 (MQ=255)
                                                               ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAg                                     >  2:2754856/1‑100 (MQ=255)
                                                               ataaATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTACACCTCATTAATTCGGCAAGTCATTAg                                     >  1:555058/1‑100 (MQ=255)
                                                                 aaaTTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGAt                                   >  1:2921311/1‑100 (MQ=255)
                                                                   aTTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGc                                 >  1:1109828/1‑100 (MQ=255)
                                                                       tgctgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGa                             <  2:1971221/100‑1 (MQ=255)
                                                                         ctgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAtt                           >  1:3483698/1‑100 (MQ=255)
                                                                         ctgTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAtt                           >  2:3449983/1‑100 (MQ=255)
                                                                           gTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGAtttt                         <  2:1089479/100‑1 (MQ=255)
                                                                             cATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTCATTTTTg                       >  2:1061881/1‑100 (MQ=255)
                                                                              atatAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGa                      >  1:2201682/1‑100 (MQ=255)
                                                                               tataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                     >  2:3472537/1‑100 (MQ=255)
                                                                               tataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                     >  1:2619984/1‑100 (MQ=255)
                                                                               tataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGaa                     >  2:1699348/1‑100 (MQ=255)
                                                                                ataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAATTTGCTACTCCTTTTTACTTTCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAg                    <  2:3236724/100‑1 (MQ=255)
                                                                                ataAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAg                    <  1:2993519/100‑1 (MQ=255)
                                                                                 taAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAgg                   <  1:643152/100‑1 (MQ=255)
                                                                                 taAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAgg                   >  2:2325631/1‑100 (MQ=255)
                                                                                 taAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAgg                   <  1:2387654/100‑1 (MQ=255)
                                                                                 taAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAgg                   <  2:3102475/100‑1 (MQ=255)
                                                                                 taAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAgg                   <  1:692797/100‑1 (MQ=255)
                                                                                  gtgcaagtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGt                  >  2:2269509/9‑100 (MQ=17)
                                                                                   aaaaTGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGta                 >  2:3561499/1‑100 (MQ=255)
                                                                                      aagtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATAc              >  1:34539/5‑100 (MQ=18)
                                                                                       tGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACa             <  2:1926192/100‑1 (MQ=255)
                                                                                       agtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACa             >  1:1463584/4‑100 (MQ=18)
                                                                                        gtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTGCTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACaa            >  2:1436380/3‑100 (MQ=18)
                                                                                        gtaaaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTTAAGGTATACa             >  2:1814158/3‑99 (MQ=18)
                                                                                        gCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACaa            >  2:1180000/1‑100 (MQ=255)
                                                                                         taaaaaGGAGTAGCAAGTTGAGTCATCTTGCTGCTCCTTTTTACTTGAACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTAGAAGGTATACaa            >  1:3186363/2‑99 (MQ=18)
                                                                                           aaaaGGAGTAGCAAGTTGAGCCATCTTGCTGCTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACagatc         >  1:99838/1‑96 (MQ=18)
                                                                                           aaaaGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAAtt         >  1:868589/1‑100 (MQ=255)
                                                                                               ggAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTc     <  2:2619965/100‑1 (MQ=255)
                                                                                                gAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCAGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCt    <  2:1284254/100‑1 (MQ=255)
                                                                                                 agcagcAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTg   <  1:3411477/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
CCGCGCATTTGCCTTTCCATGCCAATAGCGCCACAGGCAATTGCCGCGAGTATAAGACGTATAATAAATTCTGCTGTCATATAAAAATGCAAAAAGGAGCAGCAAGATGGCTCAACTTGCTACTCCTTTTTACTTGCACCTCATTAATTCGGCAAGTCATTAGATGCTTGATTTTTGAAGGTATACAAATTTTTCTGC  >  NC_000913/3655823‑3656020

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: