Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,236 C→T 6.3% intergenic (‑80/+20) 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,2360CT 6.3% 287.2 / 6.1 95intergenic (‑80/+20)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (49/40);  new base T (3/3);  total (52/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 = 5.99e-01

TCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGG  >  NC_002947/650090‑650384
                                                                                                                                                  |                                                                                                                                                     
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                                                                                                                                     <  1:96167/148‑1 (MQ=255)
                 ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                                     >  2:96167/1‑148 (MQ=255)
                 ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtc                                                                                                                                     <  4:147363/148‑1 (MQ=255)
                  cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                    <  5:208314/132‑1 (MQ=255)
                  cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                    <  5:202478/132‑1 (MQ=255)
                  cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                    >  6:208314/1‑132 (MQ=255)
                  cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                    >  6:202478/1‑132 (MQ=255)
                  cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                   <  1:214841/149‑1 (MQ=255)
                   aGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcga                                                                                                                                     <  8:248274/146‑4 (MQ=255)
                   aGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcga                                                                                                                                     >  7:248274/1‑143 (MQ=255)
                         ttCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTgg                                                                                                                            >  4:326650/1‑149 (MQ=255)
                         ttCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTgg                                                                                                                            >  2:318674/1‑149 (MQ=255)
                            cagcagCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCt                                                                                                                         >  7:81127/1‑149 (MQ=255)
                                     gATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGtt                                                                                                                  <  3:326650/147‑1 (MQ=255)
                                      aTCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAgg                                                                                                                                           <  2:143891/121‑1 (MQ=255)
                                      aTCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTAAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAgg                                                                                                                                           >  1:143891/1‑121 (MQ=255)
                                           tttCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGc                                                                                                           >  8:218943/1‑147 (MQ=255)
                                           tttCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCa                                                                                                         <  1:12540/149‑1 (MQ=255)
                                                    tttACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTc                                                                                                  >  2:222072/1‑147 (MQ=255)
                                                       aCGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAgcg                                                                                             >  5:418737/1‑149 (MQ=255)
                                                         ggATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGcgcg                                                                                           <  1:103133/149‑1 (MQ=255)
                                                         ggATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGcgcg                                                                                           <  6:418737/149‑1 (MQ=255)
                                                              tCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                   <  2:294925/105‑1 (MQ=255)
                                                              tCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                   >  1:294925/1‑105 (MQ=255)
                                                               ccATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAg                                                                                                                                                 >  4:209516/1‑90 (MQ=255)
                                                               ccATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAg                                                                                                                                                 <  3:209516/90‑1 (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:151732/1‑148 (MQ=255)
                                                                          caACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGc                                                                           >  3:63873/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                                         >  6:110432/1‑148 (MQ=255)
                                                                                                            aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                         >  1:299804/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                                          >  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                     <  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               <  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)
                                                                                                                                               gCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGaa     >  3:169102/1‑149 (MQ=255)
                                                                                                                                                  cACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACgg  <  5:189607/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGTTGTTCCAGCGCGGACTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGG  >  NC_002947/650090‑650384

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: