Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,239 T→G 8.5% intergenic (‑83/+17) 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,2390TG 8.5% 242.5 / 13.2 83intergenic (‑83/+17)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (39/36);  new base G (2/5);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

TCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGAC  >  NC_002947/650092‑650377
                                                                                                                                                    |                                                                                                                                               
tcCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACt                                                                                                                                                 >  3:190800/1‑148 (MQ=255)
  cAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACACCCAAAGGTTCGGTCAGGGGGCTAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGCAACACTTAg                                                                                                                                              >  2:331587/1‑149 (MQ=255)
        gtcgATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAgg                                                                                                                                        <  5:229626/149‑1 (MQ=255)
        gtcgATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAgg                                                                                                                                        <  8:129230/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                                                                                                                    <  5:83857/149‑1 (MQ=255)
                            gcagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtgt                                                                                                                    <  1:149075/149‑1 (MQ=255)
                            gcagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtgt                                                                                                                    >  2:248534/1‑149 (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                                                                                                                                    >  8:119614/1‑97 (MQ=255)
                                                                ttAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGttt                                                                                                                                    <  7:119614/97‑1 (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                                                                      <  8:171740/148‑1 (MQ=255)
                                                                      atcaACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAg                                                                      <  1:248534/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                                                                          <  4:256776/136‑1 (MQ=255)
                                                                              ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAggg                                                                          >  3:256776/1‑136 (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                                  >  4:61361/1‑148 (MQ=255)
                                                                                                          aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                  >  8:342761/1‑148 (MQ=255)
                                                                                                          aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                  >  3:361570/1‑148 (MQ=255)
                                                                                                          aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                  >  4:2106/1‑148 (MQ=255)
                                                                                                           aaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTa                                <  6:296338/149‑1 (MQ=255)
                                                                                                           aaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTa                                >  8:277539/1‑149 (MQ=255)
                                                                                                            aaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTAc                               <  2:67626/149‑1 (MQ=255)
                                                                                                                      tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                     >  3:259569/1‑149 (MQ=255)
                                                                                                                      tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                     >  1:170399/1‑149 (MQ=255)
                                                                                                                      tGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCg                     >  1:27363/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)
                                                                                                                                                    |                                                                                                                                               
TCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGAC  >  NC_002947/650092‑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: