Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,259 G→T 5.8% noncoding (224/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,2590GT 5.8% 226.0 / 2.6 86noncoding (224/245 nt)csrCCsrC sRNA sequesters CsrA, a carbon flux regulator
Reads supporting (aligned to +/- strand):  ref base G (48/33);  new base T (2/3);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.41e-01

TCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAA  >  NC_000913/4051112‑4051404
                                                                                                                                                   |                                                                                                                                                 
tCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAACTGTTATTAGTCAGGAAAAAAGGCGACAGAGt                                                                                                                                                  >  2:1006953/1‑149 (MQ=255)
 cgggAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTa                                                                                                                                                 <  1:299443/147‑1 (MQ=255)
  aGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTaa                                                                                                                                                >  1:713984/1‑149 (MQ=255)
   ggAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTTAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                                                               <  1:625115/149‑1 (MQ=255)
    gAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGACGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGAAAGAGTAATc                                                                                                                                              >  1:886998/1‑149 (MQ=255)
      aGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACATGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACACGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTg                                                                                                                                            >  2:716486/1‑149 (MQ=255)
          aaaCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                                                                        >  1:440650/1‑149 (MQ=255)
          aaaCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                                                                        >  1:436574/1‑149 (MQ=255)
           aaCGTCCGCCGACGTAATTAAACGGAAATCTCATCAACACGGATTGTTCCGGCTAAAGGAAAACCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGcc                                                                                                                                       <  2:102711/149‑1 (MQ=255)
              gTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTTTAAGACACGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCatt                                                                                                                                    >  1:883761/1‑149 (MQ=255)
                ccGCAGACGTAATTAACCGGGAATGGCATCAACACGGATTGTTCCGGCTAAAGGAAAACCCGGGTGAGTTGGCGGCCTGCACGGATTGTAAGACCCGATAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCttttt                                                                                                                                  <  2:217550/149‑1 (MQ=255)
                  ggAGACGTAATTAACCGGAAATGGAATCAACACCTATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTc                                                                                                                                <  1:103514/149‑1 (MQ=39)
                   gagaCGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGAATGTAAGACCCGTTAAGGGTTACGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCt                                                                                                                               >  2:240635/1‑149 (MQ=255)
                      aCGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttg                                                                                                                            <  1:9533/149‑1 (MQ=255)
                        gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGTCGACAGAGTAATCTGTCGACTTTTTTCTttgct                                                                                                                          >  1:518915/1‑149 (MQ=255)
                           aTTAAACGGAAATGGAATCAACACGGATTGTTCCGGCAAAAGG‑AAAACAGGGTGTGTTGGCGGCCTGCAAGGAATGTAAGACCCGTTAAGGGTTAAGAGTAAGGAAAAAAGGCGTCAGAGTAATCTCTCGACTTTATACTTtgctagct                                                                                                                      >  2:872424/1‑149 (MQ=255)
                            ttAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTtgcttgct                                                                                                                      <  2:63389/149‑1 (MQ=255)
                            ttAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCCCCTTTTTTCTTtgcttgct                                                                                                                      >  1:865875/1‑149 (MQ=255)
                              aaaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCttt                                                                                                                    <  2:89486/149‑1 (MQ=255)
                               accGGAAATGGATTCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                                                                   <  2:161098/147‑1 (MQ=255)
                                 cGGAAATGGAATCAACACGGATTTTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATATGTCGCCTTTTTTCTTTGCTTGCTTTCTg                                                                                                                 <  2:426211/149‑1 (MQ=255)
                                   gAAATGGTATCAACACGGATTGTTCCGTCTAAAGGAAAAAGAGGGGGTGTTGGCGGCCTGCATGGTTTTTAAGACCCGTTAAGTGTTATGAGTCAGGAAAAACGGCGACAGAGTAATCTGTCGCGTTTTTTCTTTGCTTGCTTTCTGtt                                                                                                               <  1:411194/149‑1 (MQ=255)
                                   gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                                                               <  2:1029579/149‑1 (MQ=255)
                                   gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                                                                <  2:736343/148‑1 (MQ=255)
                                   gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                                                                >  1:736343/1‑148 (MQ=255)
                                    aaaTGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                              >  1:45023/1‑149 (MQ=255)
                                      aTGGAATCAACACGGATTGTTCCGGCTAAAGGAAATACAGGGGGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                                                                            <  1:1016175/149‑1 (MQ=255)
                                            tCAACACGGATTGTTCCGGCTAAAGGGAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                                                                      >  2:394206/1‑149 (MQ=255)
                                            tCAACACGGATTGTTCCGGCTAAAGGAAAAACAGTGTGGGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                                                                      <  2:253478/149‑1 (MQ=255)
                                               acacGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                                                                   >  1:450191/1‑149 (MQ=255)
                                               acacGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                                                                   >  1:299505/1‑149 (MQ=255)
                                                  cGGATTGTTCCGGCTAAAGGAAAAACAGGGGGTGTTGGCGGCCTGCAAGGATTGTAAGAGCCGTTAAGGGTTATGAGTCCGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaat                                                                                                >  2:660135/1‑148 (MQ=255)
                                                    gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCAGTTAGATTCCGCCGCAAATc                                                                                              >  1:126136/1‑149 (MQ=255)
                                                        gTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtata                                                                                          >  1:170858/1‑149 (MQ=255)
                                                         ttCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTATGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                                                                         >  2:39910/1‑149 (MQ=255)
                                                         ttCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTACTGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                                                                         <  1:133619/149‑1 (MQ=255)
                                                          tCCGGCTTAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                                                                        <  1:155699/149‑1 (MQ=255)
                                                             ggcgcaaGGAATCACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGaa                                                                                     <  2:237146/144‑1 (MQ=255)
                                                              gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAAt                                                                                    >  2:879275/1‑149 (MQ=255)
                                                                  aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaca                                                                                >  1:896116/1‑147 (MQ=255)
                                                                    ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTTAGGGTGATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACg                                                                              >  2:176232/1‑149 (MQ=255)
                                                                     gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                                                             <  1:193021/149‑1 (MQ=255)
                                                                     gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGATTGCGTTCTTTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                                                             >  2:846788/1‑149 (MQ=255)
                                                                        aaaCAGGGTGTTTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTGTTTTCTTTTCTTGCTTTCTGTTATATTCCGCCGCAAATCGATACTGAATAAAACGCCtt                                                                          >  2:902222/1‑148 (MQ=255)
                                                                        aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                                                                          >  1:740117/1‑149 (MQ=255)
                                                                            aGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACACGTTAAGGGTTATGAGTCAGGAAAAAAGGAGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAg                                                                      <  2:822135/149‑1 (MQ=255)
                                                                                  tgtTGGCGGCCTGCAAGGATTGTAAGCCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTTGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                                                                >  2:111375/1‑149 (MQ=255)
                                                                                   gtTGGCGGCCTGCAAGGATTGTAAGACACGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGGGATCTGTCGCCTGTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCGGCAAATGTATACTGAATAAAACGGCTAAAAGTCGAAcc                                                               >  1:45191/1‑149 (MQ=255)
                                                                                    ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                                                              <  1:510615/149‑1 (MQ=255)
                                                                                       gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                                                           <  1:935135/149‑1 (MQ=255)
                                                                                        cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt                                                          >  1:891514/1‑149 (MQ=255)
                                                                                        cggcCTGCAAGGATGGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAATTCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt                                                          <  1:722409/149‑1 (MQ=255)
                                                                                         ggcCGGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGTCGACAGAGTAATCTGTCGGCTGTTTTCTTTGCTGGCTTGGTGTTAGATACCGCGCCAAATCAATAATGATTAGAACGGCTACAAGACGAGcccggagg                                                         >  2:1024342/1‑143 (MQ=255)
                                                                                              gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAGCGGCTAAAAGACGAACCATTATGAAAcc                                                    <  2:569237/149‑1 (MQ=255)
                                                                                              gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAcc                                                    <  2:479657/149‑1 (MQ=255)
                                                                                              gCAAGGATTGCAAGACCAGTAAAGGGTTATGAGTCATGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAATCCGCCGCAAATCTATACTGAATAAAACTGCTAAAAGACGAACCATTATGAAAcc                                                    <  2:684332/149‑1 (MQ=255)
                                                                                               cAAGGATTGTAAGACCCGTTAAGGGTTATGACTCAGGAAAAAAGGCGACAGAGTAATCTGTAGCATTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGACTCAAATCTATACTGAATAAAACGGATAAAAGACGAACCATTATGAAACca                                                   >  1:279423/1‑149 (MQ=255)
                                                                                                          aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAACATCATCTTCaaa                                        >  2:566013/1‑149 (MQ=255)
                                                                                                            aCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTATTTGGTTGCTTTATGTTAGATTCCGCCGCAAAAGTATACTGAATAAAACGGCTAAATGACGAACCATTATGAAACCATAATCTTCAAACt                                      >  1:932106/1‑149 (MQ=255)
                                                                                                             cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTc                                     <  1:879275/149‑1 (MQ=255)
                                                                                                                gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAATCCATCATCTTCAAACTCACg                                  <  2:432044/149‑1 (MQ=255)
                                                                                                                   aaGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgcag                               >  2:335227/1‑149 (MQ=255)
                                                                                                                    aGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgcagc                              <  1:206083/149‑1 (MQ=255)
                                                                                                                     gggTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgcagca                             >  2:847122/1‑149 (MQ=255)
                                                                                                                         tATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAACACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAgg                         >  2:129055/1‑149 (MQ=255)
                                                                                                                          aTGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGt                        <  1:716486/149‑1 (MQ=255)
                                                                                                                           tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTc                       >  1:566327/1‑149 (MQ=255)
                                                                                                                           tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTc                       >  1:512928/1‑149 (MQ=255)
                                                                                                                                 aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatc                                                 <  2:377670/117‑1 (MQ=255)
                                                                                                                                 aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcatc                                                 >  1:377670/1‑117 (MQ=255)
                                                                                                                                 aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCCACCTCACGCAGCAAAGGTCATGCaa                 >  2:368001/1‑149 (MQ=255)
                                                                                                                                 aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaa                 <  1:846788/149‑1 (MQ=255)
                                                                                                                                   gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaaa               >  2:548876/1‑149 (MQ=255)
                                                                                                                                     aaaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgc             <  2:1014180/149‑1 (MQ=255)
                                                                                                                                        aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAACAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGcgcg          <  2:440650/149‑1 (MQ=255)
                                                                                                                                        aaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGcgcg          <  2:436574/149‑1 (MQ=255)
                                                                                                                                         aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATTCTGAATAAAACGGCTAAAAGACGACCCTTTATGAACCCTTCATCTTCAAACTCACGCAGCAAAGGTCATGCAATAGCGCGt         >  1:935116/1‑149 (MQ=255)
                                                                                                                                          ggCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTc        >  2:331720/1‑149 (MQ=255)
                                                                                                                                          ggCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAcgcgcggc        >  2:971280/1‑147 (MQ=255)
                                                                                                                                           gagACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCg       <  1:240635/147‑1 (MQ=255)
                                                                                                                                           gCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGGTTTCTGTTAGATTCCGCCGCAAATCGATACTGAAGAAAACGGCTAAAAGACGAACCATTATGCAACCATCATCTTCAAACTCAGGCAGCAAAGGTCATGCAAAAGCGCGTCg       >  1:841982/1‑149 (MQ=255)
                                                                                                                                             gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCTTCTTCAACCGCACGCAGCAAAGGTCATGCAAAAGCGCGGCGaa     >  1:70594/1‑149 (MQ=255)
                                                                                                                                             gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaa     >  1:973020/1‑149 (MQ=255)
                                                                                                                                              aCAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGAAAAAGAGAGTCGaaa    >  2:120657/1‑149 (MQ=255)
                                                                                                                                              aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaa    >  1:514527/1‑149 (MQ=255)
                                                                                                                                              aCAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaa    >  1:1011246/1‑149 (MQ=255)
                                                                                                                                                agagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaaa  <  1:394206/149‑1 (MQ=255)
                                                                                                                                                agagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaaa  <  2:45023/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAA  >  NC_000913/4051112‑4051404

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑