Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,259,548 G→T 6.5% D127Y (GAT→TAT)  hemA → glutamyl tRNA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,259,5480GT 6.5% 1010.5 / 31.4 261D127Y (GAT→TAT) hemAglutamyl tRNA reductase
Reads supporting (aligned to +/- strand):  ref base G (127/117);  new base T (8/9);  total (135/126)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 262 total aligned reads displayed.

CAATCTTAATGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATC  >  CP009273/1259400‑1259697
                                                                                                                                                    |                                                                                                                                                     
caATCTTAATGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGa                                                                                                                                                      >  2:2755227/1‑150 (MQ=255)
 aaTCTTAATGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGAt                                                                                                                                                     >  1:2150629/1‑150 (MQ=255)
    cTTAATGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCg                                                                                                                                                  <  1:2801909/150‑1 (MQ=255)
       aaTGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCaa                                                                                                                                               >  2:4252815/1‑150 (MQ=255)
       aaTGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCaa                                                                                                                                               <  2:3830251/150‑1 (MQ=255)
         tGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTGAAAAAAGCGTTTGCCGATTCGCaaaa                                                                                                                                             >  2:2730356/1‑150 (MQ=255)
           aagaagATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAg                                                                                                                                           <  2:274328/150‑1 (MQ=255)
           aagaagATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAg                                                                                                                                           >  1:1508897/1‑150 (MQ=255)
              aagaTCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTc                                                                                                                                        >  1:2413738/1‑150 (MQ=255)
               agaTCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCa                                                                                                                                       <  2:1481065/150‑1 (MQ=255)
               agaTCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGGTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCa                                                                                                                                       <  2:1023954/150‑1 (MQ=255)
                gaTCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCat                                                                                                                                      >  1:1646149/1‑150 (MQ=255)
                 aTCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCTATTCGCAAAAAGGTCata                                                                                                                                     >  1:2488786/1‑150 (MQ=255)
                 aTCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCata                                                                                                                                     >  2:1262172/1‑150 (MQ=255)
                  tCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCatat                                                                                                                                    <  1:1558150/150‑1 (MQ=255)
                   cTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATg                                                                                                                                   <  2:2985284/150‑1 (MQ=255)
                    tGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGa                                                                                                                                  <  2:1691321/150‑1 (MQ=255)
                       gTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAgg                                                                                                                               >  1:2421257/1‑150 (MQ=255)
                         aaaaGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGcc                                                                                                                             >  1:2786589/1‑150 (MQ=255)
                          aaaGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCa                                                                                                                            >  1:4299278/1‑150 (MQ=255)
                            aGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTACTGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGc                                                                                                                          <  1:3069270/150‑1 (MQ=255)
                               ctctACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGaa                                                                                                                       >  1:1871659/1‑150 (MQ=255)
                                tctACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAAc                                                                                                                      <  1:3252388/150‑1 (MQ=255)
                                  tACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTg                                                                                                                    >  2:1400085/1‑150 (MQ=255)
                                    cTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGa                                                                                                                  <  1:898910/150‑1 (MQ=255)
                                    cTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGa                                                                                                                  >  2:3602860/1‑150 (MQ=255)
                                        cATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGc                                                                                                              >  2:1952972/1‑150 (MQ=255)
                                        cATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGc                                                                                                              >  1:1364113/1‑150 (MQ=255)
                                         aTCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCa                                                                                                             >  2:817003/1‑150 (MQ=255)
                                           cAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATg                                                                                                           <  1:3085466/150‑1 (MQ=255)
                                            aGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGt                                                                                                          <  2:2308994/150‑1 (MQ=255)
                                            aGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGt                                                                                                          <  1:3911800/150‑1 (MQ=255)
                                            aGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGt                                                                                                          <  1:2558307/150‑1 (MQ=255)
                                             ggATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGtt                                                                                                         <  2:1646149/150‑1 (MQ=255)
                                                    gacgCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGaaa                                                                                                  >  2:1837528/1‑150 (MQ=255)
                                                          gTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTtc                                                                                            >  2:541387/1‑150 (MQ=255)
                                                            tAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTtctc                                                                                          >  2:983861/1‑150 (MQ=255)
                                                             aGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTctct                                                                                         <  2:2017436/150‑1 (MQ=255)
                                                              gCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTg                                                                                        >  2:3132344/1‑150 (MQ=255)
                                                                cATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCTATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTc                                                                                      <  1:3864903/150‑1 (MQ=255)
                                                                 aTTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTcg                                                                                     <  1:802399/150‑1 (MQ=255)
                                                                  tttAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTcgc                                                                                    >  2:3451272/1‑150 (MQ=255)
                                                                     aaTGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGaa                                                                                 <  2:4266539/150‑1 (MQ=255)
                                                                          gtgtTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCTATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAAcgcg                                                                            <  2:1301472/150‑1 (MQ=255)
                                                                          gtgtTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAAcgcg                                                                            <  2:4308194/150‑1 (MQ=255)
                                                                                ccAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCa                                                                      >  2:742367/1‑150 (MQ=255)
                                                                                 cAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCAc                                                                     >  1:3334743/1‑150 (MQ=255)
                                                                                 cAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCAc                                                                     <  2:2820658/150‑1 (MQ=255)
                                                                                  aGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACt                                                                    <  2:2368059/150‑1 (MQ=255)
                                                                                     ggCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGaa                                                                 >  1:2568357/1‑150 (MQ=255)
                                                                                        cTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACa                                                              >  1:2257070/1‑150 (MQ=255)
                                                                                        cTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACa                                                              >  1:2257487/1‑150 (MQ=255)
                                                                                         tGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAg                                                             >  2:2724332/1‑150 (MQ=255)
                                                                                         tGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTGTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAg                                                             <  2:2923216/150‑1 (MQ=255)
                                                                                           gATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGat                                                           >  1:2668272/1‑150 (MQ=255)
                                                                                           gATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGat                                                           >  2:1698782/1‑150 (MQ=255)
                                                                                           gATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGat                                                           >  2:1698685/1‑150 (MQ=255)
                                                                                            aTTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCTATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGata                                                          <  1:4278650/150‑1 (MQ=255)
                                                                                              tCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATc                                                        <  1:702008/150‑1 (MQ=255)
                                                                                                 cTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGt                                                     <  1:3883734/150‑1 (MQ=255)
                                                                                                  tGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTg                                                    >  1:1180831/1‑150 (MQ=255)
                                                                                                  tGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTg                                                    >  2:2871444/1‑150 (MQ=255)
                                                                                                  tGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTg                                                    >  1:1002097/1‑150 (MQ=255)
                                                                                                  tGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTg                                                    >  1:747006/1‑150 (MQ=255)
                                                                                                       cTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAgc                                               <  2:1990517/150‑1 (MQ=255)
                                                                                                          ggggAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCt                                            >  1:172155/1‑150 (MQ=255)
                                                                                                            ggAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCtgt                                          >  1:677194/1‑150 (MQ=255)
                                                                                                               gCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtc                                       <  2:4164450/150‑1 (MQ=255)
                                                                                                                ccGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtct                                      <  2:3077678/150‑1 (MQ=255)
                                                                                                                ccGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtct                                      <  1:209689/150‑1 (MQ=255)
                                                                                                                ccGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtct                                      <  1:836907/150‑1 (MQ=255)
                                                                                                                 cGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtctg                                     >  2:147520/1‑150 (MQ=255)
                                                                                                                 cGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtctg                                     <  2:297455/150‑1 (MQ=255)
                                                                                                                 cGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGtgtctg                                     >  2:1409261/1‑150 (MQ=255)
                                                                                                                    aGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCg                                  >  1:3667456/1‑150 (MQ=255)
                                                                                                                    aGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCg                                  >  1:2046538/1‑150 (MQ=255)
                                                                                                                          tCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTg                            >  1:1230351/1‑150 (MQ=255)
                                                                                                                            ggTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCg                          <  1:2188801/150‑1 (MQ=255)
                                                                                                                            ggTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCg                          >  1:4292074/1‑150 (MQ=255)
                                                                                                                            ggTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCg                          <  1:4868/150‑1 (MQ=255)
                                                                                                                            ggTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTGTTGCg                          >  2:4094150/1‑150 (MQ=255)
                                                                                                                             gTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCgg                         >  1:663533/1‑150 (MQ=255)
                                                                                                                              tCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGc                        >  1:3465733/1‑150 (MQ=255)
                                                                                                                                aGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCtt                      <  1:1131866/150‑1 (MQ=255)
                                                                                                                                aGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCtt                      <  1:935763/150‑1 (MQ=255)
                                                                                                                                aGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCtt                      <  1:519358/150‑1 (MQ=255)
                                                                                                                                aGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCtt                      >  1:565058/1‑150 (MQ=255)
                                                                                                                                  gTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGt                    >  1:1040418/1‑150 (MQ=255)
                                                                                                                                   ttAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTa                   <  2:2421247/150‑1 (MQ=255)
                                                                                                                                   ttAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTa                   <  1:382862/150‑1 (MQ=255)
                                                                                                                                   ttAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTa                   <  2:263255/150‑1 (MQ=255)
                                                                                                                                     aaaaaaGCGTTTGCCTATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACg                 >  2:433406/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACg                 >  1:3379383/1‑150 (MQ=255)
                                                                                                                                       aaaaGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTATATCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGAGCCAGCGCTGTGTCTGTCGCTTTTTCGGCTTGTACGCt               >  2:4032538/1‑150 (MQ=255)
                                                                                                                                         aaGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTgg             >  1:1749971/1‑150 (MQ=255)
                                                                                                                                           gCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGgcg           <  1:3552150/150‑1 (MQ=255)
                                                                                                                                               ttGCCTATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGc       >  1:2189535/1‑150 (MQ=255)
                                                                                                                                                 gCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAg     <  1:1377575/150‑1 (MQ=255)
                                                                                                                                                   cGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGAt   >  2:3612415/1‑150 (MQ=255)
                                                                                                                                                   cGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGAt   >  2:1791875/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAATCTTAATGAAGAAGATCTGCGTAAAAGCCTCTACTGGCATCAGGATAACGACGCGGTTAGCCATTTAATGCGTGTTGCCAGCGGCCTGGATTCACTGGTTCTGGGGGAGCCGCAGATCCTCGGTCAGGTTAAAAAAGCGTTTGCCGATTCGCAAAAAGGTCATATGAAGGCCAGCGAACTGGAACGCATGTTCCAGAAATCTTTCTCTGTCGCGAAACGCGTTCGCACTGAAACAGATATCGGTGCCAGCGCTGTGTCTGTCGCTTTTGCGGCTTGTACGCTGGCGCGGCAGATC  >  CP009273/1259400‑1259697

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: