Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,445 A→T 15.1% intergenic (‑34/+25) dtpC ← / ← cadA dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,356,4450AT15.1% 582.5 / 91.4 232intergenic (‑34/+25)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (100/96);  new base T (16/19);  total (116/116)
Fisher's exact test for biased strand distribution p-value = 5.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-02

Only 100 of 235 total aligned reads displayed.

AGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATC  >  NC_000913/4356348‑4356544
                                                                                                 |                                                                                                   
 gTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTAATTATAAAAAAAAAAGGGAAGTGGCAAg                                                                                                   <  2:2136640/99‑1 (MQ=38)
  tACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGcc                                                                                                 <  1:555315/100‑1 (MQ=255)
    cTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                               >  1:1661843/1‑100 (MQ=255)
     tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACt                                                                                              >  1:1391125/1‑100 (MQ=255)
      cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAcct                                                                                             >  1:961314/1‑98 (MQ=255)
      cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                             >  1:737897/1‑100 (MQ=255)
       ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                            <  1:2906988/100‑1 (MQ=255)
        cAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAACAAGGGAAGTGGCTTGCCACTTcc                                                                                           >  1:55181/1‑100 (MQ=255)
         cttccgatctCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                          <  2:4081098/91‑1 (MQ=255)
         aGATTTGGCGAGCCACGCTATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCAGTTcct                                                                                          <  2:2353750/100‑2 (MQ=255)
         aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                          <  2:3547889/100‑1 (MQ=255)
           aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCtt                                                                                        <  1:2500656/100‑1 (MQ=255)
             ttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttt                                                                                      <  2:55181/100‑3 (MQ=255)
             ttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGt                                                                                      <  2:1729055/100‑1 (MQ=255)
              tGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttt                                                                                     <  1:990578/100‑4 (MQ=39)
              tGGATCGCCACGATAGAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                     <  2:961308/100‑1 (MQ=255)
                 aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGa                                                                                  >  2:3176172/1‑100 (MQ=255)
                  tCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAg                                                                                 <  2:1888974/100‑1 (MQ=255)
                   cGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTcccagatcggaag                                                                                >  1:4081125/1‑90 (MQ=255)
                    gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttttgctca                                                                               <  1:486274/100‑10 (MQ=37)
                        cGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatt                                                                           >  1:2215151/1‑100 (MQ=255)
                         gATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAatta                                                                          >  1:2822278/1‑100 (MQ=255)
                          atatAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAattat                                                                         >  1:3344534/1‑100 (MQ=255)
                          atatAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAattat                                                                         >  1:3313599/1‑100 (MQ=255)
                           tagaGTATATCGTGCAAGGCTGTGAGGTTGGTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttttgc                                                                                  <  2:1014249/87‑7 (MQ=21)
                             taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttt                                                                      >  1:2594197/1‑100 (MQ=255)
                               gTATATCGCTCGCGGCTGTGAGGGTGTTTGCCTGTGTTCTCCTTATGCGCAAAAAAGGGTAGTGGCAGGCCACTTGCCTTGTACGAGCTGATTAttttt                                                                     >  2:832148/1‑99 (MQ=255)
                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                  <  2:1129344/100‑1 (MQ=255)
                                  tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCt                                                                 >  2:3054185/1‑100 (MQ=255)
                                  tatCGCGCGCGGCGGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAACAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTgat                                                                 >  2:2502974/1‑98 (MQ=255)
                                     cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttc                                                              >  2:1584219/1‑100 (MQ=255)
                                     cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttc                                                              >  1:173871/1‑100 (MQ=255)
                                      cttccgatcTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttct                                                             <  1:1109083/92‑1 (MQ=255)
                                       cgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGCGCTAATTATTTTTTGCTttctt                                                            >  2:1801284/1‑100 (MQ=255)
                                        gcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttcttc                                                           <  1:162594/100‑1 (MQ=255)
                                           cgacTGTCAAGGTGCTTTCATGTGTACTCCTTATGAGCAAAATATGGATGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTATCTGCTttcttc                                                           <  2:1658835/94‑1 (MQ=255)
                                           cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                        <  1:1901291/100‑1 (MQ=255)
                                           cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                        >  1:2020548/1‑100 (MQ=255)
                                           cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                        >  1:777340/1‑100 (MQ=255)
                                            ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                       <  1:1208697/100‑1 (MQ=255)
                                            ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                       <  1:3571362/100‑1 (MQ=255)
                                            ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                       <  1:1451647/100‑1 (MQ=255)
                                            ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                       <  2:667585/100‑1 (MQ=255)
                                             gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAATTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                      <  2:2986302/100‑1 (MQ=255)
                                             gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                      >  2:3389929/1‑100 (MQ=255)
                                             gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                      <  1:2695718/100‑1 (MQ=255)
                                              cTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                     >  2:2114980/1‑100 (MQ=255)
                                               tgtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAAt                                                    <  2:2093570/100‑1 (MQ=255)
                                               tgagagaGTGATTTCATATGTTTTTCTTATGAACAAAAAAGGGAAGTGGCAAACCACTTCCATTGTAAGAGCTAATTATTTTTTGCTTTATTATTTCaa                                                     >  2:1937717/4‑99 (MQ=255)
                                                gtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAATAAGGGAAGTTGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATa                                                   <  1:2903483/100‑1 (MQ=255)
                                                gtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCTACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATa                                                   <  2:2294525/100‑1 (MQ=255)
                                                gtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttctagatcggaag                                                   >  2:1109077/1‑90 (MQ=255)
                                                  gagagggtTTTAATGTGTTCTCCTAAAAAGAAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                  <  2:3292336/93‑1 (MQ=255)
                                                  gAGTGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGTAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                 <  2:3807531/100‑1 (MQ=255)
                                                  gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                 >  1:3615082/1‑100 (MQ=255)
                                                  gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                 >  1:1048137/1‑100 (MQ=255)
                                                  gAGCGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAATGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                 <  2:2899121/100‑1 (MQ=255)
                                                    gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                               >  2:4048772/1‑100 (MQ=255)
                                                     ggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTa                                              >  1:2089780/1‑100 (MQ=255)
                                                       ttttttCATGTGTTCTCCTTATGAGCAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                           <  1:3036639/98‑1 (MQ=255)
                                                         ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTTCTTTTTTTTTTCAATACATTAACgg                                          >  2:706909/1‑100 (MQ=255)
                                                         ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACgg                                          <  1:339957/100‑1 (MQ=255)
                                                           tgcATGTTTTCGCCTTAGGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACgg                                          <  2:1360261/96‑1 (MQ=255)
                                                                gtgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCg                                   <  1:1769051/100‑1 (MQ=255)
                                                                  gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                 <  1:2349336/100‑1 (MQ=255)
                                                                  gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                 >  1:1219282/1‑100 (MQ=255)
                                                                   ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                >  2:4104593/1‑100 (MQ=255)
                                                                     ctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCAt                              <  2:1765085/100‑1 (MQ=255)
                                                                      tcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATc                             >  2:1318897/1‑100 (MQ=255)
                                                                       ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                            >  1:1729530/1‑100 (MQ=255)
                                                                         ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAAGTATTTGTTGCTTTCTTCTGTCAATACCTTAACGGTATAGCGGCCATCAGc                          >  1:3933407/1‑100 (MQ=255)
                                                                          tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                         <  2:336631/100‑1 (MQ=255)
                                                                          cgatctcAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                         <  1:3850922/94‑1 (MQ=255)
                                                                            tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                       >  1:1137162/1‑100 (MQ=255)
                                                                            tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                       >  2:2629375/1‑100 (MQ=255)
                                                                             gctcgtacaaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGa                      >  2:391725/9‑100 (MQ=37)
                                                                              aGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGAc                     <  1:2362570/100‑1 (MQ=255)
                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                    >  1:3835981/1‑100 (MQ=255)
                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                    >  2:2701798/1‑100 (MQ=255)
                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                    >  1:662782/1‑100 (MQ=255)
                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                    >  1:3995394/1‑100 (MQ=255)
                                                                                cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGTTTTTTTTTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGAc                     >  2:1487782/1‑98 (MQ=255)
                                                                                cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCAGatcg                   >  2:3850897/1‑97 (MQ=255)
                                                                                    caaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg               >  2:3572536/2‑100 (MQ=255)
                                                                                    aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg               >  1:3182763/1‑100 (MQ=255)
                                                                                     aaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGc              <  1:2287641/100‑1 (MQ=255)
                                                                                      aGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCa             >  1:2234894/1‑100 (MQ=255)
                                                                                       gggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc            <  1:3302233/100‑1 (MQ=255)
                                                                                        ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  1:466123/1‑100 (MQ=255)
                                                                                        ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  1:3105938/1‑100 (MQ=255)
                                                                                        ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  1:3594406/1‑100 (MQ=255)
                                                                                        ggAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc           >  1:4106759/1‑100 (MQ=255)
                                                                                           aGTGGCAAGCCACTTCCCTTGTACGAGCGAATTATTTTTGGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTg        <  2:3637302/100‑1 (MQ=255)
                                                                                            gtgtCAAGCAGCTTCGCTTGTACTATCTAATTATTTTTTGCTTTCTGCTTTCAATACCTGAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTg        <  2:267355/99‑1 (MQ=255)
                                                                                            gTGGCAAGCCACTTCCCTTGTACGAGCTAATTTTTTTTTCCTGTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGa       <  2:2904900/100‑1 (MQ=255)
                                                                                            gTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGa       >  1:3552174/1‑100 (MQ=255)
                                                                                             tGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATGCCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGaa      <  2:1174514/100‑1 (MQ=255)
                                                                                             tGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGaa      >  1:720002/1‑100 (MQ=255)
                                                                                              ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat     >  1:3440234/1‑100 (MQ=255)
                                                                                                 aaGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATc  <  2:2027916/100‑1 (MQ=255)
                                                                                                 |                                                                                                   
AGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATC  >  NC_000913/4356348‑4356544

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: