Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,231 C→T 6.4% intergenic (‑75/+25) 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,2310CT 6.4% 279.6 / 6.1 94intergenic (‑75/+25)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (49/39);  new base T (3/3);  total (52/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 = 6.68e-01

TTGATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACG  >  NC_002947/650086‑650378
                                                                                                                                                 |                                                                                                                                                    
ttGATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGc                                                                                                                                                   >  4:42326/1‑149 (MQ=255)
    tctcCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATAATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTAACTGAAACGTTTCTCTCGAACCCCGGCCCAc                                                                                                                                               >  1:57087/1‑149 (MQ=255)
     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                                                                                                                               <  1:96167/148‑1 (MQ=255)
                     ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                               >  2:96167/1‑148 (MQ=255)
                      cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                              <  5:208314/132‑1 (MQ=255)
                      cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                              <  5:202478/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                                                                                                                               >  7:248274/1‑143 (MQ=255)
                       aGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcga                                                                                                                               <  8:248274/146‑4 (MQ=255)
                             ttCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTgg                                                                                                                      >  2:318674/1‑149 (MQ=255)
                             ttCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTgg                                                                                                                      >  4:326650/1‑149 (MQ=255)
                                cagcagCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCt                                                                                                                   >  7:81127/1‑149 (MQ=255)
                                   cagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCg                                                                                                                                                     <  8:83942/112‑1 (MQ=255)
                                   cagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCg                                                                                                                                                     >  7:83942/1‑112 (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)
                                                                   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                                                              <  3:42326/149‑1 (MQ=255)
                                                                                    ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAAcc                                                              >  8:207861/1‑149 (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                             <  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                                                            >  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               <  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                 >  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)
                                                                                                                                                 |                                                                                                                                                    
TTGATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACG  >  NC_002947/650086‑650378

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: