Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,783,091 G→T 7.3% E134* (GAG→TAG)  aroH → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,783,0910GT 7.3% 1169.4 / 39.7 309E134* (GAG→TAG) aroH3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base G (174/112);  new base T (5/18);  total (179/130)
Fisher's exact test for biased strand distribution p-value = 3.22e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

CACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGT  >  CP009273/1782942‑1783240
                                                                                                                                                     |                                                                                                                                                     
cACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATg                                                                                                                                                       >  2:1472181/1‑150 (MQ=255)
 aCGGCTGGAAATTGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGa                                                                                                                                                      >  1:3460994/1‑150 (MQ=255)
    gCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCt                                                                                                                                                   >  1:1267295/1‑150 (MQ=255)
     cTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTg                                                                                                                                                  >  1:3539661/1‑150 (MQ=255)
       ggAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTggg                                                                                                                                                >  1:2903780/1‑150 (MQ=255)
         aaaTCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCg                                                                                                                                              <  1:3028298/150‑1 (MQ=255)
           aTCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTc                                                                                                                                            >  2:120327/1‑150 (MQ=255)
                 aTGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACTGCAGCTATCTGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGATGGGCGTCCCaaac                                                                                                                                      >  2:417013/1‑148 (MQ=255)
                 aTGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAAcc                                                                                                                                      >  2:1298607/1‑150 (MQ=255)
                    cgcACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACcgcg                                                                                                                                   >  2:3565132/1‑150 (MQ=255)
                    cgcACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACcgcg                                                                                                                                   <  1:2425407/150‑1 (MQ=255)
                    cgcACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACcgcg                                                                                                                                   >  1:2732373/1‑150 (MQ=255)
                       aCCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGAcc                                                                                                                                >  1:875836/1‑150 (MQ=255)
                       aCCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGAcc                                                                                                                                <  2:1813896/150‑1 (MQ=255)
                       aCCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGAcc                                                                                                                                >  1:1235730/1‑150 (MQ=255)
                         cTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGa                                                                                                                              <  2:2371211/150‑1 (MQ=255)
                                gAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTc                                                                                                                       >  2:965146/1‑150 (MQ=255)
                                 aaaaaCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCg                                                                                                                      <  1:3701172/150‑1 (MQ=255)
                                     aCCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGatat                                                                                                                  >  2:2512860/1‑150 (MQ=255)
                                       cACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATgg                                                                                                                >  2:583082/1‑150 (MQ=255)
                                            aCTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGAcc                                                                                                           <  1:509530/150‑1 (MQ=255)
                                             cTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCg                                                                                                          <  1:2467949/150‑1 (MQ=255)
                                             cTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCg                                                                                                          >  1:4201498/1‑150 (MQ=255)
                                             cTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCg                                                                                                          >  1:3400922/1‑150 (MQ=255)
                                             cTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCg                                                                                                          >  2:4480970/1‑150 (MQ=255)
                                                ttgtCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTc                                                                                                       <  2:1079062/150‑1 (MQ=255)
                                                gtgtCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTc                                                                                                       <  2:2751871/149‑1 (MQ=255)
                                                      gCTGGAAAGGACTAATCTCCGATCCAGATTTAACCGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTa                                                                                                 <  1:940125/150‑1 (MQ=255)
                                                         ggAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTg                                                                                              <  2:4143266/150‑1 (MQ=255)
                                                           aaaGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCt                                                                                            >  1:263908/1‑150 (MQ=255)
                                                            aaGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTg                                                                                           >  1:2580488/1‑150 (MQ=255)
                                                                aCTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGAttt                                                                                       <  1:3111557/150‑1 (MQ=255)
                                                                   aaTCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAAt                                                                                    >  1:2064460/1‑150 (MQ=255)
                                                                   aaTCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAAt                                                                                    <  2:3308751/150‑1 (MQ=255)
                                                                    aTCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATc                                                                                   >  1:4558657/1‑150 (MQ=255)
                                                                          gATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTgg                                                                             >  1:4077907/1‑150 (MQ=255)
                                                                           aTCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTggg                                                                            >  2:1610974/1‑150 (MQ=255)
                                                                           aTCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTggg                                                                            <  2:3525373/150‑1 (MQ=255)
                                                                           aTCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTggg                                                                            <  2:3452460/150‑1 (MQ=255)
                                                                            tCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTgggg                                                                           >  1:1129870/1‑150 (MQ=255)
                                                                             ccAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGgc                                                                          <  2:4323454/150‑1 (MQ=255)
                                                                                gATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGcgcg                                                                       <  2:2761381/150‑1 (MQ=255)
                                                                                gATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGcgcg                                                                       <  2:3204491/150‑1 (MQ=255)
                                                                                 aTTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGa                                                                      <  2:2690490/150‑1 (MQ=255)
                                                                                     aaaCGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGGTTATTGCTGATTTAATCAGTTGGGGCGCGATTgg                                                                  >  2:855173/1‑150 (MQ=255)
                                                                                     aaaCGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGAGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTgg                                                                  >  2:3078446/1‑150 (MQ=255)
                                                                                        cGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGgcgc                                                               >  1:2888700/1‑150 (MQ=255)
                                                                                         ggCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGTATTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCa                                                              >  2:2232616/1‑150 (MQ=255)
                                                                                               tATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtact                                                        <  1:4307363/150‑1 (MQ=255)
                                                                                               tATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtact                                                        >  2:3108315/1‑150 (MQ=255)
                                                                                                aTCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtacta                                                       <  2:3283470/150‑1 (MQ=255)
                                                                                                aTCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGtacta                                                       <  1:3697272/150‑1 (MQ=255)
                                                                                                  cGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTAcc                                                     >  2:3426852/1‑150 (MQ=255)
                                                                                                  cGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTAcc                                                     <  2:4385904/150‑1 (MQ=255)
                                                                                                  cGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTAcc                                                     >  2:313997/1‑150 (MQ=255)
                                                                                                   gggTAAATCACGGTCTTGAGCTGGCGCGAAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCg                                                    >  2:1235260/1‑150 (MQ=255)
                                                                                                   gggTAAATCACGGTCTGGAGCTGGTGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCg                                                    >  2:716898/1‑150 (MQ=255)
                                                                                                   gggTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTTATTTAATCAGTAGGTGCGCGATTGGCGCAAGTACTACCg                                                    >  2:4188595/1‑150 (MQ=255)
                                                                                                   gggTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCg                                                    >  2:4092708/1‑150 (MQ=255)
                                                                                                       aaaTCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAg                                                <  2:1680178/150‑1 (MQ=255)
                                                                                                          tCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCa                                             >  1:3423049/1‑150 (MQ=255)
                                                                                                          tCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCa                                             >  1:3329410/1‑150 (MQ=255)
                                                                                                           cACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGGTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAg                                            >  2:3162300/1‑150 (MQ=255)
                                                                                                             cGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGAt                                          >  2:299316/1‑150 (MQ=255)
                                                                                                             cGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGAt                                          >  2:1217911/1‑150 (MQ=255)
                                                                                                              ggTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATc                                         <  2:4036146/150‑1 (MQ=255)
                                                                                                              ggTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATc                                         <  2:968984/150‑1 (MQ=255)
                                                                                                                 cTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGTGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCAc                                      >  1:1611312/1‑150 (MQ=255)
                                                                                                                 cTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCAc                                      >  1:2297414/1‑150 (MQ=255)
                                                                                                                  tGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCAcc                                     <  2:2077616/150‑1 (MQ=255)
                                                                                                                  tGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCAcc                                     >  2:3834695/1‑150 (MQ=255)
                                                                                                                  tGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCAcc                                     >  1:784537/1‑150 (MQ=255)
                                                                                                                     aGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACcgcg                                  <  2:3705591/150‑1 (MQ=255)
                                                                                                                     aGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACcgcg                                  <  2:2867378/150‑1 (MQ=255)
                                                                                                                       cTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGGTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGaa                                >  2:1110429/1‑150 (MQ=255)
                                                                                                                        tGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGaaa                               >  2:2370551/1‑150 (MQ=255)
                                                                                                                        tGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGaaa                               >  2:1983958/1‑150 (MQ=255)
                                                                                                                             cgcAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCt                          <  1:165228/150‑1 (MQ=255)
                                                                                                                              gcAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCtt                         <  2:1023206/150‑1 (MQ=255)
                                                                                                                              gcAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCtt                         <  1:3929332/150‑1 (MQ=255)
                                                                                                                              gcAAATTACTTTTACAGGTAAATGAACTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCtt                         <  1:3780848/150‑1 (MQ=255)
                                                                                                                                aaaTTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCg                       <  1:1492678/150‑1 (MQ=255)
                                                                                                                                   ttACTTTTACAGGTAAATTAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCa                    <  2:3028017/150‑1 (MQ=255)
                                                                                                                                    tACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCAc                   >  2:3675680/1‑150 (MQ=255)
                                                                                                                                     aCTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCAct                  <  2:1543339/150‑1 (MQ=255)
                                                                                                                                     aCTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCAct                  >  1:1675377/1‑150 (MQ=255)
                                                                                                                                         ttACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCTGCACTCTcc              <  2:2286935/150‑1 (MQ=255)
                                                                                                                                          tACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCt             <  1:3008798/150‑1 (MQ=255)
                                                                                                                                          tACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCt             <  1:358183/150‑1 (MQ=255)
                                                                                                                                            cAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGt           >  1:484749/1‑150 (MQ=255)
                                                                                                                                             aGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTc          >  2:926165/1‑150 (MQ=255)
                                                                                                                                             aGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTc          >  1:1023611/1‑150 (MQ=255)
                                                                                                                                              ggTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTcc         >  2:424795/1‑150 (MQ=255)
                                                                                                                                              ggTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTcc         <  2:4536016/150‑1 (MQ=255)
                                                                                                                                               gTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCg        >  1:2111046/1‑150 (MQ=255)
                                                                                                                                               gTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCg        <  2:2577619/150‑1 (MQ=255)
                                                                                                                                                tAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCgg       >  2:134677/1‑150 (MQ=255)
                                                                                                                                                 aaaTGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGt      <  1:1670472/150‑1 (MQ=255)
                                                                                                                                                 aaaTGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGt      <  2:1919334/150‑1 (MQ=255)
                                                                                                                                                   aTGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAg    <  1:555877/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACGGCTGGAAATCGTAATGCGCACCTATTTTGAAAAACCACGAACTGTTGTCGGCTGGAAAGGACTAATCTCCGATCCAGATTTAAACGGCAGCTATCGGGTAAATCACGGTCTGGAGCTGGCGCGCAAATTACTTTTACAGGTAAATGAGCTGGGCGTCCCAACCGCGACCGAGTTCCTCGATATGGTGACCGGTCAGTTTATTGCTGATTTAATCAGTTGGGGCGCGATTGGCGCACGTACTACCGAAAGTCAGATCCACCGCGAAATGGCTTCGGCACTCTCCTGTCCGGTAGGT  >  CP009273/1782942‑1783240

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