Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,196,204 G→A intergenic (‑105/+128) thiC ← / ← rsd phosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,196,2040GA99.3% 507.5 / ‑4.7 144intergenic (‑105/+128)thiC/rsdphosphomethylpyrimidine synthase/stationary phase protein, binds sigma 70 RNA polymerase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (74/69);  total (74/70)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 147 total aligned reads displayed.

TCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATT  >  NC_000913/4196107‑4196300
                                                                                                 |                                                                                                
tCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTcaca                                                                                                <  2:2023186/100‑1 (MQ=255)
 ccAAAAAGTTAAGGACGTGGCTTGTCAGACGGCGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAg                                                                                               >  2:134668/1‑100 (MQ=255)
 ccAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTcaca                                                                                                <  2:2126878/99‑1 (MQ=255)
 ccAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAg                                                                                               >  1:650776/1‑100 (MQ=255)
  cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAgg                                                                                              >  1:1097116/1‑100 (MQ=255)
  cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAgg                                                                                              >  2:1960218/1‑100 (MQ=255)
  cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAgg                                                                                              >  2:3428859/1‑100 (MQ=255)
  cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAgg                                                                                              >  1:2700995/1‑100 (MQ=255)
  cAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAga                                                                                              >  1:3853633/1‑99 (MQ=255)
   aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGt                                                                                             >  1:3182141/1‑100 (MQ=255)
   aaaaaGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGt                                                                                             >  2:3067186/1‑100 (MQ=255)
         ttAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGc                                                                                       >  1:2214489/1‑100 (MQ=255)
          tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGcc                                                                                      <  1:172323/100‑1 (MQ=255)
          tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGcc                                                                                      <  2:3792046/100‑1 (MQ=255)
          tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGcc                                                                                      <  1:1611169/100‑1 (MQ=255)
          tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGcc                                                                                      <  1:1611046/100‑1 (MQ=255)
          tAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGcc                                                                                      <  1:2538754/100‑1 (MQ=255)
           aaGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCt                                                                                     <  2:3155672/100‑1 (MQ=255)
           aaGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCt                                                                                     <  1:3555113/100‑1 (MQ=255)
             ggACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGa                                                                                   >  2:2359293/1‑100 (MQ=255)
              gACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGAt                                                                                  >  1:952312/1‑100 (MQ=255)
               aCGTGGCTTGTCAGATGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATc                                                                                 <  1:3813063/100‑1 (MQ=255)
               aCGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATc                                                                                 <  2:3642586/100‑1 (MQ=255)
               aCGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATc                                                                                 <  2:409303/100‑1 (MQ=255)
               aCGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATc                                                                                 <  1:1012803/100‑1 (MQ=255)
                cGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCa                                                                                <  1:1805777/100‑1 (MQ=255)
                cGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCa                                                                                <  2:2426716/100‑1 (MQ=255)
                 gTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAg                                                                               <  2:3280306/100‑1 (MQ=255)
                  tGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGTGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAg                                                                               <  1:2164072/99‑1 (MQ=255)
                  tGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAgg                                                                              >  2:1767456/1‑100 (MQ=255)
                   ggCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGt                                                                             <  1:2709677/100‑1 (MQ=255)
                   ggCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGt                                                                             >  2:603592/1‑100 (MQ=255)
                   ggCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCATGTTCCCTTCACAGGTATTAGCCTGATCAGGt                                                                             >  1:3180860/1‑100 (MQ=255)
                       tGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCcg                                                                         <  2:928855/100‑1 (MQ=255)
                        gTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCcgc                                                                        <  2:2529356/100‑1 (MQ=255)
                         tCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCcgag                                                                       <  1:2625571/100‑3 (MQ=255)
                                 cgGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGa                                                               <  2:3209880/100‑1 (MQ=255)
                                  ggATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGaa                                                              <  1:912952/100‑1 (MQ=255)
                                   gATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAAt                                                             >  2:2456604/1‑100 (MQ=255)
                                   gATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAAt                                                             >  1:24294/1‑100 (MQ=255)
                                      gAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATaaa                                                          >  1:3702716/1‑100 (MQ=255)
                                       aaGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAAc                                                         <  2:935589/100‑1 (MQ=255)
                                       aaGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAAc                                                         <  2:3196843/100‑1 (MQ=255)
                                       aaGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAAc                                                         >  1:420924/1‑100 (MQ=255)
                                        aGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACg                                                        <  1:21133/100‑1 (MQ=255)
                                         gCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACgg                                                       <  1:2574258/100‑1 (MQ=255)
                                          cAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGt                                                      >  2:743694/1‑100 (MQ=255)
                                             gagaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTTCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtctc                                                   >  1:169442/1‑100 (MQ=255)
                                             gagaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAAGTTCCGCGGATCCCGAATAAACGGtctc                                                   >  1:130237/1‑100 (MQ=255)
                                             gagaCGATCGCCGCAGGGGAGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGtctc                                                   >  2:2101168/1‑100 (MQ=255)
                                              agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCa                                                  >  1:2422726/1‑100 (MQ=255)
                                              agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCa                                                  <  1:3735123/100‑1 (MQ=255)
                                              agaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCa                                                  <  1:3316538/100‑1 (MQ=255)
                                               gaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAg                                                 <  2:668910/100‑1 (MQ=255)
                                               gaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAg                                                 >  1:2817266/1‑100 (MQ=255)
                                               gaCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAg                                                 <  2:3172874/100‑1 (MQ=255)
                                                aCGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGc                                                <  1:1474967/100‑1 (MQ=255)
                                                 cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGcc                                               <  1:2359301/100‑1 (MQ=255)
                                                 cGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGc                                                >  2:2192066/1‑99 (MQ=255)
                                                  gATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAgatc                                              >  1:2146750/1‑97 (MQ=255)
                                                   aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAg                                             <  1:1960224/100‑1 (MQ=255)
                                                   aTCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAACAAACGGTCTCAGCCAg                                             >  2:2259056/1‑100 (MQ=255)
                                                      gccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTa                                          >  1:2633026/1‑100 (MQ=255)
                                                      gccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTa                                          <  1:2062512/100‑1 (MQ=255)
                                                      gccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTa                                          <  2:1836487/100‑1 (MQ=255)
                                                       ccgcAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTaa                                         <  2:1320043/100‑1 (MQ=255)
                                                            ggggCGATGCGATAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCAc                                    >  2:1522293/1‑100 (MQ=255)
                                                                cgatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCg                                <  2:2812734/100‑1 (MQ=255)
                                                                 gatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGa                               <  1:1263696/100‑1 (MQ=255)
                                                                 gatgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGa                               >  2:2666069/1‑100 (MQ=255)
                                                                   tgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACa                             <  1:2976856/100‑1 (MQ=255)
                                                                   tgcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACa                             <  1:2036886/100‑1 (MQ=255)
                                                                    gcgatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACaa                            >  2:614950/1‑100 (MQ=255)
                                                                      gatAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGa                          >  1:3244352/1‑100 (MQ=255)
                                                                       atAGCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGaa                         >  2:266126/1‑100 (MQ=255)
                                                                          gCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAtt                      <  2:217534/100‑1 (MQ=255)
                                                                          gCAGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAAtt                      >  1:3509231/1‑100 (MQ=255)
                                                                            aGATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAg                    <  2:3033042/100‑1 (MQ=255)
                                                                             gATTAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGc                   >  2:2889663/1‑100 (MQ=255)
                                                                               ttAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGccc                 <  1:1073114/100‑1 (MQ=255)
                                                                                tAACTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGcccc                <  1:1922228/100‑1 (MQ=255)
                                                                                 aaCTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCg               >  1:425944/1‑100 (MQ=255)
                                                                                  aCTCTTGTTTCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGc              <  2:169441/100‑1 (MQ=255)
                                                                                  aCTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGc              >  1:1423193/1‑100 (MQ=255)
                                                                                  aCTCTTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGc              >  1:3635974/1‑100 (MQ=255)
                                                                                   ctctTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCa             >  1:1752569/1‑100 (MQ=255)
                                                                                     ctTGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCaaa           >  2:3436467/1‑100 (MQ=255)
                                                                                      ttGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAAt          >  1:2383040/1‑100 (MQ=255)
                                                                                       tGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATg         <  2:2817251/100‑1 (MQ=255)
                                                                                       tGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATg         >  1:3037563/1‑100 (MQ=255)
                                                                                       tGTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATg         <  1:2948559/100‑1 (MQ=255)
                                                                                        gTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgg        >  1:532914/1‑100 (MQ=255)
                                                                                        gTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgg        >  2:1577305/1‑100 (MQ=255)
                                                                                        gTTCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgg        >  1:3250443/1‑100 (MQ=255)
                                                                                         ttCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATggg       <  2:2422717/100‑1 (MQ=255)
                                                                                          tCCCTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATgggg      >  1:675602/1‑100 (MQ=255)
                                                                                           cccTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGc     >  2:1955658/1‑100 (MQ=255)
                                                                                           cccTTCACAGGTATTAGCCTGAGCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGc     <  1:2215022/100‑1 (MQ=255)
                                                                                            ccTTCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCa    >  2:1746411/1‑100 (MQ=255)
                                                                                              ttCACAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCAtt  >  1:2911421/1‑100 (MQ=255)
                                                                                                 |                                                                                                
TCCAAAAAGTTAAGGACGTGGCTTGTCAGACGACGGATGAAGCAAGAGACGATCGCCGCAGGGGCGATGCGATAGCAGATTAACTCTTGTTCCCTTCGCAGGTATTAGCCTGATCAGGTTCCGCGGATCCCGAATAAACGGTCTCAGCCAGTTAAGGCACTCCGACAAGAAATTAGCCCCGCAAATGGGGCATT  >  NC_000913/4196107‑4196300

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: