Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 650,236 C→T 8.2% 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 8.2% 242.3 / 8.8 86intergenic (‑80/+20)accB/aroQ‑Iacetyl‑CoA carboxylase biotin carboxyl carrier subunit/type II 3‑dehydroquinate dehydratase
Reads supporting (aligned to +/- strand):  ref base C (43/35);  new base T (2/5);  total (45/41)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.43e-01

ATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGAC  >  NC_002947/650089‑650377
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                                        <  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                                                                                               >  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                                                                                       <  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                                                                          >  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)
                                                                                                                                                    |                                                                                                                                                  
ATCTCCAGTTCGTCGATGCCAGACTCTTCCAGCAGCTCGATCAGTTTCTTGACTTTACGGATATCCATTAATCATCAACTCCCAAAGGTTCGGTCAGGGGGCGAAATACAAAAAACGTATCTGAAACGTTTCTCTCGAACCCCGGCCCACTTAGGCCAGGGTTTTCATCATTTGGGCTGTGCGTTGGCAGCCAGT‑‑TGTTCCAGCGCGGAC‑‑‑TCCAGGGCCAGGCGATAACCGCTGGCGCCCAGACCGCAGATCACCCCTACGGCAACATCGGAAAAGTAGGAGTGGTGAC  >  NC_002947/650089‑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: