Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,842,161 Δ1 bp 26.6% coding (105/1044 nt) selD ← selenophosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,842,1610C.26.6% 842.8 / inf 304coding (105/1044 nt)selDselenophosphate synthase
Reads supporting (aligned to +/- strand):  ref base C (93/128);  new base . (38/43);  total (133/171)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 312 total aligned reads displayed.

GCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCC  >  CP009273/1842012‑1842310
                                                                                                                                                     |                                                                                                                                                     
gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc                                                                                                                                                       <  2:3655235/150‑1 (MQ=255)
gccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTc                                                                                                                                                       <  2:1407644/150‑1 (MQ=255)
 ccgcAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCg                                                                                                                                                      <  2:3096489/150‑1 (MQ=255)
     aaTGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTg                                                                                                                                                  >  2:2422043/1‑150 (MQ=255)
         cGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCa                                                                                                                                              <  1:680222/150‑1 (MQ=255)
          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                                                                                                                                             <  1:3226290/150‑1 (MQ=255)
          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                                                                                                                                             <  2:1254315/150‑1 (MQ=255)
          ggCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                                                                                                                                             <  1:794797/150‑1 (MQ=255)
          ggCCAAAATCGAAAGGATTATCAACGATCGACATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCAc                                                                                                                                             >  2:2327361/1‑150 (MQ=255)
           gCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACt                                                                                                                                            <  1:662425/150‑1 (MQ=255)
             cAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTAt                                                                                                                                          >  1:2120288/1‑150 (MQ=255)
                aaTCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAg                                                                                                                                      >  2:446239/1‑150 (MQ=255)
                  tCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGa                                                                                                                                    <  1:1994261/150‑1 (MQ=255)
                    gAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGAt                                                                                                                                   <  2:3090915/150‑1 (MQ=255)
                      aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATgg                                                                                                                                 >  2:206217/1‑150 (MQ=255)
                      aaGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGt                                                                                                                                <  2:3863059/150‑1 (MQ=255)
                       aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGGTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGt                                                                                                                                <  1:3126754/150‑1 (MQ=255)
                       aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGt                                                                                                                                >  1:504805/1‑150 (MQ=255)
                       aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGt                                                                                                                                >  1:3318347/1‑150 (MQ=255)
                       aGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGt                                                                                                                                >  2:1640601/1‑150 (MQ=255)
                            tATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCa                                                                                                                           <  1:3671290/150‑1 (MQ=255)
                                 aCGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACt                                                                                                                      >  2:4241953/1‑150 (MQ=255)
                                  cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCGACACtt                                                                                                                     >  2:1272806/1‑150 (MQ=255)
                                  cGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACttt                                                                                                                    <  2:2942458/150‑1 (MQ=255)
                                    aTCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACtttt                                                                                                                   >  2:936619/1‑150 (MQ=255)
                                       ggCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTggg                                                                                                                >  2:3424965/1‑150 (MQ=255)
                                         cATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGa                                                                                                              <  2:2203721/150‑1 (MQ=255)
                                           tAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGaaa                                                                                                            >  2:3017575/1‑150 (MQ=255)
                                            aaaGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAAtt                                                                                                          <  2:3289421/150‑1 (MQ=255)
                                            aaaGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAAtt                                                                                                          >  2:615567/1‑150 (MQ=255)
                                              agaagTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAAttt                                                                                                         <  2:1383257/150‑1 (MQ=255)
                                               gaagTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAAtttt                                                                                                        >  2:495461/1‑150 (MQ=255)
                                                   tCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAg                                                                                                    <  2:1826774/150‑1 (MQ=255)
                                                          tACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAAc                                                                                             >  1:4298069/1‑150 (MQ=255)
                                                           aCTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAAcc                                                                                            >  2:3858709/1‑150 (MQ=255)
                                                           aCTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAAcc                                                                                            <  2:1128434/150‑1 (MQ=255)
                                                               aTAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTc                                                                                       >  1:3990136/1‑150 (MQ=255)
                                                               aTAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTc                                                                                       >  2:2751545/1‑150 (MQ=255)
                                                               aTAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTc                                                                                       >  1:3545605/1‑150 (MQ=255)
                                                                 aaCGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTcc                                                                                      >  1:1960521/1‑150 (MQ=255)
                                                                  aCGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCg                                                                                     <  1:1891585/150‑1 (MQ=255)
                                                                    gCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTg                                                                                   <  1:4694449/150‑1 (MQ=255)
                                                                        gTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTg                                                                               <  1:1195408/150‑1 (MQ=255)
                                                                        gTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAAGTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTg                                                                               <  1:1593768/150‑1 (MQ=255)
                                                                         tGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGt                                                                              <  1:4000460/150‑1 (MQ=255)
                                                                         tGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGt                                                                              <  1:4048143/150‑1 (MQ=255)
                                                                         tGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGt                                                                              <  1:1543403/150‑1 (MQ=255)
                                                                          gCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTAt                                                                            >  1:2265130/1‑150 (MQ=255)
                                                                          gCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTAt                                                                            >  2:740028/1‑150 (MQ=255)
                                                                          gCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTAt                                                                            >  1:3451813/1‑150 (MQ=255)
                                                                          gCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTAt                                                                            <  2:1291424/150‑1 (MQ=255)
                                                                           ccATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTAt                                                                            <  1:703108/150‑1 (MQ=255)
                                                                            cATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTAtt                                                                           <  2:2110516/150‑1 (MQ=255)
                                                                               tGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTggg                                                                        >  1:2384963/1‑150 (MQ=255)
                                                                               tGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGt                                                                       >  2:4116644/1‑150 (MQ=255)
                                                                                 cccAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCa                                                                     >  2:3976337/1‑150 (MQ=255)
                                                                                 cccAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCa                                                                     >  2:1516838/1‑150 (MQ=255)
                                                                                  ccAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCa                                                                     <  1:1018525/150‑1 (MQ=255)
                                                                                    aGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCaaa                                                                   <  1:4497371/150‑1 (MQ=255)
                                                                                     gATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCaaaa                                                                  <  1:905712/150‑2 (MQ=255)
                                                                                            cacCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGa                                                           >  2:1859212/1‑150 (MQ=255)
                                                                                             acCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGt                                                         <  1:4383807/150‑1 (MQ=255)
                                                                                              ccgccgCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGt                                                         <  2:3153904/150‑1 (MQ=255)
                                                                                                 ccgcATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTtct                                                      <  2:2385049/150‑1 (MQ=255)
                                                                                                 ccgcATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTtctc                                                     >  1:1724164/1‑150 (MQ=255)
                                                                                                     aTCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCt                                                  <  2:4381215/150‑1 (MQ=255)
                                                                                                     aTCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCt                                                  >  1:3483350/1‑150 (MQ=255)
                                                                                                     aTCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCt                                                  >  2:1567963/1‑150 (MQ=255)
                                                                                                      tcgtcgCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCa                                                <  2:926976/150‑1 (MQ=255)
                                                                                                         tcgCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATgg                                             <  2:1614502/150‑1 (MQ=255)
                                                                                                             gggTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACa                                          <  2:2767995/150‑1 (MQ=255)
                                                                                                               gTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACAtc                                        >  2:2656160/1‑150 (MQ=255)
                                                                                                                tttCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCGACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACAtct                                       <  1:1272806/150‑1 (MQ=255)
                                                                                                                   cATTACCCACAAGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTg                                   >  2:2632675/1‑150 (MQ=255)
                                                                                                                             aaGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCaa                          <  2:1846251/150‑1 (MQ=255)
                                                                                                                              aGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATc                        <  2:4548392/150‑1 (MQ=255)
                                                                                                                              aGCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATc                        <  2:3284540/150‑1 (MQ=255)
                                                                                                                               gCAAATTCGGATCAACAAACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATcc                       >  1:3883880/1‑150 (MQ=255)
                                                                                                                                 aaaTTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTTGGGAAATTTTACAGCCTCAACCAGCTCCGTGGCTGTTTTGGTTCAAACGAATCGATTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCg                      <  1:1968289/150‑1 (MQ=255)
                                                                                                                                 aaaTTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCg                      <  1:275260/150‑1 (MQ=255)
                                                                                                                                   aTTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCggg                    <  2:3474720/150‑1 (MQ=255)
                                                                                                                                    ttCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGt                   >  2:1228382/1‑150 (MQ=255)
                                                                                                                                      cGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTAt                 >  2:4026905/1‑150 (MQ=255)
                                                                                                                                        gATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATgg               <  2:2684042/150‑1 (MQ=255)
                                                                                                                                        gATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATgg               <  1:2525619/150‑1 (MQ=255)
                                                                                                                                          tCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTa             <  2:4640830/150‑1 (MQ=255)
                                                                                                                                            aacaaACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAAc           <  2:2522005/150‑1 (MQ=255)
                                                                                                                                            aacaaACTT‑GCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAAcc          <  1:554957/150‑1 (MQ=255)
                                                                                                                                             acaaACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAAcc          <  1:1916364/150‑1 (MQ=255)
                                                                                                                                             acaaACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAAcc          <  1:1067898/150‑1 (MQ=255)
                                                                                                                                              caaACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAAccc         <  1:4625092/150‑1 (MQ=255)
                                                                                                                                              caaACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAAccc         <  1:2471028/150‑1 (MQ=255)
                                                                                                                                               aaaCTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCt        <  2:1976012/150‑1 (MQ=255)
                                                                                                                                                aaCTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTc       >  2:574574/1‑150 (MQ=255)
                                                                                                                                                aaCTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTc       <  1:2648033/150‑1 (MQ=255)
                                                                                                                                                aaCTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTc       >  2:1920856/1‑150 (MQ=255)
                                                                                                                                                 aCTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCa      <  1:4092695/150‑1 (MQ=255)
                                                                                                                                                  cTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCAt     >  2:4164185/1‑150 (MQ=255)
                                                                                                                                                   ttCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCAtt    <  2:4012963/150‑1 (MQ=255)
                                                                                                                                                   ttCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCAtt    <  2:3603415/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGCAATGCGGCCAAAATCGAAAGGATTATCAACGATCGGCATAAAGAAGTCGGTGGTACTGATAACGCTGGTGCCATTGCCCAGATCGTACACCGCCGCATCGTCGCGGGTTTCATTACCCACAAGCAAATTCGGATCAACAAACTTCGCCTGCTCACTATGCAGGATGGTTTCCAACACTTTTGGGGAAATTTTACAGCCGCAACCAGCTCCGTGGCTGTATTGGGTCAAACGAATCGAGTTCTCGCTCATGGACATCTCCTGTCAATGCAATCCGGGTATGGTAACCCTCATTCC  >  CP009273/1842012‑1842310

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