Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 7.2% 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 7.2% 187.5 / 3.5 79noncoding (227/245 nt)csrCCsrC sRNA sequesters CsrA, a carbon flux regulator
Reads supporting (aligned to +/- strand):  ref base A (45/28);  new base C (4/2);  total (49/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

AGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACACGC  >  NC_000913/4051114‑4051409
                                                                                                                                                    |                                                                                                                                                   
aGGAAGACAAACGTCCGGAGACGTAATTAAACTGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTaa                                                                                                                                                     >  2:245749/1‑149 (MQ=255)
 ggAAGGCAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAAt                                                                                                                                                    <  1:405426/149‑1 (MQ=255)
  gAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATc                                                                                                                                                   >  2:207090/1‑149 (MQ=255)
        aaaCGTCCGGAGACATAATTAAACGGAAATGGAATCAACACTGATTGTTCCGGCTAAAGCAAAAACAGGGTGTGTTGGCTGCCTGCAAGGATTGTAAGACCCGTTAAGTGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGc                                                                                                                                             >  1:654337/1‑149 (MQ=255)
            gTCCGGAGACGTAATTAACCGGAAATGGAATCACCACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTCGCGGCCTGAAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCttt                                                                                                                                         >  2:846145/1‑149 (MQ=255)
                     cGTAATTAAACGGAAATGTAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGTGTGTGTTGGCTGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgc                                                                                                                                <  2:940624/149‑1 (MQ=255)
                     cGTAATTAAACGGAAATGGAATCCACACGGATTGGTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgc                                                                                                                                <  1:825504/149‑1 (MQ=255)
                      gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                                                                               >  1:825165/1‑149 (MQ=255)
                      gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGGGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCCCCTTTTTTCTttgct                                                                                                                               >  2:414726/1‑149 (MQ=255)
                         aTTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgcttgc                                                                                                                            <  1:26731/149‑1 (MQ=255)
                           tAAACGGAAAAGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATAGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGCCAGAGTAATCTGTCGCCATTTTTCTTTgcttgctt                                                                                                                          >  2:305883/1‑149 (MQ=255)
                             aaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGTGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATATGTCGCCATTTTTCTTTGCTTGCATTc                                                                                                                        >  1:45306/1‑149 (MQ=255)
                               cGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTg                                                                                                                      >  1:915383/1‑149 (MQ=255)
                                 gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAATACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGGCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                                                                    <  1:158615/149‑1 (MQ=255)
                                  acaTGGAATCAACACGGCCTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCGGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCGTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                                   <  2:989932/147‑1 (MQ=255)
                                  aaaTGGAATCAACACGGATTGTTCCGGCTAACGGAAAAACAGGGGGTGTTGGCGGCCTGCAAGGATTGTAAGCCCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                                   <  1:228716/149‑1 (MQ=255)
                                  aaaTGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                                   >  1:117127/1‑149 (MQ=255)
                                    aTGGAATCAACACGGATTGTTCCGGCTAAAGTAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                                                                                 <  2:636468/149‑1 (MQ=255)
                                    aTGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGa                                                                                                                 <  1:619040/149‑1 (MQ=255)
                                     tGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGTTTGTAAGACCCCTTAAGGGTTATGAGTCAGGAAAAAAGGAGACAGAGTAATCTGTCTCCTTTTTTCTTTGATAGCTTTCTGTTTGAt                                                                                                                >  1:562695/1‑149 (MQ=255)
                                     tGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAt                                                                                                                >  2:837164/1‑149 (MQ=255)
                                      ggAATCAACACGGATTGTTCCGGCTAAAGGAAAACCAGGGGGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAtt                                                                                                               <  2:210332/149‑1 (MQ=255)
                                       gAATCAACAGGGATTTTTCCGGCTAAAGGAAAAACTGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTc                                                                                                              <  1:497811/149‑1 (MQ=255)
                                            aaCACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGGAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgccg                                                                                                         >  1:108774/1‑149 (MQ=255)
                                             acacGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCcgccgc                                                                                                        >  2:288665/1‑149 (MQ=255)
                                                       ttCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCGTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCtataa                                                                                              >  2:784941/1‑148 (MQ=255)
                                                         ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                                                                            <  1:238143/149‑1 (MQ=255)
                                                         ccGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTg                                                                                            <  1:24429/149‑1 (MQ=255)
                                                          cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGa                                                                                           >  2:783952/1‑149 (MQ=255)
                                                          cGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTGTCTGTGAGCTTCCGCAGCAAATCTATACTGa                                                                                           >  1:915858/1‑149 (MQ=255)
                                                           ggCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGaa                                                                                          <  1:261757/149‑1 (MQ=255)
                                                            gCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAAt                                                                                         >  1:547078/1‑149 (MQ=255)
                                                              tAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaa                                                                                       >  2:740058/1‑149 (MQ=255)
                                                                aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTAGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaaa                                                                                     >  2:678340/1‑149 (MQ=255)
                                                                 aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                                                                    >  2:352804/1‑149 (MQ=255)
                                                                 aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAAc                                                                                    >  2:226437/1‑149 (MQ=255)
                                                                  ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCTGGATAATAGGCGACAGAGTAATGTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACg                                                                                   >  1:71908/1‑149 (MQ=255)
                                                                  ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGAGACAGAGTAATCTGTCGCATTTTTTCTTGGAATCCTAACTGTAACATACCGCCGCAAAACTATAATGAATAAAACg                                                                                   >  2:632719/1‑149 (MQ=255)
                                                                   gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                                                                  >  1:826481/1‑149 (MQ=255)
                                                                   gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGACGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                                                                  >  2:460213/1‑149 (MQ=255)
                                                                   gAAAAACAGGGGGGGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                                                                  <  2:646057/149‑1 (MQ=255)
                                                                     aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCt                                                                                >  1:184068/1‑149 (MQ=255)
                                                                     aaaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCt                                                                                >  1:145657/1‑149 (MQ=255)
                                                                      aaaCAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTa                                                                               >  1:23154/1‑149 (MQ=255)
                                                                       aaCAGGGTGTGTTGGCGGCCTGCAAGGATAGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCCCAAATCTATACTGAATAAAACGGCTaa                                                                              >  2:955498/1‑149 (MQ=255)
                                                                                tgtTTGCGGCCTGCAAGGATTGAAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGAGTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGACGCAAATCTATAGTGAATAAAACGGCTAAAATCCGaaa                                                                     >  2:993103/1‑148 (MQ=255)
                                                                                tgtTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGTCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                                                                     <  1:173965/149‑1 (MQ=255)
                                                                                ggtTGGCGGCCTTCAAGGATTGTAAGACCCTTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAACAGACGAAc                                                                     <  2:711843/148‑1 (MQ=255)
                                                                                      cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt                                                               <  2:725762/149‑1 (MQ=255)
                                                                                      cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCGTTGCTTGCTTGCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAGACGAACCATTATg                                                              <  2:468106/149‑1 (MQ=255)
                                                                                      cggacagCAATGAATGAAAGTCCCGGTAACGGTTATTAGTCAGGAAAAAAGGCGAGACAGTTATCTGTCGCCTTTTTTCTTTGCGTGCTTTCTGGCAGAATCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt                                                               <  2:642034/143‑1 (MQ=255)
                                                                                          cTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTACGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGaaa                                                           >  1:358246/1‑149 (MQ=255)
                                                                                             cAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACca                                                        >  2:29240/1‑149 (MQ=255)
                                                                                              aaGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACcat                                                       >  1:41884/1‑149 (MQ=255)
                                                                                                  aTTGTAAGACCCTGTAAGGGTTATTAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCatcatc                                                   <  1:353960/149‑1 (MQ=255)
                                                                                                   ttGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCt                                                  >  1:68976/1‑149 (MQ=255)
                                                                                                        ataCCCGTTAATGGTTATGAGTCAGGACAAATGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGATTGCATTCTGTTAGATTCCGCCGCAAATCAATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCaaa                                             <  1:696702/147‑1 (MQ=255)
                                                                                                        aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCaaa                                             <  1:245749/149‑1 (MQ=255)
                                                                                                         gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGATTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTTAAAGACGCACCATTATGAAACCATCATCTGAAAAc                                            >  2:870779/1‑149 (MQ=255)
                                                                                                            ccGTTAAGGGTTATGAGTCAGGTAAAAAGGCGACAGAGTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCa                                         <  2:840100/149‑1 (MQ=255)
                                                                                                            ccGTTAAGGGGTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCa                                         >  1:81616/1‑149 (MQ=255)
                                                                                                             gggtatgggTTATGAGTCTGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCTAACTCAc                                        <  2:562695/143‑1 (MQ=255)
                                                                                                              gTTAATGGTTTTGTGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCCTCATCTTCAAACTCACg                                       <  1:777094/149‑1 (MQ=255)
                                                                                                                 acgggTTATGAGTCAGGAAAAACGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgcag                                    <  2:714942/147‑1 (MQ=255)
                                                                                                                    ggTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTTAATAAACCGGCTAAAAGACGAACCTTTATGAAACCATCATCTTCAAACTCACGCAGCaa                                 >  2:846298/1‑149 (MQ=255)
                                                                                                                         tGAGTCAGGAAAATAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGATCCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTc                            <  1:140298/149‑1 (MQ=255)
                                                                                                                             tCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGc                        <  2:915383/149‑1 (MQ=255)
                                                                                                                              aaagAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCa                       >  1:316088/4‑149 (MQ=255)
                                                                                                                                ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa                     >  2:370168/1‑149 (MQ=255)
                                                                                                                                 gAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAACCCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaaa                    >  1:480766/1‑149 (MQ=255)
                                                                                                                                   aaaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgc                  >  1:324725/1‑149 (MQ=255)
                                                                                                                                    aaaaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgcg                 >  2:63591/1‑149 (MQ=255)
                                                                                                                                    aaaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgcg                 >  2:485645/1‑149 (MQ=255)
                                                                                                                                     aaaGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACAATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAgcgc                >  2:216917/1‑149 (MQ=255)
                                                                                                                                      aaGGCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTAACTGTTAGCTTCCGCCGCAAATCTATACTGACTAAAACGGCGAAAAGACGAACCATTATGAACCCGGCATCTTCACACTCACGCAGCAACGGTCATGCAAAAGcgcg               >  2:65920/1‑149 (MQ=255)
                                                                                                                                           gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAATCGGCTAAAAGACGAATCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaa          <  2:45306/149‑1 (MQ=255)
                                                                                                                                           gACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGCACCATTATGAAACCATCTTCTTCAAACTCCCGCAGCAAAGGTCATCCAAAAGCGCGTCGaa          >  2:922800/1‑149 (MQ=255)
                                                                                                                                             cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGTTTCCGCCGCAAATCTATACTGAATTAAACGGCTAAAAGACGAACCATTATGAACCCATCATCGTCAAACTCACGCAGCAAAGGTCATGCCAAAGCGCGTCGaaaa        >  2:103755/1‑149 (MQ=255)
                                                                                                                                             cAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAATGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGaaaa        >  2:515126/1‑149 (MQ=255)
                                                                                                                                                agTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAaca     >  1:69586/1‑149 (MQ=255)
                                                                                                                                                   aaTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATTAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAACGCGCGTCGCAAAACAcgc  >  1:499263/1‑149 (MQ=255)
                                                                                                                                                   aaTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAACAcgc  <  1:298383/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: