Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,246 A→G 6.3% intergenic (‑90/+10) 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,2460AG 6.3% 281.0 / 5.6 94intergenic (‑90/+10)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (46/42);  new base G (3/3);  total (49/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

TCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGTTTG  >  NC_002947/650101‑650394
                                                                                                                                                 |                                                                                                                                                     
tcgATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAggg                                                                                                                                                    <  7:335878/149‑1 (MQ=255)
   aTGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGttt                                                                                                                                                 <  5:103354/149‑1 (MQ=255)
      ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtca                                                                                                                                              <  1:218684/149‑1 (MQ=255)
      ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                                               <  1:96167/148‑1 (MQ=255)
      ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                                               <  4:147363/148‑1 (MQ=255)
      ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                                               >  2:96167/1‑148 (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                                                                                                                                             <  2:294925/105‑1 (MQ=255)
                                                   tCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                             >  1:294925/1‑105 (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                                                   >  1:299804/1‑148 (MQ=255)
                                                                                                 aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                                   >  6:110432/1‑148 (MQ=255)
                                                                                                 aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                                   >  6:160011/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                                                                            >  1:196180/1‑111 (MQ=255)
                                                                                                             tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGc                                                                            <  2:196180/111‑1 (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                               <  5:260949/149‑1 (MQ=255)
                                                                                                                    tttCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGt                               <  6:363864/149‑1 (MQ=255)
                                                                                                                        tctcGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAg                             >  2:144207/1‑147 (MQ=255)
                                                                                                                          tcGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGa                           <  5:35419/147‑1 (MQ=255)
                                                                                                                           cGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAg                             >  4:355769/1‑144 (MQ=255)
                                                                                                                           cGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAg                             <  3:355769/144‑1 (MQ=255)
                                                                                                                            gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAgt                         >  5:381168/1‑147 (MQ=255)
                                                                                                                            gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAgt                         <  6:381168/147‑1 (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)
                                                                                                                                           tAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTc         <  1:144207/148‑1 (MQ=255)
                                                                                                                                              gCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGc                                             >  7:144356/1‑109 (MQ=255)
                                                                                                                                              gCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGc                                             <  8:144356/109‑1 (MQ=255)
                                                                                                                                              gCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGt     <  4:63873/149‑1 (MQ=255)
                                                                                                                                              gCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGAAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGt     <  1:406300/149‑1 (MQ=255)
                                                                                                                                              gCCAGGGTTTTCATCATTTGGGCTGTGCATTAGCAGCCAATTGTTCCAGCGCGGACTACAGGGCAAAGCAATAACCGCTTGAGCCCAGACCGCAGATAAACCCTAAGGCAACATCGGAAAAGTAGGAGTGGTGa                    >  7:258456/1‑134 (MQ=255)
                                                                                                                                                cAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGttt   <  8:200036/149‑1 (MQ=255)
                                                                                                                                                 aGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGTTtg  <  1:377774/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGTTTG  >  NC_002947/650101‑650394

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: