Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,357:1 +T 24.5% intergenic (+151/+21) ydbL → / ← feaR DUF1318 family protein/transcriptional activator for tynA and feaB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,446,3571.T24.5% 257.0 / 114.4 106intergenic (+151/+21)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base . (44/36);  new base T (11/15);  total (55/51)
Fisher's exact test for biased strand distribution p-value = 3.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.40e-01

Only 100 of 115 total aligned reads displayed.

TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCA  >  NC_000913/1446258‑1446456
                                                                                                    |                                                                                                   
ttttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCa                                                                                                      <  1:2240312/100‑1 (MQ=255)
ttttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCa                                                                                                      <  1:1600379/100‑1 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑g                                                                                                    >  1:205424/1‑100 (MQ=255)
  ttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑Gt                                                                                                   >  1:1081253/1‑100 (MQ=255)
    ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTAc                                                                                                 >  1:2076629/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCatgtac                                                                                                 >  1:2389757/1‑95 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                                >  2:652604/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATAACTGTTTTAACTAATTGGCGTTGCA‑g                                                                                                    >  2:1577225/1‑96 (MQ=255)
      tGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACAt                                                                                               >  1:812214/1‑100 (MQ=255)
       gAAACAACGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTTAAAAATGCATAGTAAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATg                                                                                              <  1:1461613/100‑1 (MQ=255)
         aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATATGTTTTAACTAATTGGCGTTGCA‑GTACATGCa                                                                                            >  2:22470/1‑100 (MQ=255)
          aCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCCGTTTTAACTAATTGGCGTTGCA‑GTACATGCaa                                                                                           <  2:1836987/100‑1 (MQ=255)
            aaGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACg                                                                                         <  1:271216/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                       <  1:699050/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                       <  2:288602/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                       <  2:1241657/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                       <  2:1500246/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                       <  2:1818160/100‑1 (MQ=38)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGcc                                                                                       >  1:1719470/1‑100 (MQ=255)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCa                                                                                      <  2:1129508/100‑1 (MQ=38)
               ccgatctttCCAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAACATGCATAGTCAAATATCTGTTTTAGCTAATTGGCGTTGCATGTAC‑TGCAACGCCa                                                                                      <  1:538278/93‑1 (MQ=25)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                     <  1:1924958/100‑1 (MQ=255)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                     <  1:844651/100‑1 (MQ=255)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                    <  2:2220686/100‑1 (MQ=38)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                    <  2:2288746/100‑1 (MQ=38)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                    <  2:327966/100‑1 (MQ=38)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAt                                                                                    <  1:2573092/100‑1 (MQ=255)
                  aaTATTCAAAAACTGCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtt                                                                                   >  1:2239186/1‑100 (MQ=255)
                  aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtt                                                                                   <  2:738304/100‑1 (MQ=255)
                  aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAAGGCATAGTAAAATATCTGGTTTAACTAATTGGCGTTGCA‑GTACATGCAACGGCAAtt                                                                                   <  1:1166691/100‑1 (MQ=255)
                    tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttag                                                                                 <  1:1465510/100‑1 (MQ=255)
                    tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttag                                                                                 <  2:317710/100‑1 (MQ=255)
                     atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttagt                                                                                <  1:1097158/100‑1 (MQ=37)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                               <  1:1204066/100‑1 (MQ=255)
                        cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                             >  2:1296646/1‑100 (MQ=255)
                        cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                             <  2:1393288/100‑1 (MQ=255)
                        cAAAAACTCCTGTAAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                             >  2:1979309/1‑100 (MQ=255)
                          aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCg                                                                           <  2:337281/100‑1 (MQ=255)
                               tCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCGGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAt                                                                       <  1:277742/99‑1 (MQ=255)
                                 cTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTa                                                                    >  2:923751/1‑100 (MQ=255)
                                   gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                  >  2:2891785/1‑100 (MQ=255)
                                   gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGCAATTTACg                                                                  >  1:2834740/1‑100 (MQ=255)
                                     cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                                <  2:2118103/100‑1 (MQ=255)
                                      aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                               <  1:791655/100‑1 (MQ=255)
                                      aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCg                                                               >  2:378344/1‑100 (MQ=255)
                                         ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                            <  2:1053524/100‑1 (MQ=255)
                                         ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                            >  2:1706478/1‑100 (MQ=255)
                                          tGCCTTTTGCCCAGAAAAATGCATAGTCAAGTATCTGGTTTAACTAATTGGCGTTACA‑GTACATGCAACGCAAATGAGTTAGCGTAATTTACATAGATa                                                            >  1:2811678/1‑99 (MQ=255)
                                           gCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACt                                                          >  2:2496428/1‑100 (MQ=255)
                                           gCCTTGTGCCCTGAAAAATGCAGAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGAt                                                             <  1:1447459/97‑1 (MQ=255)
                                             cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                        >  2:2906098/1‑100 (MQ=255)
                                                  gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                   <  2:167830/100‑1 (MQ=255)
                                                  gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                   >  2:2550698/1‑100 (MQ=255)
                                                  gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                   >  2:315653/1‑100 (MQ=255)
                                                  gCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTgg                                                   >  2:1702466/1‑100 (MQ=255)
                                                   cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGc                                                  >  2:338817/1‑100 (MQ=255)
                                                   cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑ATACATGCAACGCCAAGTAGGTAGCGGAGTTTACGTCGATACGCGCCTGGc                                                  >  1:2442367/1‑100 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                 >  1:1190817/1‑100 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAATTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                 <  2:2764080/100‑1 (MQ=255)
                                                      tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                 <  1:2118331/98‑1 (MQ=255)
                                                       gAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                               >  2:821274/1‑99 (MQ=255)
                                                         aaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATc                                            >  1:2606128/1‑100 (MQ=255)
                                                            aggCATAGTCAAATATCTGTGTAAACGATTTGGCGTTGCA‑GTACAGGCAACGCCAAGTAGTTAGCGGAATTTACGTCGATACTCGCCTTGACTCATcc                                           <  1:2146452/96‑1 (MQ=255)
                                                               cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAg                                      >  1:2438949/1‑100 (MQ=255)
                                                                aTAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGc                                     >  2:704099/1‑100 (MQ=255)
                                                                 tAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGACATCACAAAGCg                                    <  1:332931/100‑1 (MQ=255)
                                                                   gTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                  >  2:2000800/1‑100 (MQ=255)
                                                                   gTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                  >  2:1492250/1‑100 (MQ=255)
                                                                   gTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                  <  2:1393877/100‑1 (MQ=255)
                                                                    tCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTg                                 <  2:1489906/100‑1 (MQ=255)
                                                                        atatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                                <  1:983275/97‑1 (MQ=255)
                                                                         tatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaa                            >  2:2107298/1‑100 (MQ=255)
                                                                          atCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                           >  1:2354578/1‑100 (MQ=255)
                                                                           tCTGTTTTAACTAATTCGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                          <  1:1702459/100‑1 (MQ=255)
                                                                             tGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAc                        >  1:1167766/1‑100 (MQ=255)
                                                                              ccgctAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTAGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                       >  1:1948917/5‑100 (MQ=25)
                                                                                 ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGtt                    <  2:437083/100‑1 (MQ=255)
                                                                                 ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCACAAAGCGTTGCTTAAATACCGtt                    >  1:2138654/1‑100 (MQ=255)
                                                                                  tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                   <  2:347916/100‑1 (MQ=255)
                                                                                  tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                   <  2:595313/100‑1 (MQ=255)
                                                                                  tAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTg                   <  2:2624000/100‑1 (MQ=255)
                                                                                    aCTAATTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaa                 >  1:632699/1‑100 (MQ=37)
                                                                                     cTAATTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                >  2:2505977/1‑100 (MQ=37)
                                                                                     cTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa                <  2:1022855/100‑1 (MQ=255)
                                                                                       aaTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa              >  2:300746/1‑100 (MQ=38)
                                                                                       aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa              >  1:1488491/1‑100 (MQ=255)
                                                                                       aaTTGGCGTTGAATGTAC‑TGCGAGGCCAATTAGTTCACGGAATTTACAACGATACCAACCTGGAGTCATAACAAACCGCTGCATAAATACCGTTGaaa                >  1:1005298/1‑98 (MQ=255)
                                                                                         ttGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg            <  1:2529849/100‑1 (MQ=255)
                                                                                           ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          >  1:137396/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          >  2:2435473/1‑100 (MQ=38)
                                                                                           ggCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc          >  1:2077918/1‑100 (MQ=255)
                                                                                            gCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAct         >  2:10323/1‑100 (MQ=255)
                                                                                             cGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActc        >  2:2481381/1‑100 (MQ=38)
                                                                                             cGTTGCA‑GTACATGCAACGCCAAATAGTTAACGGAATTTACGTCGATACACGCCTGGCGTAATCACAAAGCGTGGCTTAAATACAGTTGAAAAATGAc          >  1:2768112/1‑98 (MQ=255)
                                                                                                tGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg     <  2:1354062/100‑1 (MQ=255)
                                                                                                 gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTgg     >  1:1438382/1‑99 (MQ=255)
                                                                                                 gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt    >  2:177815/1‑100 (MQ=255)
                                                                                                   a‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa  >  1:1718893/1‑100 (MQ=255)
                                                                                                   a‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCa  <  2:462093/100‑1 (MQ=255)
                                                                                                   a‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAACGTTGAAAAATGACTCTGGTCa  <  1:776648/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTCA  >  NC_000913/1446258‑1446456

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: ‑