Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,262 A→C 7.6% 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.6% 239.9 / 9.7 92noncoding (227/245 nt)csrCCsrC sRNA sequesters CsrA, a carbon flux regulator
Reads supporting (aligned to +/- strand):  ref base A (44/40);  new base C (3/4);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

AAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAAC  >  NC_000913/4051117‑4051405
                                                                                                                                                 |                                                                                                                                               
aaGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACCGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCt                                                                                                                                              <  2:935476/149‑1 (MQ=255)
aaGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCt                                                                                                                                              <  2:47783/149‑1 (MQ=255)
  gACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGt                                                                                                                                            >  2:218060/1‑149 (MQ=255)
   aCAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTc                                                                                                                                           >  1:394220/1‑149 (MQ=255)
     aaaCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                                                                         <  1:1038922/149‑1 (MQ=255)
     aaaCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGc                                                                                                                                         <  2:68716/149‑1 (MQ=255)
       aCGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCt                                                                                                                                       <  1:99193/149‑1 (MQ=255)
             ggtgcagTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGATTACTCTGTCGCCTTTTTTc                                                                                                                                 <  1:304045/143‑1 (MQ=255)
             ggAGACGTAATTAAACGGAAATGGACTCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTAGGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTc                                                                                                                                 <  1:92469/149‑1 (MQ=255)
               tgaCGTAATTACACGGAAAATGAATCAACACGGATTTGTACGGCTTAATGAAAAACTGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGCCCCGTCAAGGGGTATGCGGCCGGCAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCtt                                                                                                                               <  2:781822/148‑1 (MQ=255)
                  cGTAATTAAACGGAATTGGAATCAACACGGATTTTTCAGGCTAAAGGAAAAACAGGTTGTGTTGGCGTCCTGGAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGAGACATCTTTATATGTCGCCTTATTAGTttgg                                                                                                                            >  2:97193/1‑148 (MQ=255)
                  cGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgc                                                                                                                            <  1:630266/149‑1 (MQ=255)
                   gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGTTTATGATTCAGGAAAAAAGCCGACAGAGTAATCTGACGCCTTTTTTCTttgct                                                                                                                           >  1:486855/1‑149 (MQ=255)
                   gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTttgct                                                                                                                           >  1:360286/1‑149 (MQ=255)
                   gTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGAGTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATTTGTCGCCTTTTTTCTttgct                                                                                                                           >  1:241720/1‑149 (MQ=255)
                        tAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTgcttgctt                                                                                                                      <  2:46517/149‑1 (MQ=255)
                         aaaCGGAAATGTAATCAACACTGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGTCCTGCAAGGATTGTAAGACAAGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTTTAGCATTTTTTCTTTGCTTGCttt                                                                                                                     >  2:889264/1‑149 (MQ=255)
                          aaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                                                                    >  2:540003/1‑149 (MQ=255)
                          aaCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTc                                                                                                                    >  2:1071197/1‑149 (MQ=255)
                           aCGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCt                                                                                                                   >  2:980152/1‑149 (MQ=255)
                             ggAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGt                                                                                                                 >  2:1098807/1‑149 (MQ=255)
                             ggAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTATCTTTGCTTGCTTTCTGt                                                                                                                 <  1:579189/149‑1 (MQ=255)
                              gAAATGTAATCAACACGGTTTGTGCCGGCTAAAGGAAAAACCGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                                                                <  2:1079349/149‑1 (MQ=255)
                              gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                                                                <  1:1066460/149‑1 (MQ=255)
                              gAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGtt                                                                                                                <  2:1094995/149‑1 (MQ=255)
                               aaaTGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                               >  1:475204/1‑149 (MQ=255)
                               aaaTGGAATCAACACGGATTGTTACGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTa                                                                                                               <  2:959879/149‑1 (MQ=255)
                                   ggAATCAACACGGATTGTTCCGGCTAAAGTAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAtt                                                                                                           <  2:206020/149‑1 (MQ=255)
                                   ggAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACGCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAtt                                                                                                           >  1:180075/1‑149 (MQ=255)
                                    gAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAtta                                                                                                          >  1:938156/1‑148 (MQ=255)
                                       tCAACACGGATTGTTCCGGCTAAAGTAAAAACAGGGTGTGTTGGCGGCCTTCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTccgc                                                                                                       <  1:1007772/149‑1 (MQ=255)
                                          acacGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTCTAAGACACGTTAACGGTAATGAGTCAGGAAAAAAGGATACAGAGTAATCTGACTCCTGTTTTCTTTGCTTGCTTTCTGTTAGATTCggcggc                                                                                                    >  2:122371/1‑149 (MQ=255)
                                             cGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTTGCGGCATGCAAGGATTGTAAGAACCGTTAAGGGTTATGAGTGAGGAAAAATGGCGACAGAGTAATCTGTCGACTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCaaa                                                                                                 <  2:1019238/149‑1 (MQ=255)
                                               gATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATc                                                                                               >  2:1041946/1‑149 (MQ=255)
                                                    ttCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                                                                          >  1:1066892/1‑149 (MQ=255)
                                                    ttCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                                                                          <  2:11952/149‑1 (MQ=255)
                                                     tCCGGCTCACCGACAAAAAGGTTGTGTCGCCCGCCTGCAAGGACCGTAACACCCGTTACGCGTTAGGAGTCAGCAAAAAACGCCACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACt                                                                                         <  2:765935/149‑1 (MQ=255)
                                                          cTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGACTAATCTGTCGCCTATTTTCTTTGCTTGCTTTCTGTTAGCTTCCGCCGCAAATCTATACTGAATa                                                                                    >  1:1178728/1‑149 (MQ=255)
                                                             aaGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATaaaa                                                                                 <  2:1174436/149‑1 (MQ=255)
                                                              aGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAAAAAAAc                                                                                >  2:468455/1‑149 (MQ=255)
                                                               ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACg                                                                               <  2:288266/149‑1 (MQ=255)
                                                               ggAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACAAGTTAAGGGTTATGATTCAGGAAAAAAGGAGCCAGAGTAATCTGTCGACTTTTTTCTTTGATTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACg                                                                               >  2:569713/1‑149 (MQ=255)
                                                                taaaaaCAGGGTGTGTTGGCGGCATTAAAGGATTGTAAGCCCCTTTAAGGGTTATGAGTCAGGAAAAACGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAATCCGCCGCAAATCTATACTGAATAAAACgg                                                                              <  2:1146538/148‑1 (MQ=255)
                                                                gAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACgg                                                                              >  1:938405/1‑149 (MQ=255)
                                                                        gggTGTGTTGGCGGCCTGCAATGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGa                                                                      <  1:1167478/149‑1 (MQ=255)
                                                                         ggTGGGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGAc                                                                     <  2:866725/149‑1 (MQ=255)
                                                                             tgtTGGCGGCCTGCAAGGATTGTACGACCCGTTAAGGGTTATGAGGCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAc                                                                 <  2:158970/149‑1 (MQ=255)
                                                                              gtTGGCGGCCTGCATGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTACGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAAcc                                                                >  1:371040/1‑149 (MQ=255)
                                                                               ttGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCa                                                               >  1:679848/1‑149 (MQ=255)
                                                                                tGGCGGCCTGAAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAt                                                              <  1:927088/149‑1 (MQ=255)
                                                                                  gcggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGATTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTa                                                            >  1:643520/1‑149 (MQ=255)
                                                                                   cggcCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTAt                                                           <  2:624061/149‑1 (MQ=255)
                                                                                      ccTGCAAGTATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCAATATGaa                                                        >  2:866763/1‑149 (MQ=255)
                                                                                        tGCAAGGATTGTAAGACCCGGTAAGGGTTATGAGTCAGGGAAAAAGGCGACAGAGTATTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGCTTCCGCCGCAAATCTAGACTGAAGAAAACGGCGAAAAGACGACCCATTATGAAAc                                                      >  2:695149/1‑149 (MQ=255)
                                                                                         gCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAcc                                                     <  1:14485/149‑1 (MQ=255)
                                                                                               tttGTAACCCCCGTTAAGGGTTATGAGACAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCatcatc                                               <  2:700325/148‑1 (MQ=255)
                                                                                                  ttAAGACCCGTTAAGGGTTATGAGTCAGTAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTc                                            <  1:569093/148‑1 (MQ=255)
                                                                                                  gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCTCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTc                                            <  1:728640/149‑1 (MQ=255)
                                                                                                  gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg                                                          >  2:682862/1‑135 (MQ=255)
                                                                                                  gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATg                                                          <  1:682862/135‑1 (MQ=255)
                                                                                                  gTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGATTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTc                                            >  1:164548/1‑149 (MQ=255)
                                                                                                    aaGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGTCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCaa                                          <  1:1080036/149‑1 (MQ=255)
                                                                                                     aGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCaaa                                         >  2:583738/1‑149 (MQ=255)
                                                                                                      gACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAAc                                        <  1:394303/149‑1 (MQ=255)
                                                                                                        cccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAATGAATAAAACGGCTAAAAGACCAACCATTATGAAACCATCATCTTCAAACTc                                      >  2:1014289/1‑149 (MQ=255)
                                                                                                         ccGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCCGTTAGATTCCGCCGCAAATCTATACTGAATTAAACGGCTAAATGACGAACCATTATGAAACCATCATCATCCAACTcc                                     >  2:1086380/1‑148 (MQ=255)
                                                                                                           gTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACg                                   <  1:997653/149‑1 (MQ=255)
                                                                                                            ttAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgc                                  >  1:516858/1‑149 (MQ=255)
                                                                                                            ttAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgc                                  <  1:218060/149‑1 (MQ=255)
                                                                                                             tAAGGGTTATGATTCAGGAAAAAAGGAGAGAGAGTAATCTGTCGCATTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACgca                                 <  2:534358/149‑1 (MQ=255)
                                                                                                                  gTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGAATCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCaaa                            <  2:815448/149‑1 (MQ=255)
                                                                                                                     ttgagtgagGAAAAAAGGCGACAGTTTTATCTGTCGGCTTTTTTCTTTGCTTGCGTTTTTGTAGATTCCGCCGCAAATCTATACTTAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGt                         <  1:889264/148‑1 (MQ=255)
                                                                                                                     aTGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGAAAAGGt                         >  2:459504/1‑149 (MQ=255)
                                                                                                                      tGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTc                        >  2:7779/1‑149 (MQ=255)
                                                                                                                       gAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAAGAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCAAGCAGCAAAGTTCa                       >  1:658209/1‑149 (MQ=255)
                                                                                                                        aGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCAt                      <  1:914621/149‑1 (MQ=255)
                                                                                                                        aGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGACTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCAt                      <  1:1041946/149‑1 (MQ=255)
                                                                                                                         gTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATg                     >  2:1130903/1‑149 (MQ=255)
                                                                                                                          tCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGc                    >  1:162960/1‑149 (MQ=255)
                                                                                                                          tCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGc                    >  2:742815/1‑149 (MQ=255)
                                                                                                                           aaagAAAAAAGGCGACAGAGTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCa                   >  1:543127/4‑149 (MQ=255)
                                                                                                                            aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGCCGAACCTTTATGAAACCATCATCTTCAACCTCCCGCAGCAAAGGTCATGCaa                  >  2:230751/1‑149 (MQ=255)
                                                                                                                            aGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACAATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCaa                  >  1:406212/1‑149 (MQ=255)
                                                                                                                             ggAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAAACATCATCTTCAAACTCACGCAGCAAAGGTCATGCaaa                 >  2:446815/1‑149 (MQ=255)
                                                                                                                             ggAAAAAAGGCGACAGAGTAATCTGTCGCCATTTTTCTTTGCTTGCTTTCTGTTAGATTCAGCCGCAAATCTTTACTGAATAAAACGGCTAAAAGAGGAACAATTATGAAACCATCATCTTCAACATAACGCAGCAAAGGTCATGCaaa                 >  1:703924/1‑149 (MQ=255)
                                                                                                                                    aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGt          <  2:1024657/149‑1 (MQ=255)
                                                                                                                                    aGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCGAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAATCGCGt          >  2:584618/1‑149 (MQ=255)
                                                                                                                                      gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTTATAGACGAACCATTATGATACCATCATCTTCAACCTCACGCAGCAAAGCTCATGCAAACGCGCGTCg        >  1:196551/1‑149 (MQ=255)
                                                                                                                                      gCGACAGATTACTCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCg        >  1:50442/1‑149 (MQ=255)
                                                                                                                                      gCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAGCCATTATGAAAGCCTAATCTTCAACCTCACGCCGCAAAGGTCAGGGAACAGCGCTTCg        >  2:185253/1‑149 (MQ=255)
                                                                                                                                           agagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTTGATTCCGCCGCAAATCTATACTGAATAAACCGGCTAAAAGACGAACGATGATGAAACCCTCAGCTTCAACCACCCGCAGCAAAGGTCATGCAAAAGCGCGTCGaacta   >  2:645293/1‑146 (MQ=255)
                                                                                                                                            tagTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAac  <  1:793767/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAAC  >  NC_000913/4051117‑4051405

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: