Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 5.8% 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.8% 230.6 / 2.4 85noncoding (227/245 nt)csrCCsrC sRNA sequesters CsrA, a carbon flux regulator
Reads supporting (aligned to +/- strand):  ref base A (46/34);  new base C (2/3);  total (48/37)
Fisher's exact test for biased strand distribution p-value = 6.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

AGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGC  >  NC_000913/4051114‑4051409
                                                                                                                                                    |                                                                                                                                                   
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)
                         aTTAAACGGAAATGGAATCAACACGGATTGTTCCGGCAAAAGGAAAACAGGGTGTGTTGGCGGCCTGCAAGGAATGTAAGACCCGTTAAGGGTTAAGAGTAAGGAAAAAAGGCGTCAGAGTAATCTCTCGACTTTATACTTtgctagct                                                                                                                           >  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                                                                                                                     >  1:736343/1‑148 (MQ=255)
                                 gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                                                                     <  2:736343/148‑1 (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:512928/1‑149 (MQ=255)
                                                                                                                         tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTc                            >  1:566327/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       <  2:45023/149‑1 (MQ=255)
                                                                                                                                              agagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaaa       <  1:394206/149‑1 (MQ=255)
                                                                                                                                                   aaTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACAcgc  <  1:566013/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGC  >  NC_000913/4051114‑4051409

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: