Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,357,230 T→A 16.3% intergenic (‑120/+193) 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,2300TA16.3% 327.1 / 50.4 134intergenic (‑120/+193)sapA/ymjAantimicrobial peptide transport ABC transporter periplasmic binding protein/DUF2543 family protein
Reads supporting (aligned to +/- strand):  ref base T (65/46);  new base A (11/11);  total (76/58)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01

Only 100 of 135 total aligned reads displayed.

AAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCA  >  NC_000913/1357131‑1357327
                                                                                                   |                                                                                                 
  agagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCg                                                                                                 <  2:2331462/100‑1 (MQ=255)
   gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                                >  1:869751/1‑100 (MQ=255)
   gagGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGc                                                                                                <  1:189992/100‑1 (MQ=255)
       cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                            >  1:1705870/1‑100 (MQ=255)
       cTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATc                                                                                            >  2:925999/1‑100 (MQ=255)
         aTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTg                                                                                          <  1:2789493/100‑1 (MQ=255)
         aTCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTg                                                                                          <  2:1878751/100‑1 (MQ=255)
            ttacgtatatTTAACGACATTTGCCATTACCCTTTGTGTTCAGGGGTCGTAATCAGGGCCATGTCCTTGGTCGTGGAACAGGTGACGACGCCATCTGTTc                                                                                       <  2:1708098/94‑1 (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)
                tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCaa                                                                                     <  2:500744/98‑1 (MQ=37)
                tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                   <  2:628536/100‑2 (MQ=37)
                tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                   <  1:2008034/100‑2 (MQ=37)
                tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                   <  1:3624783/100‑2 (MQ=37)
                tAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGATGGCGACGTCACCTGTTCAAcg                                                                                   <  1:578544/100‑2 (MQ=37)
                  atctttAATGACATTTGCCATTAGCGTTTGTGTTGAGGGGTCGTAATGAGGGCCGTGTGGTTCGTCGTTGAACAGATGGGGACGTCACCTGTTTAAcg                                                                                   <  2:3485703/95‑2 (MQ=255)
                  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)
                     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:3157907/100‑1 (MQ=255)
                             cATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTAccc                                                                          <  2:1362496/96‑1 (MQ=255)
                             cATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGa                                                                      >  2:3300204/1‑100 (MQ=255)
                             cATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGa                                                                      >  1:1634756/1‑100 (MQ=255)
                              aTTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAt                                                                     >  2:2472852/1‑100 (MQ=255)
                              aTTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGAt                                                                     <  1:1781518/100‑1 (MQ=255)
                                    cATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGCGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTctct                                                               <  1:547425/100‑1 (MQ=255)
                                      ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctat                                                             <  2:207776/100‑1 (MQ=255)
                                      ttACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctat                                                             >  1:2949064/1‑100 (MQ=255)
                                       tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                            >  2:816706/1‑100 (MQ=255)
                                       tACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTctata                                                            <  1:1422767/100‑1 (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)
                                               gtgtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCtataactataac                                                    >  2:3205383/1‑100 (MQ=255)
                                                tgtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTataactataact                                                   >  2:1198142/1‑100 (MQ=255)
                                                 gtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTc                                                  >  2:1470702/1‑100 (MQ=255)
                                                 gtTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTc                                                  >  1:2921735/1‑100 (MQ=255)
                                                 gtTCAGGGGTCGTAATGAGGGCCGTGTGGATGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAAGATTCGTACACGTGATTTCTCtataactataacgc                                                  >  2:1363887/1‑98 (MQ=255)
                                                 gtTCAGGGGTCGTAATGAGGGCCGAGTGGTGGGTCGTTGAACAGGTGACGTCGCAATCTGTTCAAACTTCGTACCCGTGCTTTATCtagaactataac                                                    >  2:2945802/1‑98 (MQ=255)
                                                  ttCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTca                                                 >  1:1239637/1‑100 (MQ=255)
                                                   tCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcac                                                >  1:1764543/1‑100 (MQ=255)
                                                    cAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTcaca                                               <  2:2349107/100‑1 (MQ=255)
                                                       gggTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGaa                                            >  2:140625/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                                        <  1:721269/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)
                                                               aTGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAg                                    >  2:1144122/1‑100 (MQ=255)
                                                                tGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGc                                   >  1:2926246/1‑100 (MQ=255)
                                                                 gAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCg                                  >  1:1456740/1‑100 (MQ=255)
                                                                   gggCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgag                                >  2:2987429/1‑100 (MQ=255)
                                                                     gCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagga                              >  1:1474548/1‑100 (MQ=255)
                                                                     gCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCgagga                              <  2:1456732/100‑1 (MQ=255)
                                                                       cGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGc                            >  2:583722/1‑100 (MQ=255)
                                                                        gggtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCa                           <  1:909335/98‑1 (MQ=255)
                                                                         tgtgGTTGGTCTTTGAACAGGTTACGTCGCCATCTGTTCAACATGCGTGCCAGTGATTGCGCTATAACGATAACTCACAGAACAACTTATCtagg                               <  2:2904399/95‑1 (MQ=255)
                                                                         tgtgGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCag                          <  1:3142318/100‑1 (MQ=255)
                                                                            ggttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGc                       >  1:943077/1‑100 (MQ=255)
                                                                             gttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCg                      >  1:1165197/1‑100 (MQ=255)
                                                                              ttggtCGTTGAACAGGTGACGTCGCCCTCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                     <  1:1627111/100‑1 (MQ=255)
                                                                              ttggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                     <  1:2667569/100‑1 (MQ=255)
                                                                               tggtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCgg                     <  1:2343113/99‑1 (MQ=255)
                                                                                 gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc                  <  2:3212304/100‑1 (MQ=255)
                                                                                 gtCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTAc                  >  1:56895/1‑100 (MQ=255)
                                                                                  tCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACg                 <  2:2949041/100‑1 (MQ=255)
                                                                                  atgTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACg                 >  1:3513250/3‑100 (MQ=37)
                                                                                   cGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGa                >  2:2597628/1‑100 (MQ=255)
                                                                                    tccgatctggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt               <  1:1038888/92‑1 (MQ=255)
                                                                                    gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt               >  1:2281541/1‑100 (MQ=37)
                                                                                    gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAt               >  2:754167/1‑100 (MQ=37)
                                                                                    gTTGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCACTATAACTATAACTCACAGAACAACTTAGCGAAGAGCAGAGCGGTAc                  >  2:932162/1‑97 (MQ=25)
                                                                                     ttGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGATCAGAGCGGTAGGacc              >  1:3585660/1‑98 (MQ=37)
                                                                                     ttGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATc              >  1:851822/1‑100 (MQ=37)
                                                                                      tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg             <  1:2152313/100‑1 (MQ=255)
                                                                                      tGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg             <  1:3247498/100‑1 (MQ=255)
                                                                                      tGAACAGATGGCGACGTCACCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCg             >  1:3275675/1‑100 (MQ=38)
                                                                                       gAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGt            >  2:1998905/1‑100 (MQ=255)
                                                                                          cAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTgag         <  1:2347299/100‑1 (MQ=255)
                                                                                            ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGAtagatcgga       >  2:1038880/1‑92 (MQ=255)
                                                                                            ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa       >  2:2615890/1‑100 (MQ=255)
                                                                                            ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa       >  1:1980992/1‑100 (MQ=255)
                                                                                            ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa       >  2:3615750/1‑100 (MQ=255)
                                                                                            ggTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaa       >  2:2748943/1‑100 (MQ=255)
                                                                                             gTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGaaa      >  1:1607626/1‑100 (MQ=255)
                                                                                              tGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATATCTATAGCTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAAc     <  2:3595122/100‑1 (MQ=255)
                                                                                              tGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAAc     >  1:1934647/1‑100 (MQ=255)
                                                                                               gACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAAcc    >  2:2296778/1‑100 (MQ=255)
                                                                                                 cGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCa  >  1:3422745/1‑100 (MQ=255)
                                                                                                   |                                                                                                 
AAAGAGGCTATCTTACTAATATTTAATGACATTTGCCATTACCGTTTGTGTTCAGGGGTCGTAATGAGGGCCGTGTGGTTGGTCGTTGAACAGGTGACGTCGCCATCTGTTCAACATTCGTACCCGTGATTTCTCTATAACTATAACTCACAGAACAACTTAGCGAGGAGCAGAGCGGTACGATCGTGAGAAACCCA  >  NC_000913/1357131‑1357327

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: