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,2330GA95.2% 663.2 / 18.9 206S16F (TCT→TTT) hsdSspecificity determinant for hsdM and hsdR
Reads supporting (aligned to +/- strand):  ref base G (4/6);  new base A (109/87);  total (113/93)
Fisher's exact test for biased strand distribution p-value = 3.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 206 total aligned reads displayed.

TTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGAC  >  CP009273/4571088‑4571375
                                                                                                                                                 |                                                                                                                                              
 ttAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATa                                                                                                                                             <  1:1889681/148‑1 (MQ=255)
   aGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACt                                                                                                                                           <  2:1631620/148‑1 (MQ=255)
       aCAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTgggg                                                                                                                                       >  2:577079/1‑148 (MQ=255)
         aaaaaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCg                                                                                                                                     <  1:1417261/148‑1 (MQ=255)
         aaaaaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCg                                                                                                                                     >  1:2097242/1‑148 (MQ=255)
           aaaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGAt                                                                                                                                   <  2:3240647/148‑1 (MQ=255)
            aaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATa                                                                                                                                  >  1:585686/1‑148 (MQ=255)
            aaCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATa                                                                                                                                  <  2:1954956/148‑1 (MQ=255)
             aCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATaa                                                                                                                                 >  1:2711380/1‑148 (MQ=255)
             aCCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATaa                                                                                                                                 >  1:1938937/1‑148 (MQ=255)
              ccAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAac                                                                                                                                >  2:1103931/1‑148 (MQ=255)
                aaGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccc                                                                                                                              <  2:1322552/148‑1 (MQ=255)
                 aGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAaccca                                                                                                                             >  2:2730295/1‑148 (MQ=255)
                      gtagtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCt                                                                                                                        <  2:2561080/148‑1 (MQ=255)
                        agtaTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTcc                                                                                                                      <  1:470419/148‑1 (MQ=255)
                          taTCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCgg                                                                                                                    <  1:305590/148‑1 (MQ=255)
                            tCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCa                                                                                                                  <  2:939833/148‑1 (MQ=255)
                            tCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCa                                                                                                                  <  1:2188703/148‑1 (MQ=255)
                             cAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCaa                                                                                                                 <  2:1855987/148‑1 (MQ=255)
                              aaaCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAAt                                                                                                                <  1:476065/148‑1 (MQ=255)
                               aaCTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAAtt                                                                                                               >  2:1267165/1‑148 (MQ=255)
                                  ttGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTcc                                                                                                            <  1:3357393/148‑1 (MQ=255)
                                   tGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTccc                                                                                                           >  2:1936843/1‑148 (MQ=255)
                                            gAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCAt                                                                                                  <  2:1462911/148‑1 (MQ=255)
                                              ctatTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTc                                                                                                >  2:2305350/2‑148 (MQ=255)
                                                atTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCt                                                                                              >  1:1437815/1‑148 (MQ=255)
                                                  tGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTc                                                                                            <  2:1975206/148‑1 (MQ=255)
                                                  tGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTc                                                                                            >  2:2368103/1‑148 (MQ=255)
                                                   gTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCa                                                                                           <  2:3386938/148‑1 (MQ=255)
                                                   gTTCGCACGGATAAGAGGCAAATAATCATCTGTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCa                                                                                           <  1:2453983/148‑1 (MQ=255)
                                                    ttCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCACCGGCAATTTCCCCGCACTCATTCCTTCAc                                                                                          >  1:3179073/1‑148 (MQ=255)
                                                     tCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCAcc                                                                                         <  2:793067/148‑1 (MQ=255)
                                                           ggATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCccacca                                                                                   >  1:2881761/1‑148 (MQ=255)
                                                            gATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCaa                                                                                  >  2:832566/1‑148 (MQ=255)
                                                            gATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCaa                                                                                  <  2:1612350/148‑1 (MQ=255)
                                                             aTAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCaaa                                                                                 <  1:104603/148‑1 (MQ=255)
                                                              tAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCCCGCCTTCCTTCACCCCACCCAAc                                                                                >  1:3154264/1‑148 (MQ=255)
                                                                agagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGc                                                                              >  2:215727/1‑148 (MQ=255)
                                                                agagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGc                                                                              >  1:1062864/1‑148 (MQ=255)
                                                                agagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGc                                                                              <  1:3447538/148‑1 (MQ=255)
                                                                agagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGc                                                                              <  1:3447001/148‑1 (MQ=255)
                                                                agagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGc                                                                              >  1:1159348/1‑148 (MQ=255)
                                                                agagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGc                                                                              <  1:2769889/148‑1 (MQ=255)
                                                                agagGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGc                                                                              >  1:1818836/1‑148 (MQ=255)
                                                                  agGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGctt                                                                            >  1:3437859/1‑148 (MQ=255)
                                                                    gCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGcttct                                                                          >  1:2323955/1‑148 (MQ=255)
                                                                      aaaTAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCttcttc                                                                        <  1:3526477/148‑1 (MQ=255)
                                                                       aataatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTcc                                                                       <  2:671105/148‑1 (MQ=255)
                                                                        ataatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCa                                                                      >  1:911427/1‑148 (MQ=255)
                                                                          aatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGTCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGc                                                                    >  2:2569911/1‑148 (MQ=255)
                                                                          aatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGc                                                                    >  2:33574/1‑148 (MQ=255)
                                                                          aatCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGc                                                                    <  2:2785617/148‑1 (MQ=255)
                                                                           atcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTACCCGCACACATTCCATCACCCCACCCAAAGCTTCTTACAGCa                                                                   >  2:1797854/1‑148 (MQ=255)
                                                                            tcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCaa                                                                  <  1:997772/148‑1 (MQ=255)
                                                                            tcatcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCaa                                                                  <  2:570379/148‑1 (MQ=255)
                                                                             catcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAAc                                                                 <  2:247928/148‑1 (MQ=255)
                                                                              atcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACt                                                                >  2:3168560/1‑148 (MQ=255)
                                                                              atcTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTAATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACt                                                                >  2:816138/1‑148 (MQ=255)
                                                                                cTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGa                                                              >  1:2827533/1‑148 (MQ=255)
                                                                                  tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACg                                                            <  1:1561797/148‑1 (MQ=255)
                                                                                  tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACg                                                            >  1:2529319/1‑148 (MQ=255)
                                                                                  tttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACg                                                            <  1:1078325/148‑1 (MQ=255)
                                                                                   ttAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCAACCCACCCGCAAATTACCCCCACTCAATCCAAAACCCCACCAAACCCATCTACCACCAACTGAAGc                                                           >  1:2376989/1‑148 (MQ=255)
                                                                                     aGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTg                                                         <  1:1731984/148‑1 (MQ=255)
                                                                                      gATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCCGCAATTTCCCCGCACTCATTCCGTCAACCCACCAAACGCTTCTTCCACAAACTGACGCTGc                                                        >  1:3392168/1‑148 (MQ=255)
                                                                                       aTAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCa                                                       >  2:528659/1‑148 (MQ=255)
                                                                                        tAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCaa                                                      >  2:1902050/1‑148 (MQ=255)
                                                                                        tAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCaa                                                      <  2:1951978/148‑1 (MQ=255)
                                                                                           tttATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATc                                                   >  2:130949/1‑148 (MQ=255)
                                                                                             tATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCgg                                                 >  2:3009883/1‑148 (MQ=255)
                                                                                               ttGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGcc                                               <  1:264150/148‑1 (MQ=255)
                                                                                                tGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCt                                              >  2:730069/1‑148 (MQ=255)
                                                                                                tGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTCGGGCGATAACCCACCCCACAGGCAATTTCCCCGCACACATTACTTCAACCCACCAAACCCTACATCCACCAACAGACTCAGCCAATAGGCCt                                              >  1:3328194/1‑148 (MQ=255)
                                                                                                 gCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTc                                             <  1:751254/148‑1 (MQ=255)
                                                                                                  ccTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCa                                            <  2:1356817/148‑1 (MQ=255)
                                                                                                  ccTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCa                                            <  2:62765/148‑1 (MQ=255)
                                                                                                     gCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcg                                         >  1:509218/1‑148 (MQ=255)
                                                                                                      ctctTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcgc                                        >  2:1428578/1‑148 (MQ=255)
                                                                                                      ctctTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCAtcgc                                        >  2:1648953/1‑148 (MQ=255)
                                                                                                        ctTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGACCTCAtcgctc                                      >  1:2738408/1‑148 (MQ=255)
                                                                                                          tttttATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTACCCGCACTCATTCCTTCACCCCACCAAACCCTTCTTCCAGCAACAGACGCTGCAAATCGGCCTCAtcgctcgc                                    >  1:2841491/1‑148 (MQ=255)
                                                                                                           ttttATACGTTACTCCTCGGATTAGAGTTGTGACCTTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGcc                                   <  1:1832403/148‑1 (MQ=255)
                                                                                                             ttATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGcccc                                 <  1:173766/148‑1 (MQ=255)
                                                                                                             ttATACGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGcccc                                 <  1:1265353/148‑1 (MQ=255)
                                                                                                                taCGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCCAATCGGCCTCATCGCTCGCCCCCAg                              >  2:2727073/1‑148 (MQ=255)
                                                                                                                  cGTTACTCCTCGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGtt                            <  2:1723740/148‑1 (MQ=255)
                                                                                                                        tcctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCa                      >  2:2700230/1‑148 (MQ=255)
                                                                                                                        tcctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGCCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCa                      <  2:3322157/148‑1 (MQ=255)
                                                                                                                         cctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCAt                     <  1:2405904/148‑1 (MQ=255)
                                                                                                                          ctcGGATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACAGACGCTGCAAATCGGCCACATCGCTCGCCCCCAGTTCACGCATc                    >  2:2385391/1‑148 (MQ=255)
                                                                                                                              gATTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAgcg                >  1:3206114/1‑148 (MQ=255)
                                                                                                                              gATTAGAGTTGTGACCGTAAAGACTGGCGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAgcg                <  2:3416820/148‑1 (MQ=255)
                                                                                                                               aTTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAgcgc               >  1:1639462/1‑148 (MQ=255)
                                                                                                                               aTTAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTACCCGCACGCAATCCTTCAACCCACCAAAACCTTCTTCCAGCAACTGACGCTGCAAAACGGCCTCATCGCTCGCCCCCAGTTCACGCATCAgcgc               >  1:709033/1‑148 (MQ=255)
                                                                                                                                ttAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCa              >  1:1282324/1‑148 (MQ=255)
                                                                                                                                 tAGAGTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACAGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCAACAGCGCAt             >  2:583850/1‑148 (MQ=255)
                                                                                                                                  agagTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATc            <  1:566356/148‑1 (MQ=255)
                                                                                                                                  agagTTGTGACCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATc            <  2:471568/148‑1 (MQ=255)
                                                                                                                                           aCCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGa   <  1:1398278/148‑1 (MQ=255)
                                                                                                                                           aCCGTAAATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGa   <  1:3072027/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TTTAGGAACAAAAACCAAGTCCGTAGTATCAAACTTGCCATTCTGAATATTGTTCGCACGGATAAGAGGCAAATAATCATCTTTTAGATAATTTATTGCCTGCTCTTTTTTATACGTTACTCCTCGGATTAGAGTTGTGACCGTAGATACTGGGGCGATAACCCACCCCTCCGGCAATTTCCCCGCACTCATTCCTTCACCCCACCAAACGCTTCTTCCAGCAACTGACGCTGCAAATCGGCCTCATCGCTCGCCCCCAGTTCACGCATCAGCGCATCCAGTTCAGAC  >  CP009273/4571088‑4571375

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