Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,231 C→T 9.2% 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 9.2% 237.2 / 11.3 88intergenic (‑75/+25)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (42/37);  new base T (3/5);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.21e-01

CTTGATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGAC  >  NC_002947/650085‑650377
                                                                                                                                                   |                                                                                                                                                       
cttGATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAAAACAAAACACGTATCTGAAACGTTTCTCTCGAACCCCg                                                                                                                                                        <  2:120684/148‑1 (MQ=255)
cttGATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTACCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCgg                                                                                                                                                       >  1:148620/1‑149 (MQ=255)
   gATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGcc                                                                                                                                                     >  2:34656/1‑148 (MQ=255)
   gATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGcc                                                                                                                                                     >  1:16086/1‑148 (MQ=255)
    aTCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCa                                                                                                                                                   >  1:75900/1‑149 (MQ=255)
     tctcCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCAc                                                                                                                                                  >  5:296338/1‑149 (MQ=255)
       tcCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                 >  3:190800/1‑148 (MQ=255)
         cAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACACCCAAAGGTTCGGTCAGGGGGCTAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGCAACACTTAg                                                                                                                                              >  2:331587/1‑149 (MQ=255)
               gtcgATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAgg                                                                                                                                        <  8:129230/149‑1 (MQ=255)
               gtcgATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAgg                                                                                                                                        <  5:229626/149‑1 (MQ=255)
                tcgATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAggg                                                                                                                                       <  2:16086/149‑1 (MQ=255)
                    tGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttc                                                                                                                                    <  4:123940/148‑2 (MQ=255)
                      ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcga                                                                                                                                  <  5:1121/148‑4 (MQ=255)
                      ccAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtca                                                                                                                                 <  7:21692/149‑1 (MQ=255)
                       cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcgaga                                                                                                                                <  8:154363/149‑6 (MQ=255)
                       cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtcat                                                                                                                                >  2:332210/1‑149 (MQ=255)
                       cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTtca                                                                                                                                 >  7:171740/1‑148 (MQ=255)
                       cAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTAAGGCCAGGGTTTtca                                                                                                                                 >  5:180206/1‑148 (MQ=255)
                           ctctTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCAtt                                                                                                                            <  1:193081/149‑1 (MQ=255)
                                  agcagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTTAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtg                                                                                                                     >  7:312387/1‑149 (MQ=255)
                                   gcagcTCGATGAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAACGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtgt                                                                                                                    <  4:151325/149‑1 (MQ=255)
                                   gcagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtgt                                                                                                                    >  2:248534/1‑149 (MQ=255)
                                   gcagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtgt                                                                                                                    <  1:149075/149‑1 (MQ=255)
                                   gcagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtgt                                                                                                                    <  5:83857/149‑1 (MQ=255)
                                          gATCAGTTTCTTGACTTTACGGATATCCTTTAATCATCAACTCCCAAAGGTTCTGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCGACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTg                                                                                                              <  1:331587/148‑1 (MQ=255)
                                          gATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTACTCGAACCCTGGCCTAAGTGGGCCGGGGttcgaga                                                                                                                                <  1:246931/131‑6 (MQ=255)
                                          gATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcgaga                                                                                                                                >  2:246931/1‑125 (MQ=255)
                                          gATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGtt                                                                                                               <  1:47436/147‑1 (MQ=255)
                                          gATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTAAGGTAAGGGGGCGAAATACAAAAAACGTAAAAGAAACGTTTATATAGAAACCCGGCCCAAATAGGCCAAGGTTTTAAAAATTTTGGCTTTTCGTTgg                                                                                                             >  1:282799/1‑149 (MQ=255)
                                           aTCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcgag                                                                                                                                 <  4:315216/128‑5 (MQ=255)
                                           aTCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcgag                                                                                                                                 >  3:315216/1‑124 (MQ=255)
                                              aGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGc                                                                                                         <  4:190800/149‑1 (MQ=255)
                                               gTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGcc                                                                                                        <  3:192883/149‑1 (MQ=255)
                                               gTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGc                                                                                                         >  4:48847/1‑148 (MQ=255)
                                               gTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTACGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGc                                                                                                         >  1:83394/1‑148 (MQ=255)
                                                      gaagttACGTATATCCCTTACTCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAgt‑‑tgt                                                                                                <  7:260987/144‑1 (MQ=255)
                                                       aCTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTc                                                                                              <  2:318580/149‑1 (MQ=255)
                                                                    ccATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑Tc                                                                               >  7:335870/1‑148 (MQ=255)
                                                                     cATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAgt‑‑tgtt                                                                                               <  1:73355/134‑1 (MQ=255)
                                                                     cATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAgt‑‑tgtt                                                                                               >  2:73355/1‑134 (MQ=255)
                                                                     cATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑Tc                                                                               >  1:360256/1‑147 (MQ=255)
                                                                      aTTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAg                                                                            >  6:142916/1‑149 (MQ=255)
                                                                       ttAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGttt                                                                                                                                    <  7:119614/97‑1 (MQ=255)
                                                                       ttAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGttt                                                                                                                                    >  8:119614/1‑97 (MQ=255)
                                                                        tAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAgc                                                                                          <  8:341313/136‑1 (MQ=255)
                                                                        tAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCAAgc                                                                                          >  7:341313/1‑136 (MQ=255)
                                                                          atcatcAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGcgcg                                                                                       <  7:198563/137‑1 (MQ=255)
                                                                          atcatcAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGcgcg                                                                                       >  8:198563/1‑137 (MQ=255)
                                                                           tcatcaACTCCCAAAGGTTCGGTCAGGGGGAGAAATACAAAAAACGTATCTGAAACGTTTATCTCGAACCCAGGCACACTTAGGCCAGGGTTTTAATAATTTGGGCTTTGCGTTGGCAGCCAGT‑‑TGTTCAAGAGAGGAC‑‑‑TCCCGGGc                                                                         >  4:243055/1‑147 (MQ=255)
                                                                            catcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAg                                                                      >  8:26074/1‑149 (MQ=255)
                                                                             atcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTTAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAg                                                                      <  8:312387/148‑1 (MQ=255)
                                                                             atcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAgg                                                                     <  3:48847/149‑1 (MQ=255)
                                                                             atcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAg                                                                      <  1:248534/148‑1 (MQ=255)
                                                                             atcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAg                                                                      <  8:171740/148‑1 (MQ=255)
                                                                                aaCTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGa                                                                  >  8:82343/1‑149 (MQ=255)
                                                                                  cTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATa                                                                <  2:360256/149‑1 (MQ=255)
                                                                                     ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAggg                                                                          >  3:256776/1‑136 (MQ=255)
                                                                                     ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAggg                                                                          <  4:256776/136‑1 (MQ=255)
                                                                                     ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAAcc                                                             >  2:307570/1‑149 (MQ=255)
                                                                                     ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAAcc                                                             <  2:75900/149‑1 (MQ=255)
                                                                                     ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATAATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAAcc                                                             >  7:297018/1‑149 (MQ=255)
                                                                                               ggTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGCCCAGGCGATAACCGATGGCGcc                                                    >  5:286146/1‑148 (MQ=255)
                                                                                                      gggCGAAATACAAAAAACGTATCTTAAACGTTTCTCTCGATCCCCGGCCCCCTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTCGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TACAGGGCCAGGCGAGAACCGCTGGCGCCGAGACCg                                             >  5:286047/1‑148 (MQ=255)
                                                                                                      gggCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGc                                            <  8:335870/149‑1 (MQ=255)
                                                                                                          gAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCg                                             >  1:67626/1‑144 (MQ=255)
                                                                                                               aCAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAccc                                   >  7:171622/1‑149 (MQ=255)
                                                                                                                cAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                  <  5:142916/149‑1 (MQ=255)
                                                                                                                 aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                  >  3:361570/1‑148 (MQ=255)
                                                                                                                 aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                  >  8:342761/1‑148 (MQ=255)
                                                                                                                 aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                  >  4:61361/1‑148 (MQ=255)
                                                                                                                 aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                  >  4:2106/1‑148 (MQ=255)
                                                                                                                  aaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTa                                >  8:277539/1‑149 (MQ=255)
                                                                                                                  aaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTa                                <  6:296338/149‑1 (MQ=255)
                                                                                                                   aaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTAc                               <  2:67626/149‑1 (MQ=255)
                                                                                                                             tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                     >  1:170399/1‑149 (MQ=255)
                                                                                                                             tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                     >  1:27363/1‑149 (MQ=255)
                                                                                                                             tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                     >  3:259569/1‑149 (MQ=255)
                                                                                                                              gAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                     >  7:119729/1‑148 (MQ=255)
                                                                                                                               aaaCGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCgg                    >  1:299017/1‑148 (MQ=255)
                                                                                                                                 aCGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaaa                 <  3:61361/149‑1 (MQ=255)
                                                                                                                                  cGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACTTCGGaaaa                <  4:361570/149‑1 (MQ=255)
                                                                                                                                  cGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaaaa                <  7:342761/149‑1 (MQ=255)
                                                                                                                                  cGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGaaaa                <  6:278264/149‑1 (MQ=255)
                                                                                                                                   gTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAg               >  4:157889/1‑149 (MQ=255)
                                                                                                                                       ctctcGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAg            <  7:26074/148‑1 (MQ=255)
                                                                                                                                             gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAgtgg      >  7:349296/1‑149 (MQ=255)
                                                                                                                                             gAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAg         <  8:171622/146‑1 (MQ=255)
                                                                                                                                               aCCCCGGCCCACTTATGCCAGGTTTTTTATAATTTTGGCTGTGGGTTGGCAGCCATTCGTTTTCAACCGCGCCCCCTTCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGa          <  1:307570/148‑1 (MQ=255)
                                                                                                                                                 cccGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCGGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGAc  >  4:328719/1‑149 (MQ=255)
                                                                                                                                                 cccGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGGGAc  >  4:268995/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                       
CTTGATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGAC  >  NC_002947/650085‑650377

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