Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,808,987 G→T 68.7% 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,9870GT68.7% 280.5 / inf 285S271* (TCG→TAG) gshAglutamate‑cysteine ligase
Reads supporting (aligned to +/- strand):  ref base G (40/49);  new base T (105/91);  total (145/140)
Fisher's exact test for biased strand distribution p-value = 2.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

Only 100 of 287 total aligned reads displayed.

GTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTC  >  CP009273/2808838‑2809135
                                                                                                                                                     |                                                                                                                                                    
  aaCAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTAt                                                                                                                                                    >  1:3162153/1‑150 (MQ=255)
         cATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGttt                                                                                                                                             >  1:1570390/1‑150 (MQ=255)
             aGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGAt                                                                                                                                         <  1:4486698/150‑1 (MQ=255)
              gAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGAtt                                                                                                                                        <  1:3238747/150‑1 (MQ=255)
                   cGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTg                                                                                                                                   <  2:4225229/150‑1 (MQ=255)
                   cGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTg                                                                                                                                   <  1:4397212/150‑1 (MQ=255)
                      cTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGttt                                                                                                                                >  1:1302628/1‑150 (MQ=255)
                       tCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTa                                                                                                                               >  1:2365003/1‑150 (MQ=255)
                       tCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTa                                                                                                                               >  1:4059743/1‑150 (MQ=255)
                            gcTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGcc                                                                                                                          <  1:768165/150‑1 (MQ=255)
                             cTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGccc                                                                                                                         >  2:2271714/1‑150 (MQ=255)
                                 gggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTa                                                                                                                     <  1:2483602/150‑1 (MQ=255)
                                 gggTAACGCGTTTTGGACGAATCGGCGCGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTa                                                                                                                     <  1:1125241/150‑1 (MQ=255)
                                    tAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGt                                                                                                                  >  2:2506716/1‑150 (MQ=255)
                                     aaCGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTGTTGATTGCCTGTTTAAGGCCCGCTACGTa                                                                                                                 >  2:2641527/1‑150 (MQ=255)
                                     aaCGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTa                                                                                                                 >  2:1099516/1‑150 (MQ=255)
                                                gACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGAt                                                                                                      <  1:1671805/150‑1 (MQ=255)
                                                 aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATc                                                                                                     <  2:925944/150‑1 (MQ=255)
                                                 aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATc                                                                                                     <  2:530368/150‑1 (MQ=255)
                                                 aCGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATc                                                                                                     <  2:1575587/150‑1 (MQ=255)
                                                     aTCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTg                                                                                                 <  2:4435935/150‑1 (MQ=255)
                                                        ggCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAg                                                                                              <  1:3027085/150‑1 (MQ=255)
                                                         gcgcGTACAGTTAGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAgg                                                                                             <  1:2757937/150‑1 (MQ=255)
                                                          cgcgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAATGt                                                                                            >  1:291286/1‑150 (MQ=255)
                                                            cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa                                                                                          >  1:1052701/1‑150 (MQ=255)
                                                            cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa                                                                                          <  2:3087789/150‑1 (MQ=255)
                                                            cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa                                                                                          <  1:1263102/150‑1 (MQ=255)
                                                            cgTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTaa                                                                                          <  1:3299491/150‑1 (MQ=255)
                                                             gTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt                                                                                         <  1:1464924/150‑1 (MQ=255)
                                                             gTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAAt                                                                                         <  1:528356/150‑1 (MQ=255)
                                                                 aGTTCGTTTTCAATCTGCAACACGGTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCa                                                                                     <  2:1777558/150‑1 (MQ=255)
                                                                   ttCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAg                                                                                   <  2:4525596/150‑1 (MQ=255)
                                                                    tCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGa                                                                                  <  2:3293218/150‑1 (MQ=255)
                                                                     cGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAt                                                                                 >  2:1378979/1‑150 (MQ=255)
                                                                      gTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGAtt                                                                                <  1:1986794/150‑1 (MQ=255)
                                                                        tttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGc                                                                              <  2:678900/150‑1 (MQ=255)
                                                                        tttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGc                                                                              <  2:4573563/150‑1 (MQ=255)
                                                                        tttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGc                                                                              <  2:3609900/150‑1 (MQ=255)
                                                                         ttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCt                                                                             >  1:840650/1‑150 (MQ=255)
                                                                         ttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCt                                                                             <  2:3951067/150‑1 (MQ=255)
                                                                         ttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCt                                                                             >  1:2141173/1‑150 (MQ=255)
                                                                         ttCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCATCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGAAATACCAAGATTGCt                                                                             >  1:4176127/1‑150 (MQ=255)
                                                                          tCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCtt                                                                            >  2:2104975/1‑150 (MQ=255)
                                                                           cAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCttt                                                                           >  2:1471418/1‑150 (MQ=255)
                                                                             aTCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGc                                                                         >  2:2936604/1‑150 (MQ=255)
                                                                               cTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGa                                                                       <  1:353382/150‑1 (MQ=255)
                                                                                   aaCACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                   >  1:3335819/1‑150 (MQ=255)
                                                                                   aaCACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                   >  2:62777/1‑150 (MQ=255)
                                                                                   aaCACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                   >  1:4013688/1‑150 (MQ=255)
                                                                                   aaCACGTTGCTGTTGATTTGCAGCCTCTTACAGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTa                                                                   >  2:4625264/1‑150 (MQ=255)
                                                                                    acacGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTAt                                                                  <  2:3255980/150‑1 (MQ=255)
                                                                                      acGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTg                                                                >  2:3582945/1‑150 (MQ=255)
                                                                                      acGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTg                                                                <  2:3429394/150‑1 (MQ=255)
                                                                                        gTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGt                                                              >  2:78130/1‑150 (MQ=255)
                                                                                               ttGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCg                                                       >  2:2758722/1‑150 (MQ=255)
                                                                                                tGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCga                                                      >  2:1259261/1‑150 (MQ=255)
                                                                                                   tttGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGAt                                                   <  2:2457276/150‑1 (MQ=255)
                                                                                                       cAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCt                                               <  1:3620157/150‑1 (MQ=255)
                                                                                                           ctctTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCaaa                                           >  1:267341/1‑150 (MQ=255)
                                                                                                            tctTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAAc                                          <  1:651216/150‑1 (MQ=255)
                                                                                                              ttACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGGTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGa                                        >  2:140519/1‑150 (MQ=255)
                                                                                                              ttACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGa                                        >  1:3938429/1‑150 (MQ=255)
                                                                                                              ttACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGa                                        >  2:1076922/1‑150 (MQ=255)
                                                                                                              ttACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGa                                        >  2:637110/1‑150 (MQ=255)
                                                                                                               tACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGaa                                       >  2:2046021/1‑150 (MQ=255)
                                                                                                               tACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGaa                                       <  1:792507/150‑1 (MQ=255)
                                                                                                                  cGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagag                                    <  1:4410925/150‑1 (MQ=255)
                                                                                                                  cGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagag                                    >  1:3090888/1‑150 (MQ=255)
                                                                                                                  cGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAagag                                    >  1:4266761/1‑150 (MQ=255)
                                                                                                                       ttCTCAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTc                               >  1:1402586/1‑150 (MQ=255)
                                                                                                                       ttCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTc                               <  2:1905127/150‑1 (MQ=255)
                                                                                                                       ttCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTc                               <  2:1691067/150‑1 (MQ=255)
                                                                                                                        tctcAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCg                              <  2:4143172/150‑1 (MQ=255)
                                                                                                                        tctcAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCg                              <  2:1003702/150‑1 (MQ=255)
                                                                                                                         ctcAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGc                             <  1:1810042/150‑1 (MQ=255)
                                                                                                                           cAATACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCAt                           >  1:1884506/1‑150 (MQ=255)
                                                                                                                           cAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCACGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCAt                           <  2:1552446/150‑1 (MQ=255)
                                                                                                                            aaTACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATa                          >  1:2565661/1‑150 (MQ=255)
                                                                                                                             aTACCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATAc                         >  2:3816025/1‑150 (MQ=255)
                                                                                                                               aCCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACgg                       >  2:4578672/1‑150 (MQ=255)
                                                                                                                               aCCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACgg                       <  2:3747621/150‑1 (MQ=255)
                                                                                                                               aCCAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACgg                       >  2:4260942/1‑150 (MQ=255)
                                                                                                                                ccAATCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGc                      >  2:221463/1‑150 (MQ=255)
                                                                                                                                  aaTCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAg                    <  2:3520034/150‑1 (MQ=255)
                                                                                                                                  aaTCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAg                    >  2:2789344/1‑150 (MQ=255)
                                                                                                                                  aaTCTTCGCGTACTCTTCCGATGGCGTGTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAg                    >  1:1036986/1‑150 (MQ=255)
                                                                                                                                    tCTTCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGt                  <  1:586723/150‑1 (MQ=255)
                                                                                                                                     cTTCGCGTACTCTTCCTCTGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTa                 <  2:31011/150‑1 (MQ=255)
                                                                                                                                      ttCGCGTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTaa                <  2:1294365/150‑1 (MQ=255)
                                                                                                                                       tCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAat               <  2:2512474/150‑1 (MQ=255)
                                                                                                                                        cgcgTACTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAata              <  2:2667280/150‑1 (MQ=255)
                                                                                                                                           gTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacat           <  1:3751490/150‑1 (MQ=255)
                                                                                                                                             aCTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacatac         <  1:3948667/150‑1 (MQ=255)
                                                                                                                                             aCTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacatac         <  2:4357957/150‑1 (MQ=255)
                                                                                                                                             aCTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacatac         <  2:226630/150‑1 (MQ=255)
                                                                                                                                             aCTCTTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAatacatac         >  1:343756/1‑150 (MQ=255)
                                                                                                                                              ctctTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATAcc        <  2:3194084/150‑1 (MQ=255)
                                                                                                                                               tctTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCg       >  1:3383940/1‑150 (MQ=255)
                                                                                                                                               tctTCCTATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCg       <  1:4321504/150‑1 (MQ=255)
                                                                                                                                                   ccGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACt   >  2:4150868/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTAACAGCGCATCAGAAGGCGACTCGCCGCTGCGGGTAACGCGTTTTGGACGAATCGGCGCGTACAGTTCGTTTTCAATCTGCAACACGTTGCTGTTGATTTGCAGCCTCTTACCGTCTTTCTCAATACCAATCTTCGCGTACTCTTCCGATGGCGTTTTGATTGCCTGTTTAAGGCCCGCTACGTACTCGTAAAGATCGTTGAAGGTAATACCAAGATTGCTTTGCGATTTATTGGTATAGCCGAGATCGCTCAAACGAAGAGAGGTCGCATACGGCAGGTAATACATACCGCACTC  >  CP009273/2808838‑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: