Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,865,863 T→A 7.0% T166S (ACC→TCC)  dgoD ← galactonate dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,865,8630TA 7.0% 899.5 / 43.9 272T166S (ACC→TCC) dgoDgalactonate dehydratase
Reads supporting (aligned to +/- strand):  ref base T (104/149);  new base A (11/8);  total (115/157)
Fisher's exact test for biased strand distribution p-value = 2.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 283 total aligned reads displayed.

TCCGCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTGTAGG  >  CP009273/3865715‑3866012
                                                                                                                                                    |                                                                                                                                                     
tCCGCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGgtg                                                                                                                                                      <  2:2630769/150‑1 (MQ=255)
 ccGCCAGCACCGGCTCCTCAATACACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGgagt                                                                                                                                                     <  1:2100220/150‑1 (MQ=255)
 ccGCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGgtgt                                                                                                                                                     <  1:3553139/150‑1 (MQ=255)
   gCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTa                                                                                                                                                   >  2:864846/1‑150 (MQ=255)
    ccAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTaa                                                                                                                                                  <  1:3479964/150‑1 (MQ=255)
    ccAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTaa                                                                                                                                                  >  1:3960756/1‑150 (MQ=255)
     cAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAAc                                                                                                                                                 <  1:4522015/150‑1 (MQ=255)
      aGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAAcc                                                                                                                                                <  1:810512/150‑1 (MQ=255)
        cACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAAccgc                                                                                                                                              >  2:154630/1‑150 (MQ=255)
           cGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACcgccgc                                                                                                                                           >  1:2582374/1‑150 (MQ=255)
            ggCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCa                                                                                                                                          <  2:3423924/150‑1 (MQ=255)
            ggCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCa                                                                                                                                          <  2:4272088/150‑1 (MQ=255)
               tcctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCt                                                                                                                                       <  2:3298506/150‑1 (MQ=255)
               tcctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCt                                                                                                                                       <  2:2275643/150‑1 (MQ=255)
                cctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTa                                                                                                                                      <  1:3333897/150‑1 (MQ=255)
                 ctcAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTAc                                                                                                                                     <  2:3559556/150‑1 (MQ=255)
                    aaTAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACcgc                                                                                                                                  <  2:3935817/150‑1 (MQ=255)
                      tAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACcgcgc                                                                                                                                >  2:1027853/1‑150 (MQ=255)
                         aCAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCggg                                                                                                                             >  2:49167/1‑150 (MQ=255)
                          cAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGa                                                                                                                            <  2:3243643/150‑1 (MQ=255)
                          cAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGa                                                                                                                            <  2:2661138/150‑1 (MQ=255)
                              gggCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTa                                                                                                                        >  2:4128565/1‑150 (MQ=255)
                              gggCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTa                                                                                                                        >  2:2965858/1‑150 (MQ=255)
                                   aTACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAAt                                                                                                                   >  1:3296685/1‑150 (MQ=255)
                                    tACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAAtt                                                                                                                  <  1:494571/150‑1 (MQ=255)
                                    tACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAAtt                                                                                                                  <  1:3180542/150‑1 (MQ=255)
                                      cGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAg                                                                                                                >  1:3244529/1‑150 (MQ=255)
                                      cGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAg                                                                                                                <  1:2610827/150‑1 (MQ=255)
                                        gCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGcc                                                                                                              <  2:1445056/150‑1 (MQ=255)
                                            cAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAg                                                                                                          <  1:2199381/150‑1 (MQ=255)
                                             aGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGt                                                                                                         <  2:158605/150‑1 (MQ=255)
                                              gTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGtt                                                                                                        <  2:2263487/150‑1 (MQ=255)
                                                 cTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGttctt                                                                                                     <  1:446432/150‑1 (MQ=255)
                                                       tCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTATGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAAc                                                                                               >  1:1863799/1‑150 (MQ=255)
                                                          gCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGt                                                                                            <  1:4506054/150‑1 (MQ=255)
                                                           cACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGtt                                                                                           <  1:3641272/150‑1 (MQ=255)
                                                                 cGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGttt                                                                                     >  1:707644/1‑150 (MQ=255)
                                                                  gCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATCAGCCCCAGTTCTTCACAACCGTTCAGTTTg                                                                                    <  1:2378957/150‑1 (MQ=255)
                                                                   ccATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGa                                                                                   >  2:4502717/1‑150 (MQ=255)
                                                                    cATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGaa                                                                                  >  2:2894578/1‑150 (MQ=255)
                                                                     aTCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAg                                                                                 >  2:2035/1‑150 (MQ=255)
                                                                     aTCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAg                                                                                 >  1:4085129/1‑150 (MQ=255)
                                                                        ggCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTa                                                                              >  1:1820946/1‑150 (MQ=255)
                                                                         gcgcgcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTAt                                                                             <  1:1047185/150‑1 (MQ=255)
                                                                          cgcgcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATc                                                                            <  2:4598218/150‑1 (MQ=255)
                                                                          cgcgcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATc                                                                            <  2:1138946/150‑1 (MQ=255)
                                                                             gcTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGaa                                                                         <  2:4125826/150‑1 (MQ=255)
                                                                              cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAg                                                                        <  1:597351/150‑1 (MQ=255)
                                                                              cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAg                                                                        >  2:1804758/1‑150 (MQ=255)
                                                                              cTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAg                                                                        <  2:582519/150‑1 (MQ=255)
                                                                                gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGcc                                                                      <  2:3643559/150‑1 (MQ=255)
                                                                                gACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGcc                                                                      <  2:1324191/150‑1 (MQ=255)
                                                                                  cgcgACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGa                                                                    <  1:1645762/150‑1 (MQ=255)
                                                                                     gACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAttt                                                                 >  2:3137058/1‑150 (MQ=255)
                                                                                     gACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGAttt                                                                 <  2:261754/150‑1 (MQ=255)
                                                                                      aCCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTc                                                                <  2:3222740/150‑1 (MQ=255)
                                                                                       ccGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTcg                                                               <  1:3074273/150‑1 (MQ=255)
                                                                                        cGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTcgc                                                              <  2:1742630/150‑1 (MQ=255)
                                                                                           ggAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCGCGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTa                                                           >  1:2902349/1‑150 (MQ=255)
                                                                                           ggAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCGCGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTa                                                           >  1:2847135/1‑150 (MQ=255)
                                                                                           ggAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTa                                                           >  1:2971629/1‑150 (MQ=255)
                                                                                            gAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAg                                                          >  1:1670926/1‑150 (MQ=255)
                                                                                             aaaTCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGc                                                         <  1:4243746/150‑1 (MQ=255)
                                                                                              aaTCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCg                                                        >  2:3754767/1‑150 (MQ=255)
                                                                                              aaTCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCg                                                        <  1:2317920/150‑1 (MQ=255)
                                                                                              aaTCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCg                                                        >  2:20935/1‑150 (MQ=255)
                                                                                               aTCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGt                                                       >  2:913177/1‑150 (MQ=255)
                                                                                                tCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGtt                                                      >  2:3061154/1‑150 (MQ=255)
                                                                                                 cAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGttt                                                     <  2:2929733/150‑1 (MQ=255)
                                                                                                      ccAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGAGTTTTAATg                                                >  2:2846266/1‑150 (MQ=255)
                                                                                                      ccAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTAGCGTAGCGTTTTAATg                                                >  1:421448/1‑150 (MQ=255)
                                                                                                        aaaCTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGcc                                              <  2:369814/150‑1 (MQ=255)
                                                                                                        aaaCTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGcc                                              <  1:2126500/150‑1 (MQ=255)
                                                                                                         aaCTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCg                                             >  2:3999403/1‑150 (MQ=255)
                                                                                                          aCTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGt                                            <  1:896463/150‑1 (MQ=255)
                                                                                                           cTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTc                                           <  1:4398184/150‑1 (MQ=255)
                                                                                                                tCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATaa                                      <  2:1391810/150‑1 (MQ=255)
                                                                                                                  tGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACa                                    <  2:887205/150‑1 (MQ=255)
                                                                                                                    aTTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATc                                  <  1:1504148/150‑1 (MQ=255)
                                                                                                                     ttGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATcc                                 <  1:2590879/150‑1 (MQ=255)
                                                                                                                      tGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccg                                <  1:2512052/150‑1 (MQ=255)
                                                                                                                       gCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccgc                               <  2:242214/150‑1 (MQ=255)
                                                                                                                        ccAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATccgcc                              >  1:3731382/1‑150 (MQ=255)
                                                                                                                          aaaaGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCgg                            <  2:3881163/150‑1 (MQ=255)
                                                                                                                            aaGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGAc                          <  1:1307354/150‑1 (MQ=255)
                                                                                                                                ttCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATc                      <  1:336824/150‑1 (MQ=255)
                                                                                                                                ttCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATc                      <  2:1167711/150‑1 (MQ=255)
                                                                                                                                 tCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGATTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcg                     <  2:4401953/150‑1 (MQ=255)
                                                                                                                                   aCGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgcc                   <  2:1683409/150‑1 (MQ=255)
                                                                                                                                    cGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGCAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccg                  >  2:1704889/1‑150 (MQ=255)
                                                                                                                                     gAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccgc                 >  1:2025631/1‑150 (MQ=255)
                                                                                                                                     gAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccgc                 >  2:1900298/1‑150 (MQ=255)
                                                                                                                                     gAATTTGTGCCACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATcgccgc                 >  1:4267119/1‑150 (MQ=255)
                                                                                                                                      aaTTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCgccgcc                >  2:3456518/1‑150 (MQ=255)
                                                                                                                                       aTTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGcca               <  1:305638/150‑1 (MQ=255)
                                                                                                                                       aTTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGcca               <  1:1408032/150‑1 (MQ=255)
                                                                                                                                             gCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaacccaa         >  1:3458489/1‑150 (MQ=255)
                                                                                                                                              ccACGGAGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGccaacccaac        >  1:953422/1‑150 (MQ=255)
                                                                                                                                                   gtgtTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTGTAg   <  1:1726524/150‑1 (MQ=255)
                                                                                                                                                    tgtTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTGTAgg  <  2:4381802/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TCCGCCAGCACCGGCTCCTCAATAAACAGCGGGCGATACGGCTCCAGTTCTTTAATCAGCACTTTCGCCATCGGCGCGCTGACGCGACCGTGGAAATCAAGACCAAACTCAATCTGATTGCCAAAAGCTTCACGAATTTGTGCCACGGTGTTAACCGCCGCATCTACCGCGCGGGAGTTATCAATTAGCCCCAGTTCTTCACAACCGTTCAGTTTGAAGGTATCGAAGCCGATTTCGCGTAGCGTTTTAATGCCGTCGATAACATCCGCCGGACGATCGCCGCCAACCCAACTGTAGG  >  CP009273/3865715‑3866012

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: