Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,356,446 A→T 15.3% intergenic (‑35/+24) 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,4460AT15.3% 591.4 / 91.0 232intergenic (‑35/+24)dtpC/cadAdipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand):  ref base A (99/96);  new base T (16/19);  total (116/116)
Fisher's exact test for biased strand distribution p-value = 7.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.84e-01

Only 100 of 236 total aligned reads displayed.

AGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCG  >  NC_000913/4356348‑4356545
                                                                                                  |                                                                                                   
 gTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTAATTATAAAAAAAAAAGGGAAGTGGCAAg                                                                                                    <  2:2136640/99‑1 (MQ=38)
    cTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                >  1:1661843/1‑100 (MQ=255)
     tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACt                                                                                               >  1:1391125/1‑100 (MQ=255)
      cccAGATTTGGATGGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                              >  1:2802403/1‑100 (MQ=255)
      cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                              >  1:737897/1‑100 (MQ=255)
      cccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                              >  1:3256307/1‑100 (MQ=255)
       ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTc                                                                                             <  1:2906988/100‑1 (MQ=255)
         cttccgatctCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTccc                                                                                           <  2:4081098/91‑1 (MQ=255)
         aGATTTGGCGAGCCACGCTATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCAGTTcct                                                                                           <  2:2353750/100‑2 (MQ=255)
         aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                           <  1:3077532/100‑1 (MQ=255)
          gATTTGGATCGCCACGATATAGTATATCGCGCGTGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTcc                                                                                            <  1:2592016/98‑1 (MQ=255)
          gATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCt                                                                                          <  1:672449/100‑1 (MQ=255)
           aTTTGGATCGCCACGATATAGTATATCGCGCGCTGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                         <  1:437124/100‑1 (MQ=39)
           aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                         <  1:3583237/100‑1 (MQ=255)
           aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                         <  1:2184360/100‑1 (MQ=255)
           aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                         <  2:2286435/100‑1 (MQ=255)
           aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCtt                                                                                         <  1:3305553/100‑1 (MQ=255)
           aTTTGGATCGCAACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtt                                                                                         <  1:3404218/100‑1 (MQ=39)
             ttGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttt                                                                                       <  2:55181/100‑3 (MQ=255)
              tGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttt                                                                                      <  1:990578/100‑4 (MQ=39)
              tGGATCGCCACGATAGAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                      <  2:961308/100‑1 (MQ=255)
                gATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACg                                                                                    <  2:1406396/100‑1 (MQ=255)
                 aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGa                                                                                   >  2:3176172/1‑100 (MQ=255)
                  tCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACg                                                                                    <  2:1104894/98‑1 (MQ=255)
                   gcgcgcGATATACTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGc                                                                                 >  1:909318/6‑100 (MQ=255)
                   cGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTcccagatcggaag                                                                                 >  1:4081125/1‑90 (MQ=255)
                   cGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGc                                                                                 >  1:4101774/1‑100 (MQ=255)
                    gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttttgctca                                                                                <  1:486274/100‑10 (MQ=37)
                    gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                                >  2:1573810/1‑100 (MQ=255)
                    gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCt                                                                                >  2:1286372/1‑100 (MQ=255)
                     ccACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTa                                                                               <  1:3951305/100‑1 (MQ=255)
                       aCGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAat                                                                             >  1:249349/1‑100 (MQ=255)
                       aCGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAat                                                                             <  1:2354550/100‑1 (MQ=255)
                          atatAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAattat                                                                          >  1:3344534/1‑100 (MQ=255)
                           tagaGTATATCGTGCAAGGCTGTGAGGTTGGTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttttgc                                                                                   <  2:1014249/87‑7 (MQ=21)
                              aGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCTTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttt                                                                       >  1:2838892/1‑99 (MQ=255)
                               gTATATCGCTCGCGGCTGTGAGGGTGTTTGCCTGTGTTCTCCTTATGCGCAAAAAAGGGTAGTGGCAGGCCACTTGCCTTGTACGAGCTGATTAttttt                                                                      >  2:832148/1‑99 (MQ=255)
                                tataTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTg                                                                    <  1:548038/100‑1 (MQ=255)
                                 atatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGc                                                                   <  2:1129344/100‑1 (MQ=255)
                                  tatCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCt                                                                  >  2:3054185/1‑100 (MQ=255)
                                  tatCGCGCGCGGCGGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAACAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTgat                                                                  >  2:2502974/1‑98 (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)
                                          gcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTCCCCTTGTACGAGCTAATTATTTTTTGCTTTcttctt                                                          >  2:3394325/1‑100 (MQ=255)
                                           cgacTGTCAAGGTGCTTTCATGTGTACTCCTTATGAGCAAAATATGGATGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTATCTGCTttcttc                                                            <  2:1658835/94‑1 (MQ=255)
                                           cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                         >  1:2020548/1‑100 (MQ=255)
                                           cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                         <  1:1901291/100‑1 (MQ=255)
                                           cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                         >  1:777340/1‑100 (MQ=255)
                                            ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                        <  1:1451647/100‑1 (MQ=255)
                                            ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                        <  2:667585/100‑1 (MQ=255)
                                            ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                        <  1:1208697/100‑1 (MQ=255)
                                             gcggtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                       <  2:1716431/97‑1 (MQ=255)
                                             gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCa                                                       >  2:3389929/1‑100 (MQ=255)
                                              cTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                      >  2:2114980/1‑100 (MQ=255)
                                               tgtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTCCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAAt                                                     <  1:3394345/100‑1 (MQ=255)
                                               tgagagaGTGATTTCATATGTTTTTCTTATGAACAAAAAAGGGAAGTGGCAAACCACTTCCATTGTAAGAGCTAATTATTTTTTGCTTTATTATTTCaa                                                      >  2:1937717/4‑99 (MQ=255)
                                                gtgAGGGTGTTTTCATGTGTTCTCCTTATGTGCAAAAGAGGTAAGTTGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTAAATa                                                    <  2:1852879/100‑1 (MQ=255)
                                                  gAGTGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGTAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                  <  2:3807531/100‑1 (MQ=255)
                                                  gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                  >  1:1048137/1‑100 (MQ=255)
                                                  gAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                  >  1:3615082/1‑100 (MQ=255)
                                                  gAGCGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAATGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                  <  2:2899121/100‑1 (MQ=255)
                                                    gggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCtt                                                >  2:4048772/1‑100 (MQ=255)
                                                     ggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCCAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTa                                               <  1:2660006/100‑1 (MQ=255)
                                                     ggTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTa                                               >  1:2090399/1‑100 (MQ=255)
                                                         ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTTCTTTTTTTTTTCAATACATTAACgg                                           >  2:706909/1‑100 (MQ=255)
                                                           ttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGta                                         >  1:551986/1‑100 (MQ=255)
                                                           tgcATGTTTTCGCCTTAGGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACgg                                           <  2:1360261/96‑1 (MQ=255)
                                                                gtgtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCg                                    <  1:1769051/100‑1 (MQ=255)
                                                                  gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                  >  1:357104/1‑100 (MQ=255)
                                                                  gtTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                  >  1:1219282/1‑100 (MQ=255)
                                                                   ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGcc                                 >  2:4104593/1‑100 (MQ=255)
                                                                    tctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTCTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCa                                >  1:3909722/1‑100 (MQ=255)
                                                                     ctcCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCAt                               <  2:1765085/100‑1 (MQ=255)
                                                                       ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCa                             >  1:1729530/1‑100 (MQ=255)
                                                                          tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                          <  2:336631/100‑1 (MQ=255)
                                                                          tATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                          <  1:2318025/100‑1 (MQ=255)
                                                                          cgatctcAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGcc                          <  1:3850922/94‑1 (MQ=255)
                                                                            tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                        >  2:2629375/1‑100 (MQ=255)
                                                                            tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                        >  1:1137162/1‑100 (MQ=255)
                                                                             gctcgtacaaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGa                       >  2:391725/9‑100 (MQ=37)
                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                     >  1:662782/1‑100 (MQ=255)
                                                                               gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                     >  2:2701798/1‑100 (MQ=255)
                                                                                cAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCAGatcg                    >  2:3850897/1‑97 (MQ=255)
                                                                                  aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTa                  >  1:295211/1‑100 (MQ=255)
                                                                                    naaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTAt                 >  2:2949830/2‑99 (MQ=255)
                                                                                    caaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg                >  2:3572536/2‑100 (MQ=255)
                                                                                    aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg                >  1:3182763/1‑100 (MQ=255)
                                                                                    aaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATg                <  2:3257582/100‑1 (MQ=255)
                                                                                     aaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGc               <  1:2287641/100‑1 (MQ=255)
                                                                                       gggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAc             >  1:2229538/1‑100 (MQ=255)
                                                                                        ggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCAcc            >  1:3749262/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       <  2:2034254/100‑1 (MQ=255)
                                                                                                 aaGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATc   <  2:2027916/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
AGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCG  >  NC_000913/4356348‑4356545

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: