Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,182,647 G→T 100% pseudogene (1868/2445 nt) yqiG → pseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,182,6470GT98.6% 672.6 / 0.2 208pseudogene (1868/2445 nt)yqiGpseudogene; fimbrial export usher family;putative membrane; Not classified; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base T (99/106);  total (101/107)
Fisher's exact test for biased strand distribution p-value = 6.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 209 total aligned reads displayed.

GTATGACTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAACAA  >  CP009273/3182501‑3182791
                                                                                                                                                  |                                                                                                                                                
 tATGACTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGAc                                                                                                                                              >  1:3540501/1‑150 (MQ=255)
     aCTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGtt                                                                                                                                          >  2:1950003/1‑150 (MQ=255)
      ctctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGtta                                                                                                                                         <  2:459973/150‑1 (MQ=255)
      ctctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGtta                                                                                                                                         <  1:2335389/150‑1 (MQ=255)
       tctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGttat                                                                                                                                        <  1:3768359/150‑1 (MQ=255)
       tctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGttat                                                                                                                                        >  2:413803/1‑150 (MQ=255)
        ctTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGttatt                                                                                                                                       <  1:1282643/150‑1 (MQ=255)
           ccGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTattatt                                                                                                                                    >  1:507864/1‑150 (MQ=255)
             gATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATtct                                                                                                                                  <  2:3652514/150‑1 (MQ=255)
              aTCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATtctc                                                                                                                                 <  1:4689985/150‑1 (MQ=255)
              aTCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATtctc                                                                                                                                 <  1:693244/150‑1 (MQ=255)
               tCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTctct                                                                                                                                <  2:3386448/150‑1 (MQ=255)
               tCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTctct                                                                                                                                <  2:2890833/150‑1 (MQ=255)
               tCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTctct                                                                                                                                >  1:958368/1‑150 (MQ=255)
                cGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTg                                                                                                                               <  1:3791619/150‑1 (MQ=255)
                    ataatTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATg                                                                                                                           >  1:176257/1‑150 (MQ=255)
                       atTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCa                                                                                                                        <  1:2893259/150‑1 (MQ=255)
                       atTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCa                                                                                                                        <  1:4222295/150‑1 (MQ=255)
                       atTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCa                                                                                                                        >  2:3090769/1‑150 (MQ=255)
                          ccTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGtt                                                                                                                     <  1:2935612/150‑1 (MQ=255)
                           cTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTg                                                                                                                    >  2:1152543/1‑150 (MQ=255)
                            tGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTgg                                                                                                                   <  2:4525642/150‑1 (MQ=255)
                               aGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAAt                                                                                                                <  1:1829923/150‑1 (MQ=255)
                                  gTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGt                                                                                                             >  1:3356397/1‑150 (MQ=255)
                                  gTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGt                                                                                                             >  1:1788418/1‑150 (MQ=255)
                                  gTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCAAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGt                                                                                                             <  1:4511632/150‑1 (MQ=255)
                                   tttCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGtt                                                                                                            <  2:3585326/150‑1 (MQ=255)
                                       tGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcatt                                                                                                        >  2:3307484/1‑150 (MQ=255)
                                       tGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcatt                                                                                                        >  2:494351/1‑150 (MQ=255)
                                        gCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCCACCTCTCCGGTACATCCCAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcattc                                                                                                       <  2:2714121/150‑1 (MQ=255)
                                         cTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGACCGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcattca                                                                                                      >  2:3477386/1‑150 (MQ=255)
                                         cTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcattca                                                                                                      >  2:1239400/1‑150 (MQ=255)
                                         cTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGttcattca                                                                                                      >  1:1227797/1‑150 (MQ=255)
                                           tCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACt                                                                                                    >  2:1915038/1‑150 (MQ=255)
                                           tCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACt                                                                                                    <  1:1593602/150‑1 (MQ=255)
                                            cAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTg                                                                                                   >  1:4010160/1‑150 (MQ=255)
                                             aGGCGACAATGATTCATGAAAACATATGAAAGCGTCACTAACCGCCAGTTAACATCATAATACCTAGAACTGTCGAATCTACCTCTACGGTACATCACAGCTAGACAGTTATTATTCTCTGAATTCCAGTTGGAATGGTTCATTCACTGc                                                                                                  <  2:5088676/150‑1 (MQ=255)
                                              ggCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCg                                                                                                 <  1:4964161/150‑1 (MQ=255)
                                               gCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGa                                                                                                <  2:811779/150‑1 (MQ=255)
                                                   cAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCg                                                                                            <  1:2888643/150‑1 (MQ=255)
                                                        atgaATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACg                                                                                       <  2:1488708/150‑1 (MQ=255)
                                                        atgaATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACg                                                                                       <  2:4097798/150‑1 (MQ=255)
                                                         tgaATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACgg                                                                                      >  1:4395887/1‑150 (MQ=255)
                                                           aaTTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTg                                                                                    >  2:1717361/1‑150 (MQ=255)
                                                           aaTTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTg                                                                                    >  1:3486140/1‑150 (MQ=255)
                                                             ttCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGcc                                                                                  <  1:3248604/150‑1 (MQ=255)
                                                               cAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGc                                                                                <  1:4123744/150‑1 (MQ=255)
                                                                 aaGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCtt                                                                              >  1:3187410/1‑150 (MQ=255)
                                                                 aaGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCtt                                                                              <  2:507779/150‑1 (MQ=255)
                                                                    atatGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTcc                                                                           >  2:3360944/1‑150 (MQ=255)
                                                                      atGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCAc                                                                         <  1:1357940/150‑1 (MQ=255)
                                                                       tGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACg                                                                        <  2:1550417/150‑1 (MQ=255)
                                                                          aaGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACt                                                                     <  1:1291912/150‑1 (MQ=255)
                                                                          aaGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACt                                                                     >  2:2972871/1‑150 (MQ=255)
                                                                             cGTCACTACGCGCCAGTTATCAGCATAATACCGTGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACtata                                                                  <  1:1065177/150‑1 (MQ=255)
                                                                             cGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACtata                                                                  >  2:2143282/1‑150 (MQ=255)
                                                                               tCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGc                                                                <  2:4864265/150‑1 (MQ=255)
                                                                               tCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGc                                                                >  2:899253/1‑150 (MQ=255)
                                                                               tCAATACGCGCCTGTTATAATCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGc                                                                <  1:884524/150‑1 (MQ=255)
                                                                                cACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCg                                                               >  2:3632205/1‑150 (MQ=255)
                                                                                cACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCg                                                               >  2:2459114/1‑150 (MQ=255)
                                                                                  cTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCgga                                                             <  2:746148/150‑2 (MQ=255)
                                                                                  cTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCgga                                                             <  1:916937/150‑2 (MQ=255)
                                                                                  cTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCgga                                                             >  2:4582451/1‑149 (MQ=255)
                                                                                    aCGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAg                                                           >  2:3983637/1‑150 (MQ=255)
                                                                                     cgcgCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGt                                                          >  2:1987646/1‑150 (MQ=255)
                                                                                        gcCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCt                                                       >  2:4294576/1‑150 (MQ=255)
                                                                                         ccAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTg                                                      <  1:1680401/150‑1 (MQ=255)
                                                                                         ccAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTg                                                      <  1:4926248/150‑1 (MQ=255)
                                                                                               aTCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTcgc                                                <  1:4223726/150‑1 (MQ=255)
                                                                                                 cAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGt                                              >  1:3437347/1‑150 (MQ=255)
                                                                                                    cATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATTGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGtttt                                           <  2:1214141/150‑1 (MQ=255)
                                                                                                       aataCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATg                                        >  2:474600/1‑150 (MQ=255)
                                                                                                          aCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATc                                     <  2:3905636/150‑1 (MQ=255)
                                                                                                          aCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATc                                     <  2:2307781/150‑1 (MQ=255)
                                                                                                          aCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATc                                     <  2:4395773/150‑1 (MQ=255)
                                                                                                          aCCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATc                                     <  1:2835409/150‑1 (MQ=255)
                                                                                                           ccGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCg                                    <  2:1621263/150‑1 (MQ=255)
                                                                                                             gagaACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGAc                                  >  1:1566117/1‑150 (MQ=255)
                                                                                                             gagaACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGAc                                  >  1:2448297/1‑150 (MQ=255)
                                                                                                               gaACGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGc                                >  2:1066467/1‑150 (MQ=255)
                                                                                                                aaCGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCa                               >  2:3846894/1‑150 (MQ=255)
                                                                                                                 aCGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAg                              <  1:473368/150‑1 (MQ=255)
                                                                                                                  cGGTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGa                             <  1:2908342/150‑1 (MQ=255)
                                                                                                                   ggTCGACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAGAc                            <  1:2161407/150‑1 (MQ=255)
                                                                                                                       gACTCTACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCa                        <  1:1717361/150‑1 (MQ=255)
                                                                                                                         ctctACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACt                      <  2:139507/150‑1 (MQ=255)
                                                                                                                           ctACCTCTCCGTTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGa                    <  2:3273812/150‑1 (MQ=255)
                                                                                                                            tACCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGaa                   <  1:3876690/150‑1 (MQ=255)
                                                                                                                             aCCTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAg                  <  2:1330042/150‑1 (MQ=255)
                                                                                                                              ccTCTCCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGa                 >  1:4351746/1‑150 (MQ=255)
                                                                                                                              ccTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGa                 >  1:1419405/1‑150 (MQ=255)
                                                                                                                                 ctcCGGTACATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGatat              <  2:4450539/150‑1 (MQ=255)
                                                                                                                                        aCATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTg       >  1:4417925/1‑150 (MQ=255)
                                                                                                                                        aCATCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTg       >  1:392216/1‑150 (MQ=255)
                                                                                                                                          aTCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGaa     >  1:1871604/1‑150 (MQ=255)
                                                                                                                                          aTCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGaa     <  1:5139602/150‑1 (MQ=255)
                                                                                                                                           tCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAAc    >  2:1868001/1‑150 (MQ=255)
                                                                                                                                           tCACAGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGAAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAAc    <  1:2441997/150‑1 (MQ=255)
                                                                                                                                            cacaGCTAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAACa   >  2:320151/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTATGACTCTTCCGATCGAAATAATTCCTGGAGCGTTTCTGCTTCAGGCGACAATGATGAATTCAAAGATATGAAAGCGTCACTACGCGCCAGTTATCAGCATAATACCGAGAACGGTCGACTCTACCTCTCCGGTACATCACAGCGAGACAGTTATTATTCTCTGAATGCCAGTTGGAATGGTTCATTCACTGCGACTCGCCACGGTGCCGCTTTCCACGACTATAGCGGTAGTGCTGACTCGCGTTTTATGATCGACGCAGACGGCACTGAAGATATTCCGTTGAACAA  >  CP009273/3182501‑3182791

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: