Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,571,233 G→A 45.6% 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,2330GA45.6% 118.3 / inf 272S16F (TCT→TTT) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (100/48);  new base A (74/50);  total (174/98)
Fisher's exact test for biased strand distribution p-value = 2.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 273 total aligned reads displayed.

ATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCC  >  CP009273/4571085‑4571381
                                                                                                                                                    |                                                                                                                                                    
 tttttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGAt                                                                                                                                                    <  2:1791320/150‑1 (MQ=255)
 tttttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAAt                                                                                                                                                    <  2:2352575/150‑1 (MQ=255)
 tttttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAAt                                                                                                                                                    <  2:2353442/150‑1 (MQ=255)
     tAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTg                                                                                                                                                <  1:4341816/150‑1 (MQ=255)
         aaCAAAAACCAAGTCCGTAGTTTCAATCTTGCCATTCTGAATTTTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGTTTAGAGTTGTGACCGTAAATACTggggg                                                                                                                                            <  2:1387303/150‑2 (MQ=255)
         aaCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGc                                                                                                                                            >  2:3030024/1‑150 (MQ=255)
         aaCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGc                                                                                                                                            >  1:4315901/1‑150 (MQ=255)
          aCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACTTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCg                                                                                                                                           >  1:1526198/1‑150 (MQ=255)
             aaaaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATa                                                                                                                                        <  1:4190747/150‑1 (MQ=255)
              aaaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATaa                                                                                                                                       >  2:1828237/1‑150 (MQ=255)
               aaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAac                                                                                                                                      >  1:2796305/1‑150 (MQ=255)
                 ccAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccc                                                                                                                                    >  2:42306/1‑150 (MQ=255)
                   aaGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAacccac                                                                                                                                  <  1:2631032/150‑1 (MQ=255)
                    aGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAacccacc                                                                                                                                 >  1:1005575/1‑150 (MQ=255)
                      tCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCAcccc                                                                                                                               <  2:1296933/150‑1 (MQ=255)
                       ccGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCt                                                                                                                              <  1:3003930/150‑1 (MQ=255)
                          tagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCg                                                                                                                           <  1:2151782/150‑1 (MQ=255)
                             taTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCa                                                                                                                        <  1:1233272/150‑1 (MQ=255)
                              aTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCaa                                                                                                                       >  2:4017290/1‑150 (MQ=255)
                                 aaaCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAAttt                                                                                                                    >  2:168434/1‑150 (MQ=255)
                                   aCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTcc                                                                                                                  >  2:3379539/1‑150 (MQ=255)
                                       gCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCACCGGCAATTTCCCCGc                                                                                                              >  1:4640588/1‑150 (MQ=255)
                                          aTTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAAATTCCCCGcccg                                                                                                           >  2:3853961/1‑147 (MQ=255)
                                           ttCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTc                                                                                                          <  2:2271841/150‑1 (MQ=255)
                                           ttCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTc                                                                                                          >  2:237208/1‑150 (MQ=255)
                                            tCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCa                                                                                                         >  2:1472704/1‑150 (MQ=255)
                                            tCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCCGCAATTTCCCCGcccgca                                                                                                         >  1:2483642/1‑145 (MQ=255)
                                              tGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCAtt                                                                                                       <  2:2447439/150‑1 (MQ=255)
                                               gAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTc                                                                                                      <  2:1617092/150‑1 (MQ=255)
                                               gAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCGCCGGCAATTTCCCCGCACACCTTc                                                                                                      >  1:475004/1‑150 (MQ=255)
                                                  tatTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCtt                                                                                                   >  1:527181/1‑150 (MQ=255)
                                                  tatTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACACATTCCtt                                                                                                   >  1:2916610/1‑150 (MQ=255)
                                                   atTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTc                                                                                                  >  2:1039605/1‑150 (MQ=255)
                                                   atTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTc                                                                                                  >  2:490466/1‑150 (MQ=255)
                                                   atTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTc                                                                                                  <  2:1489565/150‑1 (MQ=255)
                                                    ttGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCa                                                                                                 >  1:109793/1‑150 (MQ=255)
                                                    ttGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCa                                                                                                 <  1:3802670/150‑1 (MQ=255)
                                                    ttGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACCCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCGTCa                                                                                                 >  1:2504504/1‑150 (MQ=255)
                                                         cGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCcca                                                                                            <  2:4640897/150‑1 (MQ=255)
                                                         cGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCcca                                                                                            <  2:675321/150‑1 (MQ=255)
                                                             cGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCaa                                                                                        <  2:4112214/150‑1 (MQ=255)
                                                              ggATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCCATTTCCCCGCACTCATTCCTTCACCCCACCaaa                                                                                       >  1:2252703/1‑150 (MQ=255)
                                                              ggATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCaaa                                                                                       <  2:3405376/150‑1 (MQ=255)
                                                              ggATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACGCCTTCCTTCACCCCACCcaa                                                                                       >  1:369622/1‑150 (MQ=255)
                                                               gATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAAc                                                                                      >  2:1577437/1‑150 (MQ=255)
                                                                aTAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACg                                                                                     >  2:2544773/1‑150 (MQ=255)
                                                                aTAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCACCGGCCATTTACCCGCACACATTCCGTCACCCCACCAAACg                                                                                     >  1:2369213/1‑150 (MQ=255)
                                                                   agagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGctt                                                                                  >  2:1851323/1‑150 (MQ=255)
                                                                    gagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttc                                                                                 <  2:1640634/150‑1 (MQ=255)
                                                                      ggCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttctt                                                                               >  2:4459527/1‑150 (MQ=255)
                                                                      ggCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttctt                                                                               >  2:2055755/1‑150 (MQ=255)
                                                                      ggCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttctt                                                                               >  2:1507668/1‑150 (MQ=255)
                                                                         aaaTAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCACCGGCAATTTCCCCGCACTCATTCCTTCACCCCCCCAAACGCTTCTTCCa                                                                            >  2:712947/1‑150 (MQ=255)
                                                                           ataatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGc                                                                          <  2:1127019/150‑1 (MQ=255)
                                                                             aatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCaa                                                                        >  2:3712577/1‑150 (MQ=255)
                                                                             aatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCaa                                                                        >  2:551827/1‑150 (MQ=255)
                                                                             aatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCCGCAATTTACCCCCACTCATTCCATCACCCCAACAAACGCTTCATCCAGCaa                                                                        >  2:91667/1‑150 (MQ=255)
                                                                             aatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCaa                                                                        >  1:3452087/1‑150 (MQ=255)
                                                                              atcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAAc                                                                       <  1:4431723/150‑1 (MQ=255)
                                                                                     tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTACTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGcg                                                                >  2:796772/1‑149 (MQ=255)
                                                                                        aGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCa                                                             >  1:2464058/1‑150 (MQ=255)
                                                                                          aTAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTACCCGCCCACATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACAGACGCTGCaaa                                                           >  1:3575996/1‑150 (MQ=255)
                                                                                           tAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAAt                                                          <  1:2982036/150‑1 (MQ=255)
                                                                                           tAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAAt                                                          >  2:2020265/1‑150 (MQ=255)
                                                                                             aTTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCACCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCATCTTCCAGCAACTGACGCTGCAAATCg                                                        >  2:3135680/1‑150 (MQ=255)
                                                                                               ttATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGc                                                      >  2:2526333/1‑150 (MQ=255)
                                                                                               ttATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTACTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGc                                                      >  1:416718/1‑150 (MQ=255)
                                                                                               ttATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCACCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGc                                                      >  1:600205/1‑150 (MQ=255)
                                                                                                tATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCACCCGCAAATTCCCCGCCCGCAATCCATCACCCCACCAAACCCATCATCCAGCAACAGACGCAGCAAATCGGcc                                                     >  1:796787/1‑150 (MQ=255)
                                                                                                tATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCAATCATTCCGTCACCCCACCAAACGCTTCTTCCACCCACTGACGCTGCAAATCGGcc                                                     >  1:1686520/1‑150 (MQ=255)
                                                                                                  gtGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTc                                                   <  1:3860260/149‑1 (MQ=255)
                                                                                                     ccTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtc                                                <  1:1631567/150‑1 (MQ=255)
                                                                                                     ccTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtc                                                >  2:1416669/1‑150 (MQ=255)
                                                                                                     ccTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtc                                                <  1:977252/150‑1 (MQ=255)
                                                                                                     ccTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCCATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtc                                                >  1:3506465/1‑150 (MQ=255)
                                                                                                      cTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcg                                               >  2:1001324/1‑150 (MQ=255)
                                                                                                      cTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcg                                               >  2:1970313/1‑150 (MQ=255)
                                                                                                       tGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTACCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACAGACGCTGCAAATCGGCCTCAtcgc                                              >  1:3260827/1‑150 (MQ=255)
                                                                                                        gCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcgct                                             <  1:706341/150‑1 (MQ=255)
                                                                                                         ctctTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCACCGGCAATTTCCCCGCCCTCATTACTTCACCCCACCAAACGCATCTTACAGCCACTGACGCTGCCAATAGGCCACAacgctc                                            >  1:34454/1‑150 (MQ=255)
                                                                                                               tttATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGccccc                                      <  2:2778628/150‑1 (MQ=255)
                                                                                                                ttATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCa                                     <  2:2502563/150‑1 (MQ=255)
                                                                                                                 tataCGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAg                                    <  2:521290/150‑1 (MQ=255)
                                                                                                                   taCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCGTCTTCCAGCAACTGACGCTGCAAATCGGCCTCCTCGCTCGccccccggt                                  >  2:3806652/1‑146 (MQ=255)
                                                                                                                    aCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTc                                 <  1:78214/150‑1 (MQ=255)
                                                                                                                     cGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCa                                <  2:2822865/150‑1 (MQ=255)
                                                                                                                      gTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCAc                               <  2:3194674/150‑1 (MQ=255)
                                                                                                                         aCTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCTGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCa                            <  1:4651687/150‑1 (MQ=255)
                                                                                                                             ctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAg                        <  1:1479320/150‑1 (MQ=255)
                                                                                                                              tcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTACCCGCACTCATTACTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAgc                       >  2:3574973/1‑150 (MQ=255)
                                                                                                                               cGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCCCTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAgcg                      >  1:839078/1‑150 (MQ=255)
                                                                                                                                  aTTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCAt                   <  1:2389062/150‑1 (MQ=255)
                                                                                                                                  aTTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATAGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCAt                   >  1:1239231/1‑150 (MQ=255)
                                                                                                                                    tAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCCGCGCATcc                 >  2:2895126/1‑150 (MQ=255)
                                                                                                                                          tgtgACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCCGTTCACGCATCAGCGCATCCAGTTCa           >  2:1191029/1‑150 (MQ=255)
                                                                                                                                          tgtgACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCa           >  2:4568158/1‑150 (MQ=255)
                                                                                                                                              aCCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACa       <  1:1905826/150‑1 (MQ=255)
                                                                                                                                              aCCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACa       <  2:852814/150‑1 (MQ=255)
                                                                                                                                                 gTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAgcg    >  2:3458752/1‑150 (MQ=255)
                                                                                                                                                   aGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGcc  <  2:109793/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATTTTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGACAGCGCC  >  CP009273/4571085‑4571381

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