Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 5.0% noncoding (227/245 nt) csrC → CsrC sRNA sequesters CsrA, a carbon flux regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,051,2620AC 5.0% 1764.5 / 54.4 506noncoding (227/245 nt)csrCCsrC sRNA sequesters CsrA, a carbon flux regulator
Reads supporting (aligned to +/- strand):  ref base A (244/234);  new base C (19/8);  total (263/243)
Fisher's exact test for biased strand distribution p-value = 7.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 515 total aligned reads displayed.

AGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGA  >  NC_000913/4051114‑4051411
                                                                                                                                                    |                                                                                                                                                     
aGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCTTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                                                                      >  2:2808884/1‑150 (MQ=255)
  gAAGACAAACGTACGGACACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTACAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCt                                                                                                                                                    >  2:6995304/1‑150 (MQ=255)
    aGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGCGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                                                                                  >  2:4264757/1‑150 (MQ=255)
     gACACACGTCCGGAGACGTAATTAAACGGAAATGGAACCACCACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTc                                                                                                                                                 <  2:7381946/150‑1 (MQ=255)
      aCAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCg                                                                                                                                                <  2:8675151/150‑1 (MQ=255)
        aaaCGTCCGTAGACTTAATTAAACGGAAATGTAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGTGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGcc                                                                                                                                              <  1:7359261/150‑1 (MQ=255)
        aaaCGTCCGGAGACTTAATTAAACGGAAATGTAATCAACACGGATTGTTCCGGCTAAAGTAAAAACTGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGcc                                                                                                                                              <  1:3404236/150‑1 (MQ=255)
        aaaCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGcc                                                                                                                                              >  2:4022791/1‑150 (MQ=255)
        aaaCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGcc                                                                                                                                              >  1:3328742/1‑150 (MQ=255)
             tCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttttt                                                                                                                                         >  2:8560417/1‑150 (MQ=255)
                  agaCGGAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGCAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCttt                                                                                                                                    <  2:4503431/150‑1 (MQ=255)
                   gaCGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                                                                   <  2:4565115/150‑1 (MQ=255)
                      gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgctt                                                                                                                                >  1:3664736/1‑150 (MQ=255)
                      gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgctt                                                                                                                                >  2:5631652/1‑150 (MQ=255)
                          ttAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                                                                            >  1:2368992/1‑150 (MQ=255)
                           tAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCttt                                                                                                                           >  2:4407163/1‑150 (MQ=255)
                             aaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACTGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCt                                                                                                                         <  1:5623698/150‑1 (MQ=255)
                             aaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCt                                                                                                                         >  1:7489233/1‑150 (MQ=255)
                             aaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCt                                                                                                                         >  1:163758/1‑150 (MQ=255)
                             aaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCt                                                                                                                         >  1:4184247/1‑150 (MQ=255)
                               cGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                                                                       <  1:1335020/150‑1 (MQ=255)
                               cGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                                                                       <  2:5659201/150‑1 (MQ=255)
                                 gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                                     <  1:8210669/150‑1 (MQ=255)
                                   aaTGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                                                                                   <  1:3516795/150‑1 (MQ=255)
                                      ggAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTc                                                                                                                >  1:8363747/1‑150 (MQ=255)
                                      ggAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTc                                                                                                                >  2:6496452/1‑150 (MQ=255)
                                        aaTCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                                                                              <  1:88414/150‑1 (MQ=255)
                                        aaTCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccg                                                                                                              >  2:6346452/1‑150 (MQ=255)
                                            aaCACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                                                                          >  2:1744920/1‑150 (MQ=255)
                                               acGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                                                                       >  1:4970543/1‑150 (MQ=255)
                                                 ggATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                                                                     <  2:6122646/150‑1 (MQ=255)
                                                 ggATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                                                                     >  2:664474/1‑150 (MQ=255)
                                                  gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCt                                                                                                    <  1:4060200/150‑1 (MQ=255)
                                                   aTTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCta                                                                                                   >  2:4942334/1‑150 (MQ=255)
                                                    ttGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtat                                                                                                  <  2:4110383/150‑1 (MQ=255)
                                                         ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                                                                             >  1:1830740/1‑150 (MQ=255)
                                                         ccGGATAAAGGAAAAACAGGGTGTTTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTATTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                                                                             <  1:17312/150‑1 (MQ=255)
                                                             cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaa                                                                                         >  1:6941681/1‑150 (MQ=255)
                                                               aaaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaaa                                                                                       <  2:6885592/150‑1 (MQ=255)
                                                                  ggAAAAACTGGGTGTGTTGGCGGCCTGCACGGTTTGTAAGACCCGTTACGGGTTATGAGTCTGTAAAAATGGCGACAGTGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                                                                    <  1:2655918/150‑1 (MQ=255)
                                                                   gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                                                                                   >  2:35627/1‑150 (MQ=255)
                                                                   gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                                                                                   >  1:6832044/1‑150 (MQ=255)
                                                                   gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGc                                                                                   >  1:5218295/1‑150 (MQ=255)
                                                                       aaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaa                                                                               <  2:3484603/150‑1 (MQ=255)
                                                                       aaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaa                                                                               >  1:4407543/1‑150 (MQ=255)
                                                                        aCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTaaaa                                                                              <  2:9127340/150‑1 (MQ=255)
                                                                            ggTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACg                                                                          <  2:8388404/150‑1 (MQ=255)
                                                                             gtgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                                                                         <  1:5259349/150‑1 (MQ=255)
                                                                             gtgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                                                                         >  1:1939393/1‑150 (MQ=255)
                                                                             gtgtgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGa                                                                         <  2:2691189/150‑1 (MQ=255)
                                                                               gtgtTTGCGGACTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTAATGAAAAAAGGCGACAGAGTAATCTTTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                                                                       <  2:2423927/150‑1 (MQ=255)
                                                                                            gCAAGTATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca                                                          <  1:3770756/150‑1 (MQ=255)
                                                                                            gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAATGCGACAGAGTAATCTGTCGCCTTTTTTATTTGCTTTCTTTCTGTTATATTCCGCCGCAAATCTATAATGAAAAAAACTGCTAAAAGACGAACAAAAATGAAAAca                                                          >  1:9126132/1‑150 (MQ=255)
                                                                                             cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGACCCATTATGAAAccct                                                         >  1:2067487/1‑148 (MQ=255)
                                                                                             cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTATTTGCTTGCTTTCTGTTAGATTACGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat                                                         >  1:101139/1‑150 (MQ=255)
                                                                                                 gATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCatcatc                                                     >  2:5246426/1‑150 (MQ=255)
                                                                                                  aTTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCt                                                    <  2:774116/150‑1 (MQ=255)
                                                                                                   ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCAGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCtt                                                   <  2:7163471/150‑1 (MQ=255)
                                                                                                    tGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTc                                                  <  1:7836978/150‑1 (MQ=255)
                                                                                                      tAAGACCAGTTAAGGGTTATGAGTAAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTATTCTTTGCTTACTTTCTGTAAGATTCCGCCGCAAAACTATACTGAATAAAACGGCTAAAACACGAACCATTATGAAACCATCATCTTCaa                                                >  2:5896853/1‑150 (MQ=255)
                                                                                                          aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTc                                            >  2:7652557/1‑150 (MQ=255)
                                                                                                           cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCa                                           <  1:2227009/150‑1 (MQ=255)
                                                                                                           cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCa                                           <  1:2916948/150‑1 (MQ=255)
                                                                                                           cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCa                                           <  1:4851706/150‑1 (MQ=255)
                                                                                                           cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCa                                           <  1:6610104/150‑1 (MQ=255)
                                                                                                            ccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCAc                                          >  2:1510179/1‑150 (MQ=255)
                                                                                                            ccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAACTGTCGCCATTTTTCTTTGCTTGCTTTCTGTTAGATTACGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAAGAACCAATATGAAACCATCATCTTCAAACTCAc                                          >  1:6299395/1‑150 (MQ=255)
                                                                                                             cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACg                                         <  2:4774465/150‑1 (MQ=255)
                                                                                                             cGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTAGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAAACTAAACTGAATAAAAAGGCGAAAATACGCAACATTATGAAACCATCATCTTAAAACTCACg                                         >  1:2521675/1‑150 (MQ=255)
                                                                                                               ttAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgca                                       >  2:5438721/1‑150 (MQ=255)
                                                                                                                    ggTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCGGCaag                                  >  2:8073724/1‑149 (MQ=255)
                                                                                                                     gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATATGTCGCCTTTTTTATTTGCTTGCATTCTGTTAGATTACGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAACAATATGAAAACATCAAAATAAAAATCACGCAGCAAAg                                 >  1:7577326/1‑150 (MQ=255)
                                                                                                                       tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGt                               <  2:3027643/150‑1 (MQ=255)
                                                                                                                        aTGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTATTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTc                              <  1:724155/150‑1 (MQ=255)
                                                                                                                          gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCAt                            <  2:2342402/150‑1 (MQ=255)
                                                                                                                            gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGc                          <  2:7560782/150‑1 (MQ=255)
                                                                                                                             tCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTATTTTCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAATATGAAACCATCATCTTCAAACACACGCAGCAAAGGTCATGCa                         >  1:8093257/1‑150 (MQ=255)
                                                                                                                               aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa                       >  2:1912492/1‑150 (MQ=255)
                                                                                                                                 gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAg                     <  2:163758/150‑1 (MQ=255)
                                                                                                                                  aaaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgc                    <  1:1493124/150‑1 (MQ=255)
                                                                                                                                    aaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgcgc                  <  1:7776608/150‑1 (MQ=255)
                                                                                                                                    aaaaGGCGACAGAGTAATCTGTCGACTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgtgc                  <  2:7619096/150‑1 (MQ=255)
                                                                                                                                     aaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGcgcg                 <  2:1034266/150‑1 (MQ=255)
                                                                                                                                      aaGGCGACAGAGGAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAACGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGt                <  2:8477528/150‑1 (MQ=255)
                                                                                                                                       aGGCGACAGAGTAATCTGTCTCCAATTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTc               <  1:5966639/150‑1 (MQ=255)
                                                                                                                                       aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTc               <  1:641562/150‑1 (MQ=255)
                                                                                                                                       aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCGGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCCTTATGAAACCCTCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTc               <  2:8890028/150‑1 (MQ=255)
                                                                                                                                        ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCg              >  2:6255859/1‑150 (MQ=255)
                                                                                                                                         gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGa             >  2:2103658/1‑150 (MQ=255)
                                                                                                                                         gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGa             >  2:790899/1‑150 (MQ=255)
                                                                                                                                           gACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaa           >  1:5852886/1‑150 (MQ=255)
                                                                                                                                           gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaa           >  2:8182261/1‑150 (MQ=255)
                                                                                                                                           gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaa           >  1:4491874/1‑150 (MQ=255)
                                                                                                                                            acagtctAATCTGTCGCCTTTTTTCATTTCTTGGTTTCTGTTAGATTCCGCCGCAAATCTATACTGTATAAAACGGCTATAAGTATAACATTTATGATAACATCCTATTGAAACTGTCGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaa          <  2:5036128/144‑1 (MQ=255)
                                                                                                                                            aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaa          >  2:9044535/1‑150 (MQ=255)
                                                                                                                                            aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaa          <  2:3167825/150‑1 (MQ=255)
                                                                                                                                             cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaaa         >  1:1164645/1‑150 (MQ=255)
                                                                                                                                             cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaaa         >  1:6122601/1‑150 (MQ=255)
                                                                                                                                               gagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAaca       <  1:5673575/150‑1 (MQ=255)
                                                                                                                                                  tAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACAcgc    >  2:1276511/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGCGA  >  NC_000913/4051114‑4051411

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