Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,571,233 G→A 100% S16F (TCT→TTT)  hsdS ← specificity determinant for hsdM and hsdR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,571,2330GA98.8% 754.9 / ‑2.0 230S16F (TCT→TTT) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base A (121/106);  total (122/108)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 230 total aligned reads displayed.

TTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGC  >  CP009273/4571086‑4571378
                                                                                                                                                   |                                                                                                                                                 
  tttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATAc                                                                                                                                               >  2:1086374/1‑150 (MQ=255)
  tttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACAGTAAATAc                                                                                                                                               <  1:4656998/150‑1 (MQ=255)
    tAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTg                                                                                                                                             >  2:3286371/1‑150 (MQ=255)
    tAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTg                                                                                                                                             <  2:1243758/150‑1 (MQ=255)
       gAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTCGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGAGTAGAGTTGTGACCGTAACTACTgggg                                                                                                                                          >  2:5046017/1‑150 (MQ=255)
       gAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTgggg                                                                                                                                          >  2:1479082/1‑150 (MQ=255)
            aaaaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATa                                                                                                                                     >  2:1334767/1‑150 (MQ=255)
               aCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAacc                                                                                                                                  >  2:4158936/1‑150 (MQ=255)
                ccAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccc                                                                                                                                 >  1:586774/1‑150 (MQ=255)
                ccAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccc                                                                                                                                 >  2:834368/1‑150 (MQ=255)
                ccAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccc                                                                                                                                 >  1:664332/1‑150 (MQ=255)
                 cAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccca                                                                                                                                <  2:4896896/150‑1 (MQ=255)
                 cAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccca                                                                                                                                >  1:1816227/1‑150 (MQ=255)
                  taGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAacccac                                                                                                                               >  2:2359149/2‑150 (MQ=255)
                  aaGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAacccac                                                                                                                               <  1:2471881/150‑1 (MQ=255)
                      ccGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCt                                                                                                                           <  2:3180277/150‑1 (MQ=255)
                        gtagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTcc                                                                                                                         >  1:511847/1‑150 (MQ=255)
                         tagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCg                                                                                                                        >  1:1745982/1‑150 (MQ=255)
                         tagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCg                                                                                                                        <  2:2421752/150‑1 (MQ=255)
                         tagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCg                                                                                                                        <  1:1390390/150‑1 (MQ=255)
                         tagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCg                                                                                                                        >  1:1745901/1‑150 (MQ=255)
                         tagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTCTGACCGTAAATACTGGGGCGATAACCCACCCCTCCg                                                                                                                        >  1:4506653/1‑150 (MQ=255)
                           gtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGc                                                                                                                      <  1:4845830/150‑1 (MQ=255)
                            taTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCGTCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCa                                                                                                                     <  2:4990199/150‑1 (MQ=255)
                            taTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTCGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCa                                                                                                                     >  2:1268073/1‑150 (MQ=255)
                                aaaCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAAttt                                                                                                                 <  1:3613088/150‑1 (MQ=255)
                                 aaCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTc                                                                                                                >  2:1707985/1‑150 (MQ=255)
                                  aCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTcc                                                                                                               <  2:148768/150‑1 (MQ=255)
                                  aCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTcc                                                                                                               >  1:4049663/1‑150 (MQ=255)
                                     tGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCg                                                                                                            <  2:1800349/150‑1 (MQ=255)
                                     tGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCg                                                                                                            >  2:436927/1‑150 (MQ=255)
                                     tGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCg                                                                                                            <  2:3860601/150‑1 (MQ=255)
                                         aTTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACt                                                                                                        >  1:3248496/1‑150 (MQ=255)
                                          ttCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTc                                                                                                       <  1:162276/150‑1 (MQ=255)
                                            cTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCAt                                                                                                     <  1:1118881/150‑1 (MQ=255)
                                             tGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCAtt                                                                                                    <  2:3346381/150‑1 (MQ=255)
                                             tGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCAtt                                                                                                    <  2:2755930/150‑1 (MQ=255)
                                                  atTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTc                                                                                               <  1:3280453/150‑1 (MQ=255)
                                                  atTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTc                                                                                               >  1:2451285/1‑150 (MQ=255)
                                                   ttGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCa                                                                                              <  1:2807260/150‑1 (MQ=255)
                                                    tGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCAc                                                                                             >  2:2609126/1‑150 (MQ=255)
                                                      ttCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCAccc                                                                                           >  1:4979035/1‑150 (MQ=255)
                                                            cGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCaa                                                                                     >  2:4497583/1‑150 (MQ=255)
                                                            cGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCaa                                                                                     >  2:3617197/1‑150 (MQ=255)
                                                              gATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAAc                                                                                   <  1:4622659/150‑1 (MQ=255)
                                                                 aaGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGct                                                                                >  1:290714/1‑150 (MQ=255)
                                                                   gagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttc                                                                              <  2:4333090/150‑1 (MQ=255)
                                                                      gCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCttcttc                                                                           <  2:231067/150‑1 (MQ=255)
                                                                        aaaTAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCa                                                                         <  2:4766500/150‑1 (MQ=255)
                                                                              tcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACt                                                                   <  2:1940542/150‑1 (MQ=255)
                                                                              tcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACt                                                                   <  2:386331/150‑1 (MQ=255)
                                                                                atcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGa                                                                 <  2:644071/150‑1 (MQ=255)
                                                                                   ttttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGc                                                              >  2:3663389/1‑150 (MQ=255)
                                                                                    tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCt                                                             <  1:707904/150‑1 (MQ=255)
                                                                                    tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCt                                                             >  1:1900438/1‑150 (MQ=255)
                                                                                    tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCt                                                             >  1:1511476/1‑150 (MQ=255)
                                                                                      tAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGc                                                           >  1:3132227/1‑150 (MQ=255)
                                                                                       aGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCa                                                          >  2:657229/1‑150 (MQ=255)
                                                                                          tAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAAt                                                       <  1:1165238/150‑1 (MQ=255)
                                                                                          tAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAAt                                                       <  1:1333738/150‑1 (MQ=255)
                                                                                          tAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAAAGCTTCTTCCAGCAACTGACGCTGCAAAt                                                       >  1:1031308/1‑150 (MQ=255)
                                                                                           aaTTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCAGCAAATc                                                      >  1:2975225/1‑150 (MQ=255)
                                                                                            aTTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCg                                                     >  2:4605331/1‑150 (MQ=255)
                                                                                            aTTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCg                                                     <  2:912638/150‑1 (MQ=255)
                                                                                              ttATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCCAATCGGc                                                   >  2:1926940/1‑150 (MQ=255)
                                                                                              ttATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGc                                                   >  1:4760363/1‑150 (MQ=255)
                                                                                                  tGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCa                                               <  2:1995103/150‑1 (MQ=255)
                                                                                                  tGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCACCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCCGCAACTGACGCTGCAAATCGGCCTCa                                               >  2:4536050/1‑150 (MQ=255)
                                                                                                   gCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAt                                              <  2:1767383/150‑1 (MQ=255)
                                                                                                   gCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAt                                              <  1:889564/150‑1 (MQ=255)
                                                                                                    ccTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtc                                             <  1:1557/150‑1 (MQ=255)
                                                                                                    ccTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCCGCAAATCGGCCTCAtc                                             >  2:2802069/1‑150 (MQ=255)
                                                                                                     cTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcg                                            >  1:5010554/1‑150 (MQ=255)
                                                                                                     cTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcg                                            <  1:4999762/150‑1 (MQ=255)
                                                                                                     cTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcg                                            >  2:3002767/1‑150 (MQ=255)
                                                                                                       gCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcgct                                          >  2:3800169/1‑150 (MQ=255)
                                                                                                        ctctTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTATTCCAGCAACTGACGCTGCAAATCGGCCTCAtcgctc                                         >  1:1993868/1‑150 (MQ=255)
                                                                                                                   aCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTc                              >  2:3531086/1‑150 (MQ=255)
                                                                                                                   aCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTc                              >  2:4394160/1‑150 (MQ=255)
                                                                                                                     gTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCAc                            >  2:579781/1‑150 (MQ=255)
                                                                                                                     gTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCAc                            <  1:1350441/150‑1 (MQ=255)
                                                                                                                     gTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCACCGGCAATTTACCCCCAAGCCATACGTCACCCCACCAAAAGCATATTCCACCAACATAAGCTGCAAATAGGCCCCATCGCGCGccccccggtaaa                            >  1:2605487/1‑143 (MQ=255)
                                                                                                                     gTTACTCCTCGGATTAGAGTCGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCAc                            <  1:3228494/150‑1 (MQ=255)
                                                                                                                       tACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGc                          >  1:128626/1‑150 (MQ=255)
                                                                                                                            ctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAg                     <  2:3406782/150‑1 (MQ=255)
                                                                                                                            ctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAg                     >  1:838403/1‑150 (MQ=255)
                                                                                                                                gATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCCTTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCACATCGCTCGCCCCCAGTTCACGCATCAGCGCa                 >  1:4349149/1‑150 (MQ=255)
                                                                                                                                  ttAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAACCGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATc               <  1:3417276/150‑1 (MQ=255)
                                                                                                                                  ttAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATc               >  2:3354035/1‑150 (MQ=255)
                                                                                                                                    agagTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCCGTTCACGCATCAGCGCATCCa             >  1:2426531/1‑150 (MQ=255)
                                                                                                                                    agagTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCa             <  1:4158936/150‑1 (MQ=255)
                                                                                                                                     gagTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACGCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCCTCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAg            >  1:3048551/1‑150 (MQ=255)
                                                                                                                                      agTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGt           <  1:371198/150‑1 (MQ=255)
                                                                                                                                      agTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGt           <  1:4287944/150‑1 (MQ=255)
                                                                                                                                        ttGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTc         <  1:2812825/150‑1 (MQ=255)
                                                                                                                                        ttGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTc         <  2:552480/150‑1 (MQ=255)
                                                                                                                                         tgtgACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCa        >  2:308041/1‑150 (MQ=255)
                                                                                                                                         tgtgACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCa        <  2:2121444/150‑1 (MQ=255)
                                                                                                                                           tgACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGa      >  1:2270078/1‑150 (MQ=255)
                                                                                                                                               cGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCAACAGCGCATCCAGTTCAGACAgc  >  1:511431/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
TTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGC  >  CP009273/4571086‑4571378

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: