Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,224 G→A 15.3% intergenic (‑114/+199) sapA ← / ← ymjA antimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,357,2240GA15.3% 493.4 / 64.0 195intergenic (‑114/+199)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base G (83/80);  new base A (14/16);  total (97/98)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.74e-01

Only 100 of 201 total aligned reads displayed.

TTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAA  >  NC_000913/1357125‑1357322
                                                                                                   |                                                                                                  
ttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAgg                                                                                                    <  2:2872647/100‑1 (MQ=255)
ttATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAgg                                                                                                    <  2:4200935/100‑1 (MQ=255)
 tATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGt                                                                                                   >  2:2064178/1‑100 (MQ=255)
  aTGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTg                                                                                                  >  1:3851194/1‑100 (MQ=255)
   tGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGa                                                                                                 <  2:3425713/100‑1 (MQ=255)
   tGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGa                                                                                                 >  1:2368073/1‑100 (MQ=255)
   ggTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTg                                                                                                  <  1:3477270/98‑1 (MQ=255)
    gTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGAc                                                                                                >  1:339876/1‑100 (MQ=255)
      aaaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGt                                                                                              >  2:2860762/1‑100 (MQ=255)
      aaaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGt                                                                                              >  1:3929349/1‑100 (MQ=255)
        agagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                            >  1:1396117/1‑100 (MQ=255)
         gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                           >  1:131706/1‑100 (MQ=255)
          agGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGcc                                                                                          >  1:3568060/1‑100 (MQ=255)
           ggCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCa                                                                                         >  2:4277548/1‑100 (MQ=255)
             cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGGGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                       <  1:2120535/100‑1 (MQ=255)
             cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                       <  2:1460524/100‑1 (MQ=255)
             cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                       >  2:3516329/1‑100 (MQ=255)
             cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTAATGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                       >  1:2731012/1‑100 (MQ=255)
               agcTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTGCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                       <  2:3557931/96‑1 (MQ=255)
                   tACTAATATTTAATGACGTTTGACATTACCGTTTGTGTTCAGGGGTGGTAATGAGGGGCGTGTGGTGGGTCGTTGAACAGGTGACGTCGCCATCTGtt                                                                                   <  1:241574/98‑1 (MQ=255)
                    agtagtatTTAATGACATTGTGTATTTCGGTGTGTGATGAGGGGGCGTAATGAGTGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCaa                                                                                <  1:3447506/95‑1 (MQ=9)
                     cTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAc                                                                               <  2:3917978/100‑1 (MQ=37)
                      tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACa                                                                              >  1:2198324/1‑100 (MQ=255)
                      tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACa                                                                              >  2:639820/1‑100 (MQ=255)
                      tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                              <  1:2990933/100‑2 (MQ=37)
                      tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                              <  2:723813/100‑2 (MQ=37)
                      tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCACGGGTCGTAATGAGGGCCGTTTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                              <  1:2394287/100‑2 (MQ=25)
                        atatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcgac                                                                            <  1:1937784/100‑4 (MQ=25)
                         tatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAGCATTc                                                                           <  2:809306/100‑1 (MQ=255)
                          atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGTGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGTTGACGTCGCCATCTGTTCAACATTCg                                                                          <  1:227209/100‑1 (MQ=255)
                          atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                          >  1:1959396/1‑100 (MQ=255)
                          atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                          >  1:718007/1‑100 (MQ=255)
                          atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                          <  2:4131646/100‑1 (MQ=255)
                            ttAATGACATTTGACATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTa                                                                        >  2:2992425/1‑100 (MQ=255)
                               aTGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAccc                                                                     >  2:2206050/1‑100 (MQ=255)
                                 gACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGt                                                                   >  2:2750668/1‑100 (MQ=255)
                                         ccATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtctc                                                           <  1:3018573/100‑1 (MQ=255)
                                         ccATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtctc                                                           <  1:901258/100‑1 (MQ=255)
                                          cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                          >  1:3692087/1‑100 (MQ=255)
                                             tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                       <  2:1870030/100‑1 (MQ=255)
                                               ccGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGTTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                     <  1:1616406/100‑1 (MQ=255)
                                               ccGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                     <  1:591839/100‑1 (MQ=255)
                                               ccGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGGTCAACATTCGTACCCGTGATTTCTctataac                                                     >  1:1820517/1‑100 (MQ=255)
                                                cGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataact                                                    <  2:2659770/100‑1 (MQ=255)
                                                 gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                   <  1:1100765/100‑1 (MQ=255)
                                                 gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                   >  2:3013536/1‑100 (MQ=255)
                                                  tttGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataactat                                                  <  1:1891841/100‑1 (MQ=255)
                                                      tgtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTataactataact                                              >  2:122613/1‑100 (MQ=255)
                                                       gtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTataactataact                                              >  1:3678311/1‑99 (MQ=255)
                                                         tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCAGTGATTTCTCTATAACTATAACTcac                                           <  2:1396122/100‑1 (MQ=255)
                                                           aGGGGTCGTAATGAGGGCCGTGTGGTTTGTCGTTGAACAGGTGACTTCGCTATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcaca                                          <  1:1848302/99‑1 (MQ=255)
                                                            ggggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGa                                        >  2:3227055/1‑100 (MQ=255)
                                                             gggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                       >  2:4002108/1‑100 (MQ=255)
                                                              ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaac                                      <  1:129831/100‑1 (MQ=255)
                                                                 cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaac                                   <  1:780885/100‑1 (MQ=255)
                                                                 cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaac                                   >  1:2206442/1‑100 (MQ=255)
                                                                 cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaac                                   >  2:2017207/1‑100 (MQ=255)
                                                                   tAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACtt                                 >  2:3530210/1‑100 (MQ=255)
                                                                     aTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                               >  2:3689368/1‑100 (MQ=255)
                                                                     aTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                               <  2:2225565/100‑1 (MQ=255)
                                                                      tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGc                              >  1:4031053/1‑100 (MQ=255)
                                                                       gAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCg                             <  1:2473022/100‑1 (MQ=255)
                                                                        aGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCga                            >  2:1650801/1‑100 (MQ=255)
                                                                         gggCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgag                           >  2:2127432/1‑100 (MQ=255)
                                                                          ggtcGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACTCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagg                          <  1:3058374/97‑1 (MQ=255)
                                                                            ccGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgaggag                        <  1:3268309/100‑1 (MQ=255)
                                                                            ccGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgaggag                        <  1:408974/100‑1 (MQ=255)
                                                                            ccGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGGTCAACATTCGTACCCGTGATTTCGCTATAACTAAAACTCACAGAACAACTTAGCgaggag                        <  1:550532/100‑1 (MQ=255)
                                                                             cGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGc                       >  2:3357500/1‑100 (MQ=255)
                                                                               tgtgGTTGGTCGTTGAACAGGTGACGTCGCCATCCGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCag                     <  2:259385/100‑1 (MQ=255)
                                                                                 ggtacgaatgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGCGCagag                   >  1:1343668/10‑100 (MQ=21)
                                                                                  ggttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                  <  2:2120577/100‑1 (MQ=255)
                                                                                   gttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                 <  1:262425/100‑1 (MQ=255)
                                                                                    ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATCACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                 >  1:3572433/1‑99 (MQ=255)
                                                                                     tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGt               <  2:4309001/100‑1 (MQ=255)
                                                                                     cgaatgTTGAACAGATGGCGACGTCACCTGTTCAACGTTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGt               >  2:3721568/6‑100 (MQ=21)
                                                                                      ggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTa              <  1:3071259/100‑1 (MQ=255)
                                                                                       gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc             <  2:3064488/100‑1 (MQ=255)
                                                                                       gtCCTTGAACAGTTGACTTCTCCTTCTTTTTAACATTCTTACCTGTTATTTCTCTATAACTATAACTCACCGCTCTTTTTATCGACGAGCAGAGccg                <  1:978322/97‑3 (MQ=255)
                                                                                        atgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTa              >  2:712196/3‑98 (MQ=37)
                                                                                        atgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACg            >  2:3457402/3‑100 (MQ=37)
                                                                                         tgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGa           >  2:11823/2‑100 (MQ=37)
                                                                                          gTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt          >  1:3502324/1‑100 (MQ=255)
                                                                                          gTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGGACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc             <  1:2801003/97‑1 (MQ=255)
                                                                                          gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt          >  2:281623/1‑100 (MQ=37)
                                                                                           ttGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATc         <  2:1029614/100‑1 (MQ=255)
                                                                                           ttGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATc         >  2:1061539/1‑100 (MQ=37)
                                                                                            tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg        <  2:1843795/100‑1 (MQ=255)
                                                                                            tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg        <  1:1188529/100‑1 (MQ=255)
                                                                                            tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg        <  2:4031082/100‑1 (MQ=255)
                                                                                             gAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt       >  1:3342357/1‑100 (MQ=255)
                                                                                             gAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt       >  2:3576923/1‑100 (MQ=255)
                                                                                              aaCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTg      <  2:3921610/100‑1 (MQ=255)
                                                                                              aaCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTg      >  1:4341081/1‑100 (MQ=255)
                                                                                              aaCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTg      <  2:1746858/100‑1 (MQ=255)
                                                                                              aaCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGAAGCAGAGCGGTACGATCGTg      <  1:1465392/100‑1 (MQ=255)
                                                                                               aCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTga     <  1:2713188/100‑1 (MQ=255)
                                                                                                cAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTgag    <  2:1425125/100‑1 (MQ=255)
                                                                                                  ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa  >  1:638894/1‑100 (MQ=255)
                                                                                                  ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa  <  2:1009412/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TTATGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAA  >  NC_000913/1357125‑1357322

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: