Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,883,189 C→A 5.4% V90V (GTG→GTT cysI ← sulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,883,1890CA 5.4% 1412.1 / 34.8 353V90V (GTG→GTTcysIsulfite reductase, beta subunit, NAD(P)‑binding, heme‑binding
Reads supporting (aligned to +/- strand):  ref base C (169/165);  new base A (12/7);  total (181/172)
Fisher's exact test for biased strand distribution p-value = 3.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.47e-01

Only 100 of 370 total aligned reads displayed.

ACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCTTAAAGCCGCCGGT  >  CP009273/2883040‑2883338
                                                                                                                                                     |                                                                                                                                                     
accgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATc                                                                                                                                                       <  2:2621201/150‑1 (MQ=255)
accgaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATc                                                                                                                                                       <  2:4164973/150‑1 (MQ=255)
   gaGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcacc                                                                                                                                                    <  1:1459574/150‑1 (MQ=255)
     gTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcacccc                                                                                                                                                  <  1:2442405/150‑1 (MQ=255)
       gcagcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcaccccac                                                                                                                                                >  1:2248060/1‑150 (MQ=255)
        cagcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcaccccacc                                                                                                                                               <  2:1821954/150‑1 (MQ=255)
        cagcaTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATcaccccacc                                                                                                                                               >  2:4584934/1‑150 (MQ=255)
           caTCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCgg                                                                                                                                            >  1:2203238/1‑150 (MQ=255)
             tCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCa                                                                                                                                          <  2:821244/150‑1 (MQ=255)
                   gCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAgacgac                                                                                                                                    <  2:2496662/150‑1 (MQ=255)
                    cACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACa                                                                                                                                   <  2:595282/150‑1 (MQ=255)
                    cACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACa                                                                                                                                   <  2:4888222/150‑1 (MQ=255)
                       cGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCg                                                                                                                                >  1:1363937/1‑150 (MQ=255)
                        ggTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGa                                                                                                                               <  2:1555634/150‑1 (MQ=255)
                         gTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGaa                                                                                                                              >  2:1495713/1‑150 (MQ=255)
                           ttGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAagc                                                                                                                            >  2:1066060/1‑150 (MQ=255)
                            tGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAagca                                                                                                                           >  1:43412/1‑150 (MQ=255)
                             gACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAagcag                                                                                                                          >  2:4004065/1‑150 (MQ=255)
                               cGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAgcagca                                                                                                                        <  1:4818336/150‑1 (MQ=255)
                               cGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAgcagca                                                                                                                        >  1:3999729/1‑150 (MQ=255)
                                 ttCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATc                                                                                                                      >  2:2822777/1‑150 (MQ=255)
                                      tttCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTg                                                                                                                 >  2:2238205/1‑150 (MQ=255)
                                        tCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCAGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGgc                                                                                                               >  1:522756/1‑150 (MQ=255)
                                                ccGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTcc                                                                                                       <  1:2730327/150‑1 (MQ=255)
                                                 cGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCa                                                                                                      <  1:1989362/150‑1 (MQ=255)
                                                 cGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCa                                                                                                      <  1:1654580/150‑1 (MQ=255)
                                                  gTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAg                                                                                                     <  2:2059542/150‑1 (MQ=255)
                                                     gAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATAACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCtt                                                                                                  >  1:3279957/1‑150 (MQ=255)
                                                      aaCTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATAACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTc                                                                                                 >  1:3259562/1‑150 (MQ=255)
                                                        cTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTg                                                                                               <  2:3720261/150‑1 (MQ=255)
                                                        cTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTg                                                                                               >  1:2057000/1‑150 (MQ=255)
                                                         tGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGt                                                                                              >  1:4099013/1‑150 (MQ=255)
                                                            aaCGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCa                                                                                           >  1:4243919/1‑150 (MQ=255)
                                                             aCGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAg                                                                                          <  2:1789510/150‑1 (MQ=255)
                                                             aCGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAg                                                                                          <  2:4067173/150‑1 (MQ=255)
                                                               gTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCa                                                                                        <  1:1790913/150‑1 (MQ=255)
                                                                tCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCAc                                                                                       <  1:181884/150‑1 (MQ=255)
                                                                  tGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGt                                                                                     >  1:1730089/1‑150 (MQ=255)
                                                                    gCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTc                                                                                   <  1:4294662/150‑1 (MQ=255)
                                                                     cGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCg                                                                                  <  2:2005625/150‑1 (MQ=255)
                                                                      ggTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCgg                                                                                 <  2:2131016/150‑1 (MQ=255)
                                                                       gTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGgc                                                                                >  2:4713306/1‑150 (MQ=255)
                                                                         tGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGgcgc                                                                              >  1:3392121/1‑150 (MQ=255)
                                                                            tCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGa                                                                           <  1:4713398/150‑1 (MQ=255)
                                                                                gCGAATGCTGCCCGAGCTGGTGTTGTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCCCCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTc                                                                       <  2:3688040/150‑1 (MQ=255)
                                                                                gCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTc                                                                       <  2:2285960/150‑1 (MQ=255)
                                                                                gCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTc                                                                       <  1:548324/150‑1 (MQ=255)
                                                                                gCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATAACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTc                                                                       <  1:2594226/150‑1 (MQ=255)
                                                                                 cGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTcg                                                                      <  2:4040027/150‑1 (MQ=255)
                                                                                      gctgcCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTc                                                                 >  1:1783762/1‑150 (MQ=255)
                                                                                      gctgcCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTc                                                                 <  2:3879970/150‑1 (MQ=255)
                                                                                      gctgcCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTc                                                                 >  2:3312852/1‑150 (MQ=255)
                                                                                        tgcCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCAt                                                               <  1:2917636/150‑1 (MQ=255)
                                                                                         gcCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATc                                                              <  2:2373831/150‑1 (MQ=255)
                                                                                          ccATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATcc                                                             >  1:123930/1‑150 (MQ=255)
                                                                                           cATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCct                                                            <  1:1893845/150‑1 (MQ=255)
                                                                                           cATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCct                                                            <  2:1842639/150‑1 (MQ=255)
                                                                                           cATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCct                                                            >  1:2545584/1‑150 (MQ=255)
                                                                                            aTAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctg                                                           >  2:1859204/1‑150 (MQ=255)
                                                                                             tAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgc                                                          <  2:2537555/150‑1 (MQ=255)
                                                                                              aGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct                                                         >  1:1917830/1‑150 (MQ=255)
                                                                                              aGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCctgct                                                         <  2:3175764/150‑1 (MQ=255)
                                                                                                aTGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGa                                                       >  2:3477218/1‑150 (MQ=255)
                                                                                                    tgtTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACa                                                   <  2:4317320/150‑1 (MQ=255)
                                                                                                     gtTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACAt                                                  >  2:202636/1‑150 (MQ=255)
                                                                                                     gtTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACAt                                                  <  1:4490217/150‑1 (MQ=255)
                                                                                                      ttttCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATg                                                 >  1:2617771/1‑150 (MQ=255)
                                                                                                       tttCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGc                                                <  1:1066060/150‑1 (MQ=255)
                                                                                                        ttCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGcc                                               >  2:3599091/1‑150 (MQ=255)
                                                                                                         tCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCg                                              <  1:3709727/150‑1 (MQ=255)
                                                                                                          cACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGt                                             <  2:4405461/150‑1 (MQ=255)
                                                                                                          cACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGt                                             >  1:3650336/1‑150 (MQ=255)
                                                                                                           cccGGCAAATTTGGCGATCGCCTGCCACTGTTTAGTGGTAATAACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTg                                            >  1:4595144/2‑150 (MQ=255)
                                                                                                             cGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGa                                          <  1:381993/150‑1 (MQ=255)
                                                                                                             cGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGa                                          <  2:4477703/150‑1 (MQ=255)
                                                                                                                cAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGc                                       <  2:3983861/150‑1 (MQ=255)
                                                                                                                 aaaTTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCg                                      <  1:603913/150‑1 (MQ=255)
                                                                                                                  aaTTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGa                                     >  1:4910054/1‑150 (MQ=255)
                                                                                                                  aaTTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGa                                     <  2:1604703/150‑1 (MQ=255)
                                                                                                                     ttGTCGATCGCCTGCCACTGTTGAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATc                                  <  2:285023/150‑1 (MQ=255)
                                                                                                                         cGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATcagca                              <  1:870194/150‑1 (MQ=255)
                                                                                                                         cGATCGCCTGCCACTGTTTAGTGGTAATAACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATcagca                              >  2:3503880/1‑150 (MQ=255)
                                                                                                                           aTCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAgg                            >  1:397207/1‑150 (MQ=255)
                                                                                                                            tCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGa                           <  1:1787910/150‑1 (MQ=255)
                                                                                                                              gCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAg                         <  2:4672793/150‑1 (MQ=255)
                                                                                                                              gCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAg                         >  2:5066311/1‑150 (MQ=255)
                                                                                                                                cTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGtt                       >  1:3496873/1‑150 (MQ=255)
                                                                                                                                cTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGtt                       >  2:1253939/1‑150 (MQ=255)
                                                                                                                                cTGCCACTGTTTAGTGGTAATAACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGtt                       <  2:4971636/150‑1 (MQ=255)
                                                                                                                                 tGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGtcg                      <  2:3097976/150‑3 (MQ=255)
                                                                                                                                    cACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCg                   >  1:4479593/1‑150 (MQ=255)
                                                                                                                                      cTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGcc                 >  1:2561176/1‑150 (MQ=255)
                                                                                                                                       tGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGccc                >  2:118354/1‑150 (MQ=255)
                                                                                                                                         tttAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCtt              <  1:2372032/150‑1 (MQ=255)
                                                                                                                                          ttAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCTTa             >  1:3037266/1‑150 (MQ=255)
                                                                                                                                            aGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCTTaaa           >  2:2751798/1‑150 (MQ=255)
                                                                                                                                               ggTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCTTAAAgcc        <  1:285456/150‑1 (MQ=255)
                                                                                                                                                gTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCTTAAAgccg       <  1:4546548/150‑1 (MQ=255)
                                                                                                                                                 tAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCTTAAAgccgc      <  1:357532/150‑1 (MQ=255)
                                                                                                                                                   aTCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCTTAAAGccgccg    <  1:4654578/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGAGTGCAGCATCTGGTGCACCGGTTTGACGTTCTTTTTCAGAATGCCGTGGAACTGAAACGTCTGGCGGTTGGTCAGGCGAATGCTGCCATAGATGGTGTTTTCACCGGCAAATTTGTCGATCGCCTGCCACTGTTTAGTGGTAATCACCCCACCCGGCAGACGACAGCGAAGCAGCATCGCGTGGCGCGGCTCCAGCTTCTGTTCAGCACGTTCGGCGCGGATGTCGCGGTCATCCTGCTGATACATGCCGTGGAAGCGAATCAGCAGGAAGTTGTCGCCCTTAAAGCCGCCGGT  >  CP009273/2883040‑2883338

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: