Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,983,298 G→T 29.9% R277S (CGC→AGC)  hemC ← hydroxymethylbilane synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,983,2980GT29.9% 557.9 / 209.6 287R277S (CGC→AGC) hemChydroxymethylbilane synthase
Reads supporting (aligned to +/- strand):  ref base G (101/100);  new base T (49/37);  total (150/137)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 301 total aligned reads displayed.

TCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGC  >  CP009273/3983151‑3983447
                                                                                                                                                   |                                                                                                                                                     
tctcCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCgctgc                                                                                                                                                     <  1:4744150/150‑1 (MQ=255)
    cAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTc                                                                                                                                                 <  2:795648/150‑1 (MQ=255)
     aGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCa                                                                                                                                                <  2:2120775/150‑1 (MQ=255)
     aGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCa                                                                                                                                                >  1:1623952/1‑150 (MQ=255)
      gCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCAc                                                                                                                                               <  2:4709002/150‑1 (MQ=255)
        gggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcg                                                                                                                                             >  2:1051908/1‑150 (MQ=255)
        gggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcg                                                                                                                                             >  2:56831/1‑150 (MQ=255)
        gggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcg                                                                                                                                             >  1:4530704/1‑150 (MQ=255)
         ggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcgc                                                                                                                                            <  1:2300449/150‑1 (MQ=255)
         ggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcgc                                                                                                                                            >  1:4309777/1‑150 (MQ=255)
         ggAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACcgc                                                                                                                                            >  1:3655269/1‑150 (MQ=255)
          gagaCGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACcgcg                                                                                                                                           >  2:3674679/1‑150 (MQ=255)
             aCGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGaat                                                                                                                                        <  1:304191/150‑1 (MQ=255)
              cgggcgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGaata                                                                                                                                       <  2:1837000/150‑1 (MQ=255)
                ggcgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGaataat                                                                                                                                     <  1:1234964/150‑1 (MQ=255)
                  cgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCt                                                                                                                                   <  2:2556977/150‑1 (MQ=255)
                  cgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGaataatat                                                                                                                                   >  1:3240433/1‑148 (MQ=255)
                  cgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCt                                                                                                                                   <  2:647967/150‑1 (MQ=255)
                  cgggTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCt                                                                                                                                   <  1:4274462/150‑1 (MQ=255)
                        aCAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAAc                                                                                                                             <  1:4763425/150‑1 (MQ=255)
                           aGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGt                                                                                                                          >  2:416892/1‑150 (MQ=255)
                           aGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGt                                                                                                                          <  1:2533684/150‑1 (MQ=255)
                            ggATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTc                                                                                                                         >  1:1494608/1‑150 (MQ=255)
                             gATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTcc                                                                                                                        <  1:4863506/150‑1 (MQ=255)
                              aTACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCg                                                                                                                       >  1:3812746/1‑150 (MQ=255)
                              aTACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCg                                                                                                                       <  2:1794367/150‑1 (MQ=255)
                              aTACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCg                                                                                                                       <  1:1996042/150‑1 (MQ=255)
                                       ccGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGa                                                                                                              >  1:4316674/1‑150 (MQ=255)
                                         ggggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGAcc                                                                                                            <  2:1813209/150‑1 (MQ=255)
                                         ggggCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAAACCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGccc                                                                                                            >  2:2310766/1‑150 (MQ=255)
                                           ggCGTCACCGTTATAGACTTCATCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATATGTTCGGCAACTTTCGGCGCACCTCGGCGTTCACCGCGAATAATCTGCGAACCGTCCCGCTCGCCGACCAg                                                                                                          >  1:4976288/1‑150 (MQ=255)
                                            gCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgc                                                                                                         <  2:2367454/150‑1 (MQ=255)
                                             cGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcg                                                                                                        <  1:4998520/150‑1 (MQ=255)
                                              gTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc                                                                                                       >  1:2670270/1‑150 (MQ=255)
                                              gTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAgcgc                                                                                                       <  1:916878/150‑1 (MQ=255)
                                                ctcCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                                                     <  2:489452/150‑1 (MQ=255)
                                                ctcCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                                                     <  1:1276636/150‑1 (MQ=255)
                                                ctcCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                                                     >  2:3466370/1‑150 (MQ=255)
                                                ctcCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                                                     >  2:3912306/1‑150 (MQ=255)
                                                ctcCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGcgcac                                                                                                     >  2:3942396/1‑150 (MQ=255)
                                                    gTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAg                                                                                                 >  2:3151124/1‑150 (MQ=255)
                                                     ttATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                                                >  2:2051610/1‑150 (MQ=255)
                                                     ttATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCAACGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGc                                                                                                >  1:4586045/1‑150 (MQ=255)
                                                         aGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGa                                                                                            >  1:3135118/1‑150 (MQ=255)
                                                           aCTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAtt                                                                                          >  2:3859287/1‑150 (MQ=255)
                                                            cTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGAttt                                                                                         >  2:5127559/1‑150 (MQ=255)
                                                                aGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGcc                                                                                     >  1:1818198/1‑150 (MQ=255)
                                                                     ggATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCaa                                                                                >  2:3448557/1‑150 (MQ=255)
                                                                      gATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAAt                                                                               <  2:3784328/150‑1 (MQ=255)
                                                                       aTCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATa                                                                              >  2:2250480/1‑150 (MQ=255)
                                                                       aTCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCAATTTCATCGTCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCAGCGAACCGTCCTGCGCGCCGACCAGCGCACGCAGACAGATTTCGCCATCAATa                                                                              >  1:4252498/1‑150 (MQ=255)
                                                                         ctcGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAg                                                                            <  1:7580/150‑1 (MQ=255)
                                                                           cgcgcgcgCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCt                                                                          <  1:54612/150‑1 (MQ=255)
                                                                             cgcgcgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCTTCCGTCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCg                                                                        <  2:1352174/150‑1 (MQ=255)
                                                                             cgcgcgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCg                                                                        <  2:1215431/150‑1 (MQ=255)
                                                                               cgcgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGc                                                                      >  2:2449308/1‑150 (MQ=255)
                                                                               cgcgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGc                                                                      >  1:5314304/1‑150 (MQ=255)
                                                                                 cgcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGt                                                                    >  2:2613252/1‑150 (MQ=255)
                                                                                  gcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTa                                                                   <  1:662546/150‑1 (MQ=255)
                                                                                  gcCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTa                                                                   >  1:361270/1‑150 (MQ=255)
                                                                                    cGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGc                                                                 <  2:1557126/150‑1 (MQ=255)
                                                                                    cGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGc                                                                 >  1:515867/1‑150 (MQ=255)
                                                                                            aGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCCGATTTCGCCATCACTAAGCTCGGCGTAGCTACCAAta                                                         >  2:2524554/1‑149 (MQ=255)
                                                                                                gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCa                                                     >  1:5221900/1‑150 (MQ=255)
                                                                                                gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCa                                                     <  1:2848701/150‑1 (MQ=255)
                                                                                                gCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCa                                                     >  1:4372367/1‑150 (MQ=255)
                                                                                                 ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGATCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCTATTGGCAc                                                    >  2:5045110/1‑150 (MQ=255)
                                                                                                 ctctTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCAc                                                    >  1:3513177/1‑150 (MQ=255)
                                                                                                   cttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCt                                                  >  1:2148373/1‑150 (MQ=255)
                                                                                                    ttctGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTg                                                 >  2:4663636/1‑150 (MQ=255)
                                                                                                              gAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCt                                       <  2:5113890/150‑1 (MQ=255)
                                                                                                               aaaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCtt                                      >  2:692152/1‑150 (MQ=255)
                                                                                                                aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCTCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTCACAATTGTCACCTGACATCCGCCTTc                                     >  2:292992/1‑150 (MQ=255)
                                                                                                                aaTCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTc                                     <  2:5226715/150‑1 (MQ=255)
                                                                                                                  tCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCga                                   <  1:1103960/150‑1 (MQ=255)
                                                                                                                   ccccATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgag                                  >  2:782496/1‑150 (MQ=255)
                                                                                                                    cccATTTGTTCGGCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCgaga                                 >  2:3085068/1‑150 (MQ=255)
                                                                                                                      cATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACg                               >  2:270817/1‑150 (MQ=255)
                                                                                                                         ttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGt                            <  2:2330760/150‑1 (MQ=255)
                                                                                                                         ttGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGt                            <  1:2613252/150‑1 (MQ=255)
                                                                                                                           gTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTAt                          <  1:416892/150‑1 (MQ=255)
                                                                                                                            ttCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTAtt                         <  2:4291884/150‑1 (MQ=255)
                                                                                                                              cGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCa                       >  1:4684397/1‑150 (MQ=255)
                                                                                                                               ggCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCAt                      <  2:4957293/150‑1 (MQ=255)
                                                                                                                               ggCATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCAt                      >  2:1474277/1‑150 (MQ=255)
                                                                                                                                 cATCTTGCGGCGCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATgg                    >  2:1947700/1‑150 (MQ=255)
                                                                                                                                 cATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATgg                    >  1:3369734/1‑150 (MQ=255)
                                                                                                                                 cATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATgg                    <  2:3238371/150‑1 (MQ=255)
                                                                                                                                  aTCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGgc                   >  1:370289/1‑150 (MQ=255)
                                                                                                                                       gcggcgCACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCGGGCGCGGCGAGCCGCGCACGCAGCCAGCTTTCGCCATCCATAAGGTCGGCGTAGCTACCAATTGGCACCTGACTTCCGCCTTCGAGACGGGTATTCATGGCGCGtt              <  2:1519181/150‑1 (MQ=255)
                                                                                                                                         ggcgCACCGCGGCGGTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCGACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCt            <  1:2926690/150‑1 (MQ=255)
                                                                                                                                          gcgcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTg           <  1:768737/150‑1 (MQ=255)
                                                                                                                                           cgcACCGCTGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGc          >  1:3288316/1‑150 (MQ=255)
                                                                                                                                            gcACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCg         <  1:3942396/150‑1 (MQ=255)
                                                                                                                                              aCCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGt       <  1:3327602/150‑1 (MQ=255)
                                                                                                                                                cgctGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTaa     >  1:1433122/1‑150 (MQ=255)
                                                                                                                                                cgcgGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTaa     >  2:5157127/1‑150 (MQ=255)
                                                                                                                                                 gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAAc    <  1:56831/150‑1 (MQ=255)
                                                                                                                                                 gcggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAAc    <  2:1623952/150‑1 (MQ=255)
                                                                                                                                                   ggcgTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAAcgc  >  2:2874939/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TCTCCAGCGGGAGACGGGCGGGTGACAAGGATACTCATGCCGGGGCGTCTCCGTTATAGACTTCAGCGAGGATCTCGCGCGCGCCGTTATTCAGTAGCTCTTCTGCCAGCGAAATCCCCATTTGTTCGGCATCTTGCGGCGCACCGCGGCGTTCACCGCGAATAATCTGCGAACCGTCCGGCGCGCCGACCAGCGCACGCAGCCAGATTTCGCCATCAATAAGCTCGGCGTAGCTACCAATTGGCACCTGACATCCGCCTTCGAGACGGGTATTCATGGCGCGTTCTGCGGTAACGC  >  CP009273/3983151‑3983447

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: ‑