Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,987 G→T 74.9% S271* (TCG→TAG)  gshA ← glutamate‑cysteine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,808,9870GT74.9% 410.0 / 233.5 271S271* (TCG→TAG) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (41/27);  new base T (99/104);  total (140/131)
Fisher's exact test for biased strand distribution p-value = 1.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 271 total aligned reads displayed.

AACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTC  >  CP009273/2808840‑2809135
                                                                                                                                                   |                                                                                                                                                    
aaCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTAt                                                                                                                                                    <  2:954709/150‑1 (MQ=255)
 aCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATg                                                                                                                                                   >  1:1462495/1‑150 (MQ=255)
  cAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATgg                                                                                                                                                  <  1:2647266/150‑1 (MQ=255)
   aGCGCATCCGCAGGCGACTCGCCGCTGCGGGGAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGc                                                                                                                                                 <  2:2567033/150‑1 (MQ=255)
    gcgcATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCg                                                                                                                                                >  2:2164582/1‑150 (MQ=255)
            gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGTATCTGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTTAtt                                                                                                                                        <  1:4143846/150‑1 (MQ=255)
             aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTg                                                                                                                                       >  2:4133814/1‑150 (MQ=255)
             aaGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTg                                                                                                                                       >  2:3438117/1‑150 (MQ=255)
               ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGcc                                                                                                                                     >  1:3877658/1‑150 (MQ=255)
               ggCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGcc                                                                                                                                     >  2:3956879/1‑150 (MQ=255)
                gCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCt                                                                                                                                    <  2:2415073/150‑1 (MQ=255)
                gCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCt                                                                                                                                    >  2:1254433/1‑150 (MQ=255)
                gCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCt                                                                                                                                    <  2:2415166/150‑1 (MQ=255)
                 cGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTg                                                                                                                                   <  1:3133616/150‑1 (MQ=255)
                 cGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTg                                                                                                                                   >  2:2929572/1‑150 (MQ=255)
                  gACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGt                                                                                                                                  >  2:282157/1‑150 (MQ=255)
                  gACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGt                                                                                                                                  >  2:3636543/1‑150 (MQ=255)
                    cTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGttt                                                                                                                                <  1:2227338/150‑1 (MQ=255)
                     tCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTa                                                                                                                               >  1:4071162/1‑150 (MQ=255)
                        ccgcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAgg                                                                                                                            <  1:1300434/150‑1 (MQ=255)
                         cgcTGCGGGTAACGAGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGc                                                                                                                           <  1:3545521/150‑1 (MQ=255)
                          gcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGcc                                                                                                                          <  2:2873729/150‑1 (MQ=255)
                            tGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCg                                                                                                                        <  2:2093480/150‑1 (MQ=255)
                              cGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCt                                                                                                                      >  1:2857751/1‑150 (MQ=255)
                                ggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGAGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCTCTAc                                                                                                                    <  1:65086/150‑1 (MQ=255)
                                   aaCGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTGTTGATTGCCTGTTTAAGGCCCGCTACGTa                                                                                                                 >  1:4038659/1‑150 (MQ=255)
                                   aaCGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTa                                                                                                                 <  2:3952883/150‑1 (MQ=255)
                                     cgcgTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACt                                                                                                               <  1:318978/150‑1 (MQ=255)
                                      gcgTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTc                                                                                                              <  2:606067/150‑1 (MQ=255)
                                          tttGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaa                                                                                                          <  1:2954232/150‑1 (MQ=255)
                                           ttGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTaaa                                                                                                         <  1:3581317/150‑1 (MQ=255)
                                             ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                                                                                                       >  2:489698/1‑150 (MQ=255)
                                             ggACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGa                                                                                                       <  2:3848062/150‑1 (MQ=255)
                                                 gAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGt                                                                                                   >  1:3418695/1‑150 (MQ=255)
                                                 gAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGt                                                                                                   >  1:364912/1‑150 (MQ=255)
                                                  aaTCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGtt                                                                                                  <  1:2134077/150‑1 (MQ=255)
                                                      ggCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAg                                                                                              <  2:1588019/150‑1 (MQ=255)
                                                         gcgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTa                                                                                           >  2:1424921/1‑150 (MQ=255)
                                                          cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa                                                                                          <  2:1089522/150‑1 (MQ=255)
                                                             aCAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCCTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATAc                                                                                       >  2:3066050/1‑150 (MQ=255)
                                                               aGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCa                                                                                     <  1:3802812/150‑1 (MQ=255)
                                                                 ttCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAg                                                                                   >  1:4212066/1‑150 (MQ=255)
                                                                 ttCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAg                                                                                   >  1:4212575/1‑150 (MQ=255)
                                                                   cGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAt                                                                                 >  1:2319574/1‑150 (MQ=255)
                                                                    gTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAtt                                                                                <  1:2615304/150‑1 (MQ=255)
                                                                     atttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCGTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTg                                                                               <  2:1634277/149‑1 (MQ=255)
                                                                      tttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGc                                                                              >  1:3059291/1‑150 (MQ=255)
                                                                        tCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCtt                                                                            <  1:1107226/150‑1 (MQ=255)
                                                                         cAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCttt                                                                           <  2:4345925/150‑1 (MQ=255)
                                                                         cAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCttt                                                                           >  2:4576896/1‑150 (MQ=255)
                                                                          aaTCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTg                                                                          >  1:1517233/1‑150 (MQ=255)
                                                                          aaTCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTg                                                                          >  1:1723306/1‑150 (MQ=255)
                                                                              tGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAt                                                                      >  2:2721285/1‑150 (MQ=255)
                                                                               gCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAtt                                                                     >  2:3518781/1‑150 (MQ=255)
                                                                               gCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAtt                                                                     <  1:1272063/150‑1 (MQ=255)
                                                                                cAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAttt                                                                    <  1:3767194/150‑1 (MQ=255)
                                                                                cAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAttt                                                                    >  2:3134159/1‑150 (MQ=255)
                                                                                cAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGAttt                                                                    >  2:4116805/1‑150 (MQ=255)
                                                                                 aaCACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                   >  1:1044340/1‑150 (MQ=255)
                                                                                   cacGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTAtt                                                                 >  2:4106220/1‑150 (MQ=255)
                                                                                    acGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTg                                                                <  2:4499369/150‑1 (MQ=255)
                                                                                      gTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGt                                                              >  2:3996846/1‑150 (MQ=255)
                                                                                      gTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGt                                                              <  1:1820083/150‑1 (MQ=255)
                                                                                      gTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGt                                                              <  1:1213832/150‑1 (MQ=255)
                                                                                      gTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGt                                                              >  1:2427999/1‑150 (MQ=255)
                                                                                      gTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGt                                                              >  1:1982170/1‑150 (MQ=255)
                                                                                         gCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGtata                                                           <  1:756888/150‑1 (MQ=255)
                                                                                           tgatGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTAAAGGCCCGCAACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                                                         >  1:2049364/4‑150 (MQ=255)
                                                                                           tGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                                                         >  2:4052941/1‑150 (MQ=255)
                                                                                           tGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                                                         >  1:4514118/1‑150 (MQ=255)
                                                                                           tGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGc                                                         >  1:4408507/1‑150 (MQ=255)
                                                                                            gTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGcc                                                        <  1:2929572/150‑1 (MQ=255)
                                                                                            gTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGcc                                                        >  2:3181756/1‑150 (MQ=255)
                                                                                             ttGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCg                                                       <  1:476257/150‑1 (MQ=255)
                                                                                             ttGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCg                                                       <  2:4300288/150‑1 (MQ=255)
                                                                                               gATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCgag                                                     <  2:406317/150‑1 (MQ=255)
                                                                                                aTTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCgaga                                                    >  2:2652483/1‑150 (MQ=255)
                                                                                                      aGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTc                                              <  2:1914643/150‑1 (MQ=255)
                                                                                                           ctTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACg                                         >  2:1116575/1‑150 (MQ=255)
                                                                                                            ttACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGa                                        >  1:907020/1‑150 (MQ=255)
                                                                                                              aCCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAag                                      >  2:2442870/1‑150 (MQ=255)
                                                                                                               ccGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAaga                                     >  2:3720351/1‑150 (MQ=255)
                                                                                                                cGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagag                                    >  1:3286773/1‑150 (MQ=255)
                                                                                                                 gTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagaga                                   >  2:1307000/1‑150 (MQ=255)
                                                                                                                     ttCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTc                               >  2:1164094/1‑150 (MQ=255)
                                                                                                                      tctcAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCg                              <  2:2296987/150‑1 (MQ=255)
                                                                                                                       ctcAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGc                             >  2:4438032/1‑150 (MQ=255)
                                                                                                                            tACCACTCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACg                        <  1:1163233/150‑1 (MQ=255)
                                                                                                                              ccAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGc                      >  1:1191131/1‑150 (MQ=255)
                                                                                                                               cAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCa                     <  1:558201/150‑1 (MQ=255)
                                                                                                                                    ttCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTaa                <  2:701270/150‑1 (MQ=255)
                                                                                                                                        cgTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAataca            >  1:4355157/1‑150 (MQ=255)
                                                                                                                                        cgTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAataca            >  2:4064720/1‑150 (MQ=255)
                                                                                                                                         gTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacat           <  1:282157/150‑1 (MQ=255)
                                                                                                                                         gTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacat           >  2:1798457/1‑150 (MQ=255)
                                                                                                                                          tACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacata          <  1:1115621/150‑1 (MQ=255)
                                                                                                                                           aCTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacatac         >  1:2570515/1‑150 (MQ=255)
                                                                                                                                           aCTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacatac         <  1:1511051/150‑1 (MQ=255)
                                                                                                                                            ctctTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATAcc        <  2:3009936/150‑1 (MQ=255)
                                                                                                                                                tCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCAc    <  2:1805827/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTC  >  CP009273/2808840‑2809135

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: