Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,239 T→G 6.6% intergenic (‑83/+17) accB ← / ← aroQ‑I acetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947650,2390TG 6.6% 270.9 / 8.3 89intergenic (‑83/+17)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (42/39);  new base G (3/3);  total (47/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTC  >  NC_002947/650091‑650387
                                                                                                                                                    |                                                                                                                                                     
ctcCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                       >  2:103133/1‑149 (MQ=255)
  ccAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                       >  2:12540/1‑147 (MQ=255)
  ccAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTACGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                       >  6:103354/1‑147 (MQ=255)
     gTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGc                                                                                                                                                  <  6:91365/149‑1 (MQ=255)
          tcgATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAggg                                                                                                                                             <  7:335878/149‑1 (MQ=255)
             aTGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGttt                                                                                                                                          <  5:103354/149‑1 (MQ=255)
                ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtca                                                                                                                                       <  1:218684/149‑1 (MQ=255)
                ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                                        <  4:147363/148‑1 (MQ=255)
                ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                                        >  2:96167/1‑148 (MQ=255)
                ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                                        <  1:96167/148‑1 (MQ=255)
                 cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                       <  5:202478/132‑1 (MQ=255)
                 cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                       <  5:208314/132‑1 (MQ=255)
                 cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                       >  6:202478/1‑132 (MQ=255)
                 cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                       >  6:208314/1‑132 (MQ=255)
                 cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                      <  1:214841/149‑1 (MQ=255)
                  aGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcga                                                                                                                                        <  8:248274/146‑4 (MQ=255)
                  aGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcga                                                                                                                                        >  7:248274/1‑143 (MQ=255)
                        ttCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTgg                                                                                                                               >  4:326650/1‑149 (MQ=255)
                        ttCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTgg                                                                                                                               >  2:318674/1‑149 (MQ=255)
                           cagcagCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCt                                                                                                                            >  7:81127/1‑149 (MQ=255)
                                    gATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGtt                                                                                                                     <  3:326650/147‑1 (MQ=255)
                                     aTCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAgg                                                                                                                                              <  2:143891/121‑1 (MQ=255)
                                     aTCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTAAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAgg                                                                                                                                              >  1:143891/1‑121 (MQ=255)
                                          tttCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGc                                                                                                              >  8:218943/1‑147 (MQ=255)
                                          tttCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCa                                                                                                            <  1:12540/149‑1 (MQ=255)
                                                   tttACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTc                                                                                                     >  2:222072/1‑147 (MQ=255)
                                                      aCGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAgcg                                                                                                >  5:418737/1‑149 (MQ=255)
                                                        ggATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGcgcg                                                                                              <  1:103133/149‑1 (MQ=255)
                                                        ggATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGcgcg                                                                                              <  6:418737/149‑1 (MQ=255)
                                                             tCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                      <  2:294925/105‑1 (MQ=255)
                                                             tCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                      >  1:294925/1‑105 (MQ=255)
                                                              ccATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAg                                                                                                                                                    <  3:209516/90‑1 (MQ=255)
                                                              ccATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAg                                                                                                                                                    >  4:209516/1‑90 (MQ=255)
                                                              ccATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTc                                                                                         >  3:325917/1‑148 (MQ=255)
                                                               cATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTc                                                                                         >  7:200036/1‑147 (MQ=255)
                                                                aTTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAg                                                                                      >  8:112517/1‑149 (MQ=255)
                                                                    atcatcAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGc                                                                                   >  6:189607/1‑148 (MQ=255)
                                                                        tcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGc                                                                              >  8:204984/1‑149 (MQ=255)
                                                                        tcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGc                                                                              >  4:103725/1‑149 (MQ=255)
                                                                         caACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGc                                                                              >  3:63873/1‑148 (MQ=255)
                                                                         caACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGc                                                                              >  3:151732/1‑148 (MQ=255)
                                                                          aaCTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGa                                                                            >  6:35419/1‑149 (MQ=255)
                                                                           aCTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcg                                                                                                                                         >  2:174074/1‑86 (MQ=37)
                                                                           aCTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcg                                                                                                                                         <  1:174074/88‑3 (MQ=37)
                                                                               ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATaa                                                                         <  3:103725/147‑1 (MQ=255)
                                                                               ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATaa                                                                         <  4:151732/147‑1 (MQ=255)
                                                                               ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAAcc                                                                       >  8:207861/1‑149 (MQ=255)
                                                                               ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAAcc                                                                       <  3:42326/149‑1 (MQ=255)
                                                                               ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGATGTGCGTTGGCAGCCAGTTGTTACAGCGCGGACTCCAGGGCCAGGCGATAAcc                                                                       >  3:260781/1‑149 (MQ=255)
                                                                               ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGATGTGCGTTGGCAGCCAGTAGTTCCAGCGCGGACTCAAGGGCCAGGCGATaa                                                                         >  5:205847/1‑147 (MQ=255)
                                                                                  aaGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCt                                                                    <  3:315753/149‑1 (MQ=255)
                                                                                      ttCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGgcg                                                                <  7:204984/149‑1 (MQ=255)
                                                                                         ggTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGcc                                                              >  3:117717/1‑148 (MQ=255)
                                                                                          gTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGTCCTCCAGGGCCAGGCGATACCCGCTGgcgac                                                              >  8:286541/1‑145 (MQ=255)
                                                                                             aGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGACAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGa                                                          <  1:318674/149‑1 (MQ=255)
                                                                                                 ggCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCa                                                     <  4:117717/149‑1 (MQ=255)
                                                                                                     aaaTACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGAt                                                  >  7:40756/1‑148 (MQ=255)
                                                                                                        tACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcc                                              <  8:81127/149‑1 (MQ=255)
                                                                                                           aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                            >  2:308240/1‑148 (MQ=255)
                                                                                                           aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                            >  1:299804/1‑148 (MQ=255)
                                                                                                           aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                            >  6:160011/1‑148 (MQ=255)
                                                                                                           aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                            >  6:110432/1‑148 (MQ=255)
                                                                                                           aaaaaaCGTATCTGAAACGTTTCTCTAGAACCCCGGCCCACTAAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGTCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATAAcccc                                            >  2:406300/1‑148 (MQ=255)
                                                                                                               aaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACgg                                       <  2:299804/149‑1 (MQ=255)
                                                                                                                aCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGc                                      <  1:308240/149‑1 (MQ=255)
                                                                                                                aCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGc                                      <  7:218943/149‑1 (MQ=255)
                                                                                                                   tATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAccc                                             >  2:48799/1‑139 (MQ=255)
                                                                                                                   tATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAccc                                             <  1:48799/139‑1 (MQ=255)
                                                                                                                      cTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATc                                >  1:374157/1‑149 (MQ=255)
                                                                                                                       tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGc                                                                     <  2:196180/111‑1 (MQ=255)
                                                                                                                       tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGc                                                                     >  1:196180/1‑111 (MQ=255)
                                                                                                                       tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                               >  6:260949/1‑149 (MQ=255)
                                                                                                                        gAAACGTTTCTCTTGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                               >  1:103261/1‑148 (MQ=255)
                                                                                                                        gAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACgg                                       <  3:103212/140‑1 (MQ=255)
                                                                                                                        gAAACGTTTCTCTCGAACCCCGGCCCACATAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACgg                                       >  4:103212/1‑140 (MQ=255)
                                                                                                                          aaCGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaa                            <  2:57087/149‑1 (MQ=255)
                                                                                                                            cGTTTCTCTCGAACCCCTGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaaaa                          <  6:24714/149‑1 (MQ=255)
                                                                                                                             gTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGCTCACCCCTACGGCAACATCGGAAAAg                         >  1:48325/1‑149 (MQ=255)
                                                                                                                             gTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAg                         >  2:377774/1‑149 (MQ=255)
                                                                                                                              tttCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGt                        <  6:363864/149‑1 (MQ=255)
                                                                                                                              tttCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGt                        <  5:260949/149‑1 (MQ=255)
                                                                                                                                  tctcGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAg                      >  2:144207/1‑147 (MQ=255)
                                                                                                                                    tcGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGa                    <  5:35419/147‑1 (MQ=255)
                                                                                                                                     cGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAg                      <  3:355769/144‑1 (MQ=255)
                                                                                                                                     cGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAg                      >  4:355769/1‑144 (MQ=255)
                                                                                                                                      gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAgt                  <  6:381168/147‑1 (MQ=255)
                                                                                                                                      gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAgt                  >  5:381168/1‑147 (MQ=255)
                                                                                                                                      gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAg                   <  4:325917/146‑1 (MQ=255)
                                                                                                                                      aaaCCCTGGCCTAAGTGGGCCGGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaaaa                          <  3:224151/138‑1 (MQ=255)
                                                                                                                                      aaaCCCTGGCCTAAGTGGGCCGGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaaaa                          >  4:224151/2‑139 (MQ=255)
                                                                                                                                         ccccGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGa             <  5:160011/149‑1 (MQ=255)
                                                                                                                                         ccccGGCCCACATAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGATGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGAGGTGa             >  1:20413/1‑149 (MQ=255)
                                                                                                                                           ccGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACg           <  5:110432/149‑1 (MQ=255)
                                                                                                                                              gCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGaa        >  3:169102/1‑149 (MQ=255)
                                                                                                                                                 cACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACgg     <  5:189607/149‑1 (MQ=255)
                                                                                                                                                   cTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGtt   >  7:42668/1‑149 (MQ=255)
                                                                                                                                                    ttAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTc  >  7:336140/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTC  >  NC_002947/650091‑650387

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