Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,246 A→G 8.1% intergenic (‑90/+10) accB ← / ← aroQ‑I acetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947650,2460AG 8.1% 243.0 / 9.3 86intergenic (‑90/+10)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (40/39);  new base G (2/5);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGTTTG  >  NC_002947/650100‑650394
                                                                                                                                                   |                                                                                                                                                         
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                                                                                                                                     <  5:83857/149‑1 (MQ=255)
                    gcagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtgt                                                                                                                                     >  2:248534/1‑149 (MQ=255)
                    gcagcTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCtgt                                                                                                                                     <  1:149075/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                                                                                                                                                  >  3:315216/1‑124 (MQ=255)
                            aTCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCTGGCCTAAGTGGGCCGGGGttcgag                                                                                                                                                  <  4:315216/128‑5 (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                                                                                                                >  2:73355/1‑134 (MQ=255)
                                                      cATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAgt‑‑tgtt                                                                                                                <  1:73355/134‑1 (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                                                                                                        >  8:198563/1‑137 (MQ=255)
                                                           atcatcAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGcgcg                                                                                                        <  7:198563/137‑1 (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:75900/149‑1 (MQ=255)
                                                                      ccAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAAcc                                                                              >  2:307570/1‑149 (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                                                   >  4:2106/1‑148 (MQ=255)
                                                                                                  aaaaaaCGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCAcccc                                                   >  3:361570/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                                      >  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)
                                                                                                                                                gCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCATACGGCAACATCGGAAAAGTAGGAGTGGTGACGGACAGGTTCGCGt     >  7:63122/1‑149 (MQ=255)
                                                                                                                                                  gaGGGTTTTCTTCATTTGGGCTGTGCGTTGGCAGGCTGT‑‑TGTTCCAGCGCGGAC‑‑‑TGCAGGGCCAGGCGATAACCGCTGGCGCCCAGACGGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGttt   <  6:180206/148‑1 (MQ=255)
                                                                                                                                                  cAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGttt   <  3:2106/149‑1 (MQ=255)
                                                                                                                                                  cAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGttt   <  2:27363/149‑1 (MQ=255)
                                                                                                                                                  cAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGttt   <  2:170399/149‑1 (MQ=255)
                                                                                                                                                   aGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGTTtg  <  7:82343/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                         
GTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGACGGAACGGTTCGCGTTTG  >  NC_002947/650100‑650394

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: