Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 747,068 A→T 100% I186I (ATT→ATA ybgQ ← putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273747,0680AT97.5% 794.2 / 6.4 243I186I (ATT→ATAybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base A (4/2);  new base T (132/105);  total (136/107)
Fisher's exact test for biased strand distribution p-value = 6.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 243 total aligned reads displayed.

TCGCCGAGGGTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTG  >  CP009273/746925‑747211
                                                                                                                                               |                                                                                                                                               
tCGCCGAGGGTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTc                                                                                                                                             <  2:2853427/148‑1 (MQ=255)
   ccGAGGGTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCcgg                                                                                                                                          >  1:1028184/1‑148 (MQ=255)
     gAGGGTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCcggcg                                                                                                                                        >  2:3364744/1‑148 (MQ=255)
     gAGGGTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCcggcg                                                                                                                                        >  2:2920244/1‑148 (MQ=255)
      aGGGTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCcggcgg                                                                                                                                       >  1:784813/1‑148 (MQ=255)
       gggTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCTCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCCGCGTTTATTCCGGCGGt                                                                                                                                      >  2:904627/1‑148 (MQ=255)
       gggTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGt                                                                                                                                      <  1:1034901/148‑1 (MQ=255)
       gggTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGt                                                                                                                                      >  2:904581/1‑148 (MQ=255)
        ggTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTa                                                                                                                                     <  2:3422020/148‑1 (MQ=255)
        ggTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTa                                                                                                                                     <  2:2255062/148‑1 (MQ=255)
        ggTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTa                                                                                                                                     <  2:2255037/148‑1 (MQ=255)
           tAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCg                                                                                                                                  >  1:1918392/1‑148 (MQ=255)
             aCTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTa                                                                                                                                <  2:813791/148‑1 (MQ=255)
             aCTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTa                                                                                                                                >  1:1477929/1‑148 (MQ=255)
              cTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAg                                                                                                                               >  2:3047511/1‑148 (MQ=255)
              cTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAg                                                                                                                               <  1:2219686/148‑1 (MQ=255)
               ttAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAgg                                                                                                                              >  1:2849477/1‑148 (MQ=255)
                   agCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCAtt                                                                                                                          >  2:2110301/1‑148 (MQ=255)
                   agCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCAtt                                                                                                                          >  2:3422834/1‑148 (MQ=255)
                      cTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAg                                                                                                                       <  2:1508139/148‑1 (MQ=255)
                       tAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAgg                                                                                                                      <  1:2602913/148‑1 (MQ=255)
                       tAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAgg                                                                                                                      >  2:2299618/1‑148 (MQ=255)
                        aaTTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGt                                                                                                                     >  1:961588/1‑148 (MQ=255)
                           tGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAg                                                                                                                  >  1:741648/1‑148 (MQ=255)
                           tGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAg                                                                                                                  >  1:1242346/1‑148 (MQ=255)
                             cGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGtg                                                                                                                <  2:1690464/148‑1 (MQ=255)
                              ggTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGtgc                                                                                                               <  1:2200582/148‑1 (MQ=255)
                                 aaTGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGtgctgc                                                                                                            >  1:2556410/1‑148 (MQ=255)
                                    ggACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGtgctgctgc                                                                                                         <  2:2004517/148‑1 (MQ=255)
                                         aaaaaGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCg                                                                                                    >  2:1354553/1‑148 (MQ=255)
                                          aaaaGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGt                                                                                                   >  1:1755249/1‑148 (MQ=255)
                                          aaaaGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGt                                                                                                   >  1:3089720/1‑148 (MQ=255)
                                           aaaGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTc                                                                                                  <  1:1883416/148‑1 (MQ=255)
                                            aaGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTcc                                                                                                 >  2:3016043/1‑148 (MQ=255)
                                             aGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCt                                                                                                >  1:3211818/1‑148 (MQ=255)
                                              gATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGCCATTCAGGTTAGTGCTGCTGCTGCCGTGCtg                                                                                               >  1:3001710/1‑148 (MQ=255)
                                                  ggTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGgggg                                                                                           >  1:2506505/1‑148 (MQ=255)
                                                   gTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGc                                                                                          <  1:2199156/148‑1 (MQ=255)
                                                    tGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCg                                                                                         <  2:2049259/148‑1 (MQ=255)
                                                     gcgcGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCgg                                                                                        >  2:1924736/1‑148 (MQ=255)
                                                      cgcgATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGt                                                                                       <  2:631567/148‑1 (MQ=255)
                                                        cgATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAg                                                                                     >  2:1985083/1‑148 (MQ=255)
                                                        cgATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAg                                                                                     >  1:826428/1‑148 (MQ=255)
                                                        cgATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAg                                                                                     <  1:1496098/148‑1 (MQ=255)
                                                         gATAGTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAgc                                                                                    <  2:3289648/148‑1 (MQ=255)
                                                           tatTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAgctg                                                                                  >  1:2776267/1‑148 (MQ=255)
                                                           tatTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAgctg                                                                                  >  1:1695847/1‑148 (MQ=255)
                                                           tatTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAgctg                                                                                  >  1:433276/1‑148 (MQ=255)
                                                           tatTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAgctg                                                                                  >  1:2422726/1‑148 (MQ=255)
                                                             ttCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAgctgct                                                                                >  2:116427/1‑148 (MQ=255)
                                                              tCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGctgctg                                                                               >  2:3344880/1‑148 (MQ=255)
                                                                cGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGc                                                                             >  2:1517400/1‑148 (MQ=255)
                                                                    tGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTACGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAAc                                                                         >  1:655342/1‑148 (MQ=255)
                                                                     gACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACa                                                                        <  2:2389268/148‑1 (MQ=255)
                                                                     gACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACa                                                                        <  1:3213001/148‑1 (MQ=255)
                                                                     gACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACa                                                                        >  2:1864217/1‑148 (MQ=255)
                                                                      aCTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACaa                                                                       >  2:603485/1‑148 (MQ=255)
                                                                      aCTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACaa                                                                       >  1:2260808/1‑148 (MQ=255)
                                                                        tGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGt                                                                     >  2:2678689/1‑148 (MQ=255)
                                                                            cATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTa                                                                 <  1:589484/148‑1 (MQ=255)
                                                                                  tATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCAt                                                           <  2:2152819/148‑1 (MQ=255)
                                                                                   aTCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATc                                                          <  1:2395928/148‑1 (MQ=255)
                                                                                   aTCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATc                                                          <  2:1042538/148‑1 (MQ=255)
                                                                                   aTCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCAGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATc                                                          <  2:2454302/148‑1 (MQ=255)
                                                                                         gCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGaca                                                    <  2:2638/148‑1 (MQ=255)
                                                                                           tATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACAcc                                                  >  2:1572456/1‑148 (MQ=255)
                                                                                            aTCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCg                                                 >  2:384504/1‑148 (MQ=255)
                                                                                              cGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGc                                               <  2:366197/148‑1 (MQ=255)
                                                                                               ggTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCa                                              <  2:1133479/148‑1 (MQ=255)
                                                                                                 tCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAac                                            >  1:1334211/1‑148 (MQ=255)
                                                                                                    tattaAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACAcc                                         >  1:3384290/1‑148 (MQ=255)
                                                                                                      ttaAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACctt                                       >  2:2341262/1‑148 (MQ=255)
                                                                                                         aGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTCATCCATCAGGACACCGGCAACACcttctt                                    >  2:323758/1‑148 (MQ=255)
                                                                                                         aGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACcttctt                                    >  1:1989015/1‑148 (MQ=255)
                                                                                                            tGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTcc                                 >  1:2640245/1‑148 (MQ=255)
                                                                                                             ggTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCa                                >  1:2901387/1‑148 (MQ=255)
                                                                                                             ggTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCa                                <  1:367770/148‑1 (MQ=255)
                                                                                                             ggTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTACGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCa                                >  1:2347566/1‑148 (MQ=255)
                                                                                                              gTAATCACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCAt                               <  2:1931892/148‑1 (MQ=255)
                                                                                                                   cACTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCNTCCTGTGGGCGGGAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTaga                          >  1:1784558/1‑148 (MQ=255)
                                                                                                                    aCTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTagag                         <  2:2938139/148‑1 (MQ=255)
                                                                                                                     cTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgg                        <  2:3171356/148‑1 (MQ=255)
                                                                                                                     cTGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAgg                        <  1:1510354/148‑1 (MQ=255)
                                                                                                                      tGCGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAggg                       >  2:3117439/1‑148 (MQ=255)
                                                                                                                        cGTAAGCGCCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAggggg                     <  1:1145025/148‑1 (MQ=255)
                                                                                                                              cgcCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACCAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGccc               >  2:2424613/1‑148 (MQ=255)
                                                                                                                              cgcCATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGccc               >  1:2570262/1‑148 (MQ=255)
                                                                                                                                ccATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAg             <  1:1389463/148‑1 (MQ=255)
                                                                                                                                ccATGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAg             <  1:1388842/148‑1 (MQ=255)
                                                                                                                                  aTGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGCGCCCAGtt           >  1:2447104/1‑148 (MQ=255)
                                                                                                                                  aTGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGGTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGtt           >  1:740326/1‑148 (MQ=255)
                                                                                                                                   tGCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGttt          >  1:2848789/1‑148 (MQ=255)
                                                                                                                                    gCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGtttt         >  1:2043538/1‑148 (MQ=255)
                                                                                                                                    gCCCCGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGtttt         <  2:2422726/148‑1 (MQ=255)
                                                                                                                                        cGGCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGa     <  2:1274457/148‑1 (MQ=255)
                                                                                                                                         ggCGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAg    <  2:1314054/148‑1 (MQ=255)
                                                                                                                                           cGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTg  <  2:1124888/148‑1 (MQ=255)
                                                                                                                                           cGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTg  >  2:1457876/1‑148 (MQ=255)
                                                                                                                                           cGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTg  >  2:2700197/1‑148 (MQ=255)
                                                                                                                                           cGTTTATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTg  <  2:305893/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TCGCCGAGGGTTAACTTAGAGCCTAATTGCGGTAATGGACGAAAAAGATAGGTGCGCGATATTCCGCCTGACTGGTCATGGTTATCTTCGCTATCGGTCTTATTAAGCTGGTAATCACTGCGTAAGCGCCATGCCCCGGCGTTAATTCCGGCGGTACCGTAGGCATTCAGGTTAGTGCTGCTGCTGCCGTCCTGTGGGCGGTAGCTGCTGGCAAACAAGTTGTAATCCATCAGGACACCGGCAACACCTTCTTTCCATGTAGAGGGGGGAGCCCAGTTTTCTGAGTG  >  CP009273/746925‑747211

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