Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,051,259 G→T 8.1% 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 8.1% 234.2 / 6.4 91noncoding (224/245 nt)csrCCsrC sRNA sequesters CsrA, a carbon flux regulator
Reads supporting (aligned to +/- strand):  ref base G (44/38);  new base T (3/5);  total (48/43)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAAC  >  NC_000913/4051112‑4051405
                                                                                                                                                   |                                                                                                                                                  
tCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGt                                                                                                                                                   >  2:14485/1‑149 (MQ=255)
     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                                                                                          <  2:11952/149‑1 (MQ=255)
                                                         ttCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATAc                                                                                          >  1:1066892/1‑149 (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)
                                                                                                                                                   |                                                                                                                                                  
TCAGGAAGACAAACGTCCGGAGACGTAATTAAACGGAAATGGAATCAACACGGATTGTTCCGGCTAAAGGAAAAACAGGGTGTGTTGGCGGCCTGCAAGGATTGTAAGACCCGTTAAGGGTTATGAGTCAGGAAAAAAGGCGACAGAGTAATCTGTCGCCTTTTTTCTTTGCTTGCTTTCTGTTAGATTCCGCCGCAAATCTATACTGAATAAAACGGCTAAAAGACGAACCATTATGAAACCATCATCTTCAAACTCACGCAGCAAAGGTCATGCAAAAGCGCGTCGAAAAAC  >  NC_000913/4051112‑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: