Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,241 A→G 6.5% intergenic (‑85/+15) 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,2410AG 6.5% 273.3 / 5.7 91intergenic (‑85/+15)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (45/40);  new base G (3/3);  total (48/43)
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

GTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTC  >  NC_002947/650096‑650387
                                                                                                                                                 |                                                                                                                                                   
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)
            cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                      <  1:214841/149‑1 (MQ=255)
             aGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcga                                                                                                                                        <  8:248274/146‑4 (MQ=255)
             aGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcga                                                                                                                                        >  7:248274/1‑143 (MQ=255)
                   ttCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTgg                                                                                                                               >  2:318674/1‑149 (MQ=255)
                   ttCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTgg                                                                                                                               >  4:326650/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                                                                                              <  6:418737/149‑1 (MQ=255)
                                                   ggATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGcgcg                                                                                              <  1:103133/149‑1 (MQ=255)
                                                        tCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                      >  1:294925/1‑105 (MQ=255)
                                                        tCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                      <  2:294925/105‑1 (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                                                                              >  4:103725/1‑149 (MQ=255)
                                                                   tcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGc                                                                              >  8:204984/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                                                                         <  4:151732/147‑1 (MQ=255)
                                                                          ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATaa                                                                         <  3:103725/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                                            >  6:160011/1‑148 (MQ=255)
                                                                                                      aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                            >  6:110432/1‑148 (MQ=255)
                                                                                                      aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                            >  1:299804/1‑148 (MQ=255)
                                                                                                      aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                            >  2:308240/1‑148 (MQ=255)
                                                                                                      aaaaaaCGTATCTGAAACGTTTCTCTAGAACCCCGGCCCACTAAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGTCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATAAcccc                                            >  2:406300/1‑148 (MQ=255)
                                                                                                          aaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACgg                                       <  2:299804/149‑1 (MQ=255)
                                                                                                           aCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGc                                      <  7:218943/149‑1 (MQ=255)
                                                                                                           aCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGc                                      <  1:308240/149‑1 (MQ=255)
                                                                                                              tATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAccc                                             <  1:48799/139‑1 (MQ=255)
                                                                                                              tATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAccc                                             >  2:48799/1‑139 (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                  >  5:381168/1‑147 (MQ=255)
                                                                                                                                 gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAgt                  <  6:381168/147‑1 (MQ=255)
                                                                                                                                 gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAg                   <  4:325917/146‑1 (MQ=255)
                                                                                                                                 aaaCCCTGGCCTAAGTGGGCCGGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaaaa                          >  4:224151/2‑139 (MQ=255)
                                                                                                                                 aaaCCCTGGCCTAAGTGGGCCGGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaaaa                          <  3:224151/138‑1 (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)
                                                                                                                                                tAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTc  <  1:144207/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTC  >  NC_002947/650096‑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: