Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,446,357:1 +T 16.3% 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.T16.3% 441.4 / 92.0 135intergenic (+151/+21)ydbL/feaRDUF1318 family protein/transcriptional activator for tynA and feaB
Reads supporting (aligned to +/- strand):  ref base . (49/64);  new base T (11/11);  total (60/75)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

Only 100 of 142 total aligned reads displayed.

TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTC  >  NC_000913/1446258‑1446455
                                                                                                    |                                                                                                  
 tttCCTGAAACAAGGTGAATATTCAAAATCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑g                                                                                                   <  2:1384428/100‑1 (MQ=255)
 tttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑g                                                                                                   <  1:460281/100‑1 (MQ=255)
  ttCGTGAAATCAGGTGAATAATCAAACACTCCTGTCAAATTGCGTTTTGCACTGAAAAATGCGTAATAAGATATCTGTTTTAACTAATTTTCGTTGCA‑g                                                                                                   <  2:3155738/99‑1 (MQ=255)
  ttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑Gt                                                                                                  <  1:3209/100‑1 (MQ=255)
  ttCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑Gt                                                                                                  <  1:857461/100‑1 (MQ=255)
    ccTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTAc                                                                                                >  2:3292687/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGTCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                               >  2:54658/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTAc                                                                                                >  2:1466258/1‑99 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                               >  1:2772538/1‑100 (MQ=255)
     cTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACa                                                                                               <  1:2144255/100‑1 (MQ=255)
     atGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTAc                                                                                                >  1:949700/2‑99 (MQ=255)
         aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCa                                                                                           >  1:2115290/1‑100 (MQ=255)
         aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCa                                                                                           >  2:2265425/1‑100 (MQ=255)
         aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCa                                                                                           >  1:3400048/1‑100 (MQ=255)
         aaCAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCa                                                                                           <  1:870982/100‑1 (MQ=255)
           cAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAAc                                                                                         <  1:1969363/100‑1 (MQ=255)
           cAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAAc                                                                                         <  1:3596005/100‑1 (MQ=255)
           cAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAAc                                                                                         <  2:3456288/100‑1 (MQ=255)
            aaGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACg                                                                                        <  2:463241/100‑1 (MQ=255)
             aGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGc                                                                                       <  1:321985/100‑1 (MQ=255)
              ggTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                      <  1:1902704/100‑1 (MQ=38)
               gTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCa                                                                                     <  1:1030449/100‑1 (MQ=255)
                tGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCaa                                                                                    >  1:1658590/1‑100 (MQ=255)
                tGAATAGTCAAAGACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCaa                                                                                    <  2:2263794/100‑1 (MQ=37)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                   <  1:473401/100‑1 (MQ=38)
                 gAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAt                                                                                   <  1:1402320/100‑1 (MQ=38)
                  aaTATTCAAAATCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAt                                                                                   <  2:2996087/99‑1 (MQ=255)
                  aaTATTCAAAAACTCCTGTTAAATTGCCTAAAAAAAAAAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCa                                                                                     <  2:1246916/97‑1 (MQ=16)
                  aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTACTTGGCGTTGCA‑GTACATGCAACGCCAAtt                                                                                  >  1:75987/1‑100 (MQ=255)
                  aaTATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAtt                                                                                  >  2:2777211/1‑100 (MQ=255)
                  aaTATTCAAAAACTCCTGTCAAATTCGCTTTTGCCCTGAAAAATGCATAGTCAAAAATCTGTCTTAACTAATTGGCGTTGCATGTAC‑TGCAACGcc                                                                                      <  2:1215089/96‑1 (MQ=25)
                    tatTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttag                                                                                <  1:1772670/100‑1 (MQ=255)
                     atTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagt                                                                               <  1:1474118/100‑1 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                              >  1:1069711/1‑100 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                              <  1:1664809/100‑1 (MQ=255)
                      ttCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtt                                                                              <  1:1949176/100‑1 (MQ=255)
                        cAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagttag                                                                            >  1:1368380/1‑100 (MQ=255)
                        cAAAAACTCCTGTCAAATTGCCTTTTGCACTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAttagtta                                                                             >  2:1988115/1‑99 (MQ=255)
                         aaaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGc                                                                           >  1:1966453/1‑100 (MQ=255)
                          aaaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttagttaaaa                                                                          <  1:387419/100‑4 (MQ=25)
                           aaaCTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCATGTAC‑TGCAACGCCAAttagttaaaac                                                                         <  2:3153426/100‑5 (MQ=25)
                                ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                    >  1:1244419/1‑100 (MQ=255)
                                ccTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAAttt                                                                    >  1:2410451/1‑100 (MQ=255)
                                  tGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTAc                                                                  >  2:2698045/1‑100 (MQ=255)
                                   gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                 <  1:2566210/100‑1 (MQ=255)
                                   gTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                 <  1:215039/100‑1 (MQ=255)
                                     cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAGTTTACGTc                                                               >  2:1153884/1‑100 (MQ=255)
                                     cAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                               <  1:1664973/100‑1 (MQ=255)
                                     cAAATTGACTTGTGCCTTGAAAAATGCATCGTCAATTATCTGTTTTAACTGATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACg                                                                 <  2:3568248/98‑1 (MQ=38)
                                      aaaTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTc                                                               <  2:3448077/99‑1 (MQ=255)
                                        aTTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCTTTGCA‑GTACATGCAACGCCAACTAGTTAGCGGACTTTACGTCGAt                                                            >  2:1046246/1‑100 (MQ=255)
                                         ttGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                           <  1:3295420/100‑1 (MQ=255)
                                             cTTTTGCCCTGAAAAATGTATAGGTCAATATCGGTTTTAACTAATTGGCGTTGCA‑TTACATTCAACGCCAATTAGTTAGCGGAATTTACGTCGATa                                                           <  2:1692633/96‑1 (MQ=37)
                                             cTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCg                                                       >  2:3565831/1‑100 (MQ=255)
                                                 tGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTg                                                   <  2:1252780/100‑1 (MQ=255)
                                                   cccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGc                                                 >  1:1228874/1‑100 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  2:437764/100‑1 (MQ=255)
                                                    ccTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCg                                                <  2:2115277/100‑1 (MQ=255)
                                                     cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                               >  1:2023490/1‑100 (MQ=255)
                                                     cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                               <  2:1452359/100‑1 (MQ=255)
                                                     cTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAAATAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGt                                               >  1:2553165/1‑100 (MQ=255)
                                                      tGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTc                                              >  1:915483/1‑100 (MQ=255)
                                                       gAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCa                                             >  2:111828/1‑100 (MQ=255)
                                                        aaaaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCAt                                            <  2:834737/100‑1 (MQ=255)
                                                          aaaTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATcc                                          >  2:570407/1‑100 (MQ=255)
                                                            aTGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCa                                        >  1:1508657/1‑100 (MQ=255)
                                                              gCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                      >  1:2804801/1‑100 (MQ=255)
                                                               cATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCaaa                                      >  2:2283305/1‑99 (MQ=255)
                                                                  aGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCACTTAGTTAGCGGAATGTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGt                                  <  2:907823/100‑1 (MQ=255)
                                                                  aGTCAAATATCTGTTTTAACGAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGt                                  <  1:3367719/100‑1 (MQ=255)
                                                                   gTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                 >  2:631328/1‑100 (MQ=255)
                                                                   gTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGtt                                 <  1:2396026/100‑1 (MQ=255)
                                                                     cAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGc                               <  1:838663/100‑1 (MQ=255)
                                                                      aaaTATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCt                              >  1:3409837/1‑100 (MQ=255)
                                                                       aaTATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCtt                             <  1:1321381/100‑1 (MQ=255)
                                                                        atatCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTa                            <  1:905947/100‑1 (MQ=255)
                                                                          atCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTaaa                          >  1:2182589/1‑100 (MQ=255)
                                                                           tCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCAAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAAt                         >  2:2241064/1‑100 (MQ=255)
                                                                            cTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGGTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATa                        <  2:1780239/100‑1 (MQ=255)
                                                                              gTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATAcc                      >  2:252438/1‑100 (MQ=255)
                                                                               ttttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCg                     >  1:2636502/1‑100 (MQ=255)
                                                                                 ttAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGt                    <  1:3628660/99‑1 (MQ=255)
                                                                                     cTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaa               >  1:1435661/1‑100 (MQ=255)
                                                                                       aaTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa             >  2:2083876/1‑100 (MQ=38)
                                                                                       aaTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGaaaaa             >  1:3056318/1‑100 (MQ=255)
                                                                                        aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            <  1:1786745/100‑1 (MQ=38)
                                                                                        aTTGGCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAAt            >  2:1786737/1‑100 (MQ=38)
                                                                                        aTTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAACTACCGTTGAAAAAt            >  2:2308006/1‑100 (MQ=255)
                                                                                         ttGGCGTTGCA‑GTTCATGCAACGTCAATTAGTTAGCGGAATTTACGTAGATCTTGGCATGGCGTCATCCCAAAGTGTTACTTATATACTTTTGCAACATg           >  2:88346/1‑100 (MQ=255)
                                                                                         ttGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATg           <  2:2513306/100‑1 (MQ=255)
                                                                                          tGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGa          <  1:111828/100‑1 (MQ=255)
                                                                                           ggCGTTGCATGTAC‑TGCAACGCCAATTGGTTAGCGGAATGTACGTCGATACTCGCCTGGCGTCAGCACAAAGCGTTGCTTAAATACCGTTGAAAAATg           >  2:898013/1‑98 (MQ=37)
                                                                                           ggCGTTGCATGTAC‑TGCAACGCCAATTCGTTAGCGGAATTTACGTCGATACTCGCCTGGCGCCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  2:1573552/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  1:2616278/1‑100 (MQ=38)
                                                                                           ggCGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATACCAAAGCGTTGCTTAAATACCGTTGaaaaa             >  1:1974696/1‑96 (MQ=38)
                                                                                            gCGTTGCATGTAC‑TGCAACGCTAATTAGTTAGCGGAATTTACGTCGATACTCGCGTGGCGTCATCTCAAAACATTGCTTAAATAAAGTTGAACAATGa          >  2:2160706/1‑98 (MQ=21)
                                                                                             cGTTGCATGTAC‑TGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGAc         >  1:2546538/1‑98 (MQ=38)
                                                                                              gTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGActct      >  1:1435205/1‑100 (MQ=255)
                                                                                                 gCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGt   <  2:1824265/100‑1 (MQ=255)
                                                                                                  cA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTc  <  1:801606/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
TTTTCCTGAAACAAGGTGAATATTCAAAAACTCCTGTCAAATTGCCTTTTGCCCTGAAAAATGCATAGTCAAATATCTGTTTTAACTAATTGGCGTTGCA‑GTACATGCAACGCCAATTAGTTAGCGGAATTTACGTCGATACTCGCCTGGCGTCATCCCAAAGCGTTGCTTAAATACCGTTGAAAAATGACTCTGGTC  >  NC_000913/1446258‑1446455

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