Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,233 C→T 15.7% intergenic (‑123/+190) 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,2330CT15.7% 470.8 / 67.1 191intergenic (‑123/+190)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base C (77/84);  new base T (14/16);  total (91/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 198 total aligned reads displayed.

GAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAA  >  NC_000913/1357134‑1357332
                                                                                                   |                                                                                                   
gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                                     >  1:131706/1‑100 (MQ=255)
 agGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGcc                                                                                                    <  1:3136485/100‑1 (MQ=255)
 agGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGcc                                                                                                    >  1:3568060/1‑100 (MQ=255)
  ggCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCa                                                                                                   >  1:3064471/1‑100 (MQ=255)
  ggCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCa                                                                                                   >  1:1534654/1‑100 (MQ=255)
  ggCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCa                                                                                                   >  2:4277548/1‑100 (MQ=255)
   gCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCAt                                                                                                  <  1:1115223/100‑1 (MQ=255)
    cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                                 >  2:3516329/1‑100 (MQ=255)
    cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                                 >  2:3319384/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)
            cTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAc                                                                                         <  2:3147937/100‑1 (MQ=37)
             tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACa                                                                                        >  2:2713203/1‑100 (MQ=255)
             tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACa                                                                                        >  2:639820/1‑100 (MQ=255)
             tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACa                                                                                        <  1:1884267/100‑1 (MQ=255)
             tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                        <  2:3059601/100‑2 (MQ=37)
             tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                        <  2:1089419/100‑2 (MQ=37)
             tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                        <  2:723813/100‑2 (MQ=37)
              aaTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACAt                                                                                       >  1:3876380/1‑100 (MQ=255)
              aaTATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGCTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACAt                                                                                       >  1:2120569/1‑100 (MQ=255)
               atatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcgac                                                                                      <  2:2124539/100‑4 (MQ=25)
                tatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTTATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACa                                                                                        <  1:271838/97‑1 (MQ=255)
                tatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAGCATTc                                                                                     <  2:809306/100‑1 (MQ=255)
                tatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTc                                                                                     <  2:971557/100‑1 (MQ=255)
                 atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                                    >  1:718007/1‑100 (MQ=255)
                 atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                                    <  2:4131646/100‑1 (MQ=255)
                  tttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGt                                                                                   <  1:985064/100‑1 (MQ=255)
                      aTGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAccc                                                                               >  2:1188534/1‑100 (MQ=255)
                       ttaCATTTGCCATTACCGTTGGGGTTCAGGGGTCGTAATGAGCGCCGTGTTTTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAcc                                                                                <  1:1707796/96‑1 (MQ=255)
                       tGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCg                                                                              >  2:1422486/1‑100 (MQ=255)
                              tGCCGTTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGTTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtc                                                                       <  1:3350021/100‑1 (MQ=255)
                                ccATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtctc                                                                     >  2:103376/1‑100 (MQ=255)
                                ccATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTtctc                                                                     <  1:901258/100‑1 (MQ=255)
                                 cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                                    >  1:3692087/1‑100 (MQ=255)
                                    tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                                 <  2:1870030/100‑1 (MQ=255)
                                       cGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataact                                                              <  2:2659770/100‑1 (MQ=255)
                                        gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGATGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                             <  2:2144472/100‑1 (MQ=255)
                                        gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                             <  1:1100765/100‑1 (MQ=255)
                                          ttGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                             <  1:3044068/98‑1 (MQ=255)
                                           tgtgTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataactataa                                                          >  1:2377324/1‑100 (MQ=255)
                                              nttCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTc                                                       >  2:1013662/2‑100 (MQ=255)
                                              gtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCTTACCCGTGATTTCTCTATAACTATAACTc                                                       >  1:449956/1‑100 (MQ=255)
                                              gtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTataactataact                                                        >  1:3678311/1‑99 (MQ=255)
                                               ttCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTca                                                      >  2:1227975/1‑100 (MQ=255)
                                                 cAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcaca                                                    >  2:1790621/1‑100 (MQ=255)
                                                  aGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAg                                                   >  2:720086/1‑100 (MQ=255)
                                                    gggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                                 >  2:4002108/1‑100 (MQ=255)
                                                     ggtggtAATGAGGGCCGTGTAGTTGGTCATTGGACAGGGAACGTTGCAATAAGATAAACATTCGTAACCGTGAATTCGCTATAATTCTAGCTCACATaa                                                 <  1:2303632/99‑1 (MQ=255)
                                                     ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaac                                                <  1:129831/100‑1 (MQ=255)
                                                      gTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaaca                                               >  2:3360605/1‑100 (MQ=255)
                                                        cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaac                                             >  2:2017207/1‑100 (MQ=255)
                                                        cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaac                                             <  1:780885/100‑1 (MQ=255)
                                                         gTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACGCGTGATTTCTCTATAACTATAACTCACAGAACAACt                                            >  1:3269358/1‑100 (MQ=255)
                                                          tAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACtt                                           >  2:3530210/1‑100 (MQ=255)
                                                          tAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACtt                                           >  1:1946075/1‑100 (MQ=255)
                                                           aaTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                          >  1:2603963/1‑100 (MQ=255)
                                                           aaTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTa                                          <  1:2699929/100‑1 (MQ=255)
                                                            aTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                         >  2:1449829/1‑100 (MQ=255)
                                                            aTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                         >  2:3689368/1‑100 (MQ=255)
                                                             tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGc                                        >  1:4031053/1‑100 (MQ=255)
                                                              gAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCg                                       <  1:2473022/100‑1 (MQ=255)
                                                                gggCCGTGTGGTTGGTCGTTGAACAGGTGACATCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgag                                     >  1:221585/1‑100 (MQ=255)
                                                                  gCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCATAGAACAACTTAGCgagga                                   >  2:2215689/1‑100 (MQ=255)
                                                                   ccGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgaggag                                  <  1:408974/100‑1 (MQ=255)
                                                                     gtgtGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCa                                >  1:552458/1‑100 (MQ=255)
                                                                      tgtgGTTGGTCGTTGAACAGGTGACGTCGCCATCCGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCag                               <  2:259385/100‑1 (MQ=255)
                                                                        ggtacgaatgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGCGCagag                             >  1:1343668/10‑100 (MQ=21)
                                                                         ggttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                            >  2:292098/1‑100 (MQ=255)
                                                                         ggttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                            >  2:2994376/1‑100 (MQ=255)
                                                                           ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATCACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                           >  1:3572433/1‑99 (MQ=255)
                                                                           ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                          >  2:154265/1‑100 (MQ=255)
                                                                           ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                          <  2:165905/100‑1 (MQ=255)
                                                                            tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGt                         <  2:4309001/100‑1 (MQ=255)
                                                                            cgaatgTTGAACAGATGGCGACGTCACCTGTTCAACGTTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGt                         >  2:3721568/6‑100 (MQ=21)
                                                                              gtCCTTGAACAGTTGACTTCTCCTTCTTTTTAACATTCTTACCTGTTATTTCTCTATAACTATAACTCACCGCTCTTTTTATCGACGAGCAGAGccg                          <  1:978322/97‑3 (MQ=255)
                                                                               atgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACg                      >  2:3457402/3‑100 (MQ=37)
                                                                                 gTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt                    >  1:3502324/1‑100 (MQ=255)
                                                                                 gTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGGACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc                       <  1:2801003/97‑1 (MQ=255)
                                                                                 gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt                    >  1:1417263/1‑100 (MQ=37)
                                                                                   tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg                  <  1:2206042/100‑1 (MQ=255)
                                                                                   tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg                  <  2:4031082/100‑1 (MQ=255)
                                                                                    gAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt                 >  1:1831889/1‑100 (MQ=255)
                                                                                    gAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt                 >  2:3576923/1‑100 (MQ=255)
                                                                                     aaCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTg                <  1:1892005/100‑1 (MQ=255)
                                                                                     aaCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTg                <  2:3921610/100‑1 (MQ=255)
                                                                                      aCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTga               <  2:2768645/100‑1 (MQ=255)
                                                                                      aCAGGCGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTga               <  2:2088537/100‑1 (MQ=255)
                                                                                         ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa            >  1:2074840/1‑100 (MQ=255)
                                                                                         ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa            >  1:638894/1‑100 (MQ=255)
                                                                                            gACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAAcc         >  2:508041/1‑100 (MQ=255)
                                                                                              cGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGCACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCa       >  1:612333/1‑100 (MQ=255)
                                                                                               gTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACTACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCAt      <  1:2399696/100‑1 (MQ=255)
                                                                                               gTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCAt      >  1:121131/1‑100 (MQ=255)
                                                                                                tCGCCATCTGTTCAACATTCGTACCCTTGATTTCTCTATAACTATAACTCACAGAACAACATAGCGAGGAGCAGATCGGTACGATCGTGAGAAACCCaat     <  1:2916356/100‑3 (MQ=255)
                                                                                                tCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCAtt     <  2:2459624/100‑1 (MQ=255)
                                                                                                 cGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTa    <  1:3013519/100‑1 (MQ=255)
                                                                                                   ccATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaaa  >  1:3476113/1‑100 (MQ=255)
                                                                                                   ccATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTaaa  >  2:4169312/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCATTAAA  >  NC_000913/1357134‑1357332

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: