Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,227 A→G 16.3% intergenic (‑117/+196) 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,2270AG16.3% 298.5 / 46.3 133intergenic (‑117/+196)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base A (64/47);  new base G (12/10);  total (76/57)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.45e-01

Only 100 of 136 total aligned reads displayed.

TGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACC  >  NC_000913/1357128‑1357325
                                                                                                   |                                                                                                  
tGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGa                                                                                                    >  1:3620450/1‑100 (MQ=255)
   aaaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTTAGGTGTCGTAATGAGTGCCGTGTGGTTGGTCGTTGAAGAGGTGagg                                                                                                  >  2:2326396/1‑97 (MQ=255)
   aaaGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGAGTCGTAATGAGGGCCGTGTGTTTGGTCGTTGAACAGGTGACGt                                                                                                 >  2:2222400/1‑100 (MQ=255)
     agagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                               <  2:2331462/100‑1 (MQ=255)
      gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                              >  1:869751/1‑100 (MQ=255)
      gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                              <  1:189992/100‑1 (MQ=255)
       agGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGcc                                                                                             >  1:295788/1‑100 (MQ=255)
          cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                          >  2:925999/1‑100 (MQ=255)
          cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                          >  1:1705870/1‑100 (MQ=255)
            aTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTg                                                                                        <  2:1878751/100‑1 (MQ=255)
            aTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTg                                                                                        <  1:2789493/100‑1 (MQ=255)
               ttacgtatatTTAACGACATTTGCCATTACCCTTTGTGTTCAGGGGTCGTAATCAGGGCCATGTCCTTGGTCGTGGAACAGGTGACGACGCCATCTGTTc                                                                                     <  2:1708098/94‑1 (MQ=255)
               ttACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTc                                                                                     >  1:675280/1‑100 (MQ=255)
                tACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTcn                                                                                    <  1:1508216/100‑2 (MQ=255)
                tACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCa                                                                                    >  1:3389799/1‑100 (MQ=255)
                tACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCa                                                                                    >  1:2992759/1‑100 (MQ=255)
                  cTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAc                                                                                  <  1:2162360/100‑1 (MQ=255)
                  cTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAc                                                                                  <  1:822287/100‑1 (MQ=37)
                  cTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAc                                                                                  <  1:822154/100‑1 (MQ=37)
                   tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCaa                                                                                   <  2:500744/98‑1 (MQ=37)
                   tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                 <  2:628536/100‑2 (MQ=37)
                   tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                 <  1:2008034/100‑2 (MQ=37)
                     aacggtAATGGCAAATGTCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTAGTCGTCGAACAGGTGACGTCGCCATCTGCTCAACAt                                                                                >  2:150230/6‑99 (MQ=255)
                      tatTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcgacc                                                                              <  1:1509521/100‑5 (MQ=25)
                       atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGTGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAAGAGGTGACGTCGCCATCTGTTCAACATTCg                                                                             >  2:1596526/1‑100 (MQ=255)
                       atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                             <  1:2762093/100‑1 (MQ=255)
                       atTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCg                                                                             <  2:1517442/100‑1 (MQ=255)
                        tttAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGt                                                                            >  2:422837/1‑100 (MQ=255)
                          tAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAc                                                                          >  2:2152299/1‑100 (MQ=255)
                           aaTGACATTTGCCATTACCGTTTGTGTTCAGTGGTCATAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAcc                                                                         <  1:707541/100‑1 (MQ=255)
                           aaTGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAcc                                                                         >  2:1164278/1‑100 (MQ=255)
                               aCATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTg                                                                     >  1:1979073/1‑100 (MQ=255)
                                cATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAccc                                                                        <  2:1362496/96‑1 (MQ=255)
                                cATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGa                                                                    >  1:1634756/1‑100 (MQ=255)
                                cATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGa                                                                    >  2:3300204/1‑100 (MQ=255)
                                 aTTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAt                                                                   >  2:2472852/1‑100 (MQ=255)
                                       cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGCGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                             <  1:547425/100‑1 (MQ=255)
                                         ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctat                                                           <  2:207776/100‑1 (MQ=255)
                                          tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                          <  1:1422767/100‑1 (MQ=255)
                                          tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                          >  1:1787599/1‑100 (MQ=255)
                                          tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                          >  2:816706/1‑100 (MQ=255)
                                             cGTTTGTGTTCAGGTGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataact                                                       >  2:831849/1‑100 (MQ=255)
                                             cGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataact                                                       >  1:3212329/1‑100 (MQ=255)
                                             cGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataact                                                       >  2:352366/1‑100 (MQ=255)
                                              gTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGGCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataacta                                                      <  2:621492/100‑1 (MQ=255)
                                                ttGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctataac                                                        <  2:1418668/96‑1 (MQ=255)
                                                    gtTCAGGGGTCGTAATGAGGGCCGTGTGGATGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAGATTCGTACACGTGATTTCTCtataactataacgc                                                >  2:1363887/1‑98 (MQ=255)
                                                     ttCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTca                                               >  1:1239637/1‑100 (MQ=255)
                                                      tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcac                                              >  1:1764543/1‑100 (MQ=255)
                                                        aGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAg                                            <  1:422839/100‑1 (MQ=255)
                                                           ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaac                                         >  1:3322288/1‑100 (MQ=255)
                                                           ggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaac                                         <  1:1917777/100‑1 (MQ=255)
                                                            gTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaaca                                        >  1:2923654/1‑100 (MQ=255)
                                                             tCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaa                                       <  1:715276/100‑1 (MQ=255)
                                                              cGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaacaac                                      <  2:1175186/100‑1 (MQ=255)
                                                               gTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGAAGTCGACATCTGTTCAACATTCGTACCCGTGATTTCTCTAGAACTATAACTCACAGaacaac                                      >  1:1762377/1‑99 (MQ=255)
                                                                  aTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                  >  2:1143646/1‑100 (MQ=255)
                                                                   tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGc                                 >  1:2926246/1‑100 (MQ=255)
                                                                    gAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCg                                >  1:1456740/1‑100 (MQ=255)
                                                                      gggCCGTGTGGTTGGTCGTTGAACAGGTGGCGTGGCCATCTGTTCAACATTCGTACCTGTGATTTCTCTATAACTATTACTCACAGAACAACTTAGCgag                              >  2:1073204/1‑100 (MQ=255)
                                                                      gggCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgag                              <  2:3433746/100‑1 (MQ=255)
                                                                      gggCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgag                              >  2:2987429/1‑100 (MQ=255)
                                                                        gCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagga                            <  2:1456732/100‑1 (MQ=255)
                                                                        gCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagga                            >  1:1474548/1‑100 (MQ=255)
                                                                          cGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGc                          >  2:583722/1‑100 (MQ=255)
                                                                           gtgtGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCa                         <  1:2922797/100‑1 (MQ=255)
                                                                           gggtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCa                         <  1:909335/98‑1 (MQ=255)
                                                                            tgtgGTTGGTCTTTGAACAGGTTACGTCGCCATCTGTTCAACATGCGTGCCAGTGATTGCGCTATAACGATAACTCACAGAACAACTTATCtagg                             <  2:2904399/95‑1 (MQ=255)
                                                                            tgtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCag                        <  1:3142318/100‑1 (MQ=255)
                                                                                tacgaatgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                    >  1:2712837/8‑100 (MQ=21)
                                                                                gttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                    >  1:1165197/1‑100 (MQ=255)
                                                                                 ttggtCGTTGAACAGGTGACGTCGCCCTCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                   <  1:1627111/100‑1 (MQ=255)
                                                                                  tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                   <  1:2343113/99‑1 (MQ=255)
                                                                                    gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc                >  1:56895/1‑100 (MQ=255)
                                                                                    gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc                <  2:2919757/100‑1 (MQ=255)
                                                                                    gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc                <  2:3212304/100‑1 (MQ=255)
                                                                                     tCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACg               <  2:2949041/100‑1 (MQ=255)
                                                                                     atgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACg               >  1:3513250/3‑100 (MQ=37)
                                                                                      cGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGa              >  2:2597628/1‑100 (MQ=255)
                                                                                      cGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGAACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGa              >  2:1940210/1‑100 (MQ=255)
                                                                                       tccgatctggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt             <  1:1038888/92‑1 (MQ=255)
                                                                                       gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt             >  2:754167/1‑100 (MQ=37)
                                                                                       gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt             >  1:2281541/1‑100 (MQ=37)
                                                                                       gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCACTATAACTATAACTCACAGAACAACTTAGCGAAGAGCAGAGCGGTAc                >  2:932162/1‑97 (MQ=25)
                                                                                        ttGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGATCAGAGCGGTAGGacc            >  1:3585660/1‑98 (MQ=37)
                                                                                        ttGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATc            >  2:3101756/1‑100 (MQ=37)
                                                                                        ttGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATc            >  2:1234085/1‑100 (MQ=37)
                                                                                         tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg           <  1:3247498/100‑1 (MQ=255)
                                                                                         tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg           <  1:2152313/100‑1 (MQ=255)
                                                                                          gAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt          >  2:1171937/1‑100 (MQ=255)
                                                                                           aaCAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTg         <  1:3160434/100‑1 (MQ=255)
                                                                                             cAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTgag       <  1:2347299/100‑1 (MQ=255)
                                                                                               ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAtagatcgga     >  2:1038880/1‑92 (MQ=255)
                                                                                               ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa     >  2:3615750/1‑100 (MQ=255)
                                                                                               ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa     >  2:2748943/1‑100 (MQ=255)
                                                                                               ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa     >  2:2615890/1‑100 (MQ=255)
                                                                                               ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa     >  1:1980992/1‑100 (MQ=255)
                                                                                               ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAAATTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa     >  1:892686/1‑100 (MQ=255)
                                                                                                gTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaaa    >  1:1607626/1‑100 (MQ=255)
                                                                                                 tGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATATCTATAGCTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAAc   <  2:3595122/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
TGTAAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACC  >  NC_000913/1357128‑1357325

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