Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 321,223 C→A 100% intergenic (‑148/+65) ykgH ← / ← betA putative inner membrane protein/choline dehydrogenase, a flavoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273321,2230CA97.7% 890.8 / 9.7 260intergenic (‑148/+65)ykgH/betAputative inner membrane protein/choline dehydrogenase, a flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (4/2);  new base A (142/112);  total (146/114)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 260 total aligned reads displayed.

ATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATC  >  CP009273/321077‑321369
                                                                                                                                                  |                                                                                                                                                  
 tatatCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTc                                                                                                                                                <  2:616921/150‑1 (MQ=255)
 tatatCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTc                                                                                                                                                >  2:2698866/1‑150 (MQ=255)
  atatCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTcc                                                                                                                                               <  2:3729305/150‑1 (MQ=255)
   tatCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTccc                                                                                                                                              >  1:3317313/1‑150 (MQ=255)
   tatCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTccc                                                                                                                                              >  2:485030/1‑150 (MQ=255)
    atCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTcccc                                                                                                                                             >  2:3898006/1‑150 (MQ=255)
    atCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTcccc                                                                                                                                             >  1:2616379/1‑150 (MQ=255)
        agagAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATTTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACATCAGTGTTTATATAAGGATTAGGTCCACctct                                                                                                                                         >  1:407454/1‑150 (MQ=255)
          agaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctct                                                                                                                                       <  2:3237053/150‑1 (MQ=255)
          agaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctct                                                                                                                                       <  2:3213817/150‑1 (MQ=255)
             aTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGa                                                                                                                                    >  2:1721558/1‑150 (MQ=255)
                aaaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACg                                                                                                                                 >  2:419242/1‑150 (MQ=255)
                aaaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACg                                                                                                                                 >  1:1934972/1‑150 (MQ=255)
                aaaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACg                                                                                                                                 >  1:1282619/1‑150 (MQ=255)
                aaaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTACCCCTCTCTGGAACg                                                                                                                                 >  2:3763106/1‑150 (MQ=255)
                      gTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGTTCCCACTCTCTGGAACGGTAAct                                                                                                                           >  1:259324/1‑150 (MQ=255)
                         gCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACtctc                                                                                                                        >  1:2264254/1‑150 (MQ=255)
                              tAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCt                                                                                                                   <  1:3801655/150‑1 (MQ=255)
                                aGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGa                                                                                                                 >  1:3434256/1‑150 (MQ=255)
                                    tGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCgg                                                                                                             >  2:448275/1‑150 (MQ=255)
                                    tGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCgg                                                                                                             >  2:918637/1‑150 (MQ=255)
                                      tATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGtt                                                                                                           >  1:820908/1‑150 (MQ=255)
                                      tATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGtt                                                                                                           >  1:3136702/1‑150 (MQ=255)
                                      tATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGtt                                                                                                           >  1:925497/1‑150 (MQ=255)
                                       aTTAAGAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTc                                                                                                          <  2:86828/150‑1 (MQ=255)
                                        ttAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTcc                                                                                                         <  2:956712/150‑1 (MQ=255)
                                        ttAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTcc                                                                                                         <  2:2794492/150‑1 (MQ=255)
                                          aaTAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTg                                                                                                       <  2:3246185/150‑1 (MQ=255)
                                          aaTAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTg                                                                                                       <  2:754093/150‑1 (MQ=255)
                                           aTAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGc                                                                                                      >  2:567319/1‑150 (MQ=255)
                                            tAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCg                                                                                                     <  1:888798/150‑1 (MQ=255)
                                             aGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGt                                                                                                    >  1:1428920/1‑150 (MQ=255)
                                               tCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTa                                                                                                  <  1:1794101/150‑1 (MQ=255)
                                                 tGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGt                                                                                                <  2:2581368/150‑1 (MQ=255)
                                                   tATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTc                                                                                              <  2:226607/150‑1 (MQ=255)
                                                      cAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTcaca                                                                                           <  2:1213046/150‑1 (MQ=255)
                                                      cAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTcaca                                                                                           >  1:529836/1‑150 (MQ=255)
                                                      cAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTcaca                                                                                           >  2:1132689/1‑150 (MQ=255)
                                                       aGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACAt                                                                                          >  2:4130795/1‑150 (MQ=255)
                                                       aGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACAt                                                                                          >  2:3533158/1‑150 (MQ=255)
                                                       aGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACAt                                                                                          >  2:3684645/1‑150 (MQ=255)
                                                         gAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCa                                                                                        <  1:3155565/150‑1 (MQ=255)
                                                         gAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCa                                                                                        >  2:389153/1‑150 (MQ=255)
                                                            tGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGa                                                                                     >  1:3905911/1‑150 (MQ=255)
                                                                 gggTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAt                                                                                >  2:1001745/1‑150 (MQ=255)
                                                                 gggTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAt                                                                                >  1:4130908/1‑150 (MQ=255)
                                                                    ttAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAtttt                                                                             <  1:1862114/150‑1 (MQ=255)
                                                                     tAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAttttt                                                                            >  2:2674773/1‑150 (MQ=255)
                                                                     tAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAttttt                                                                            >  2:165751/1‑150 (MQ=255)
                                                                     tAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAttttt                                                                            <  1:559007/150‑1 (MQ=255)
                                                                        atatCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCg                                                                         <  1:2518515/150‑1 (MQ=255)
                                                                        atatCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCg                                                                         <  2:1772457/150‑1 (MQ=255)
                                                                          atCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGct                                                                       <  2:1282619/150‑1 (MQ=255)
                                                                           tCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGctc                                                                      >  1:2932230/1‑150 (MQ=255)
                                                                             aGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCtctc                                                                    <  2:909137/150‑1 (MQ=255)
                                                                             aGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCtctc                                                                    >  2:4040247/1‑150 (MQ=255)
                                                                                aaaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGTACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAcc                                                                 >  1:1508758/1‑150 (MQ=255)
                                                                                aaaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAcc                                                                 >  1:1763869/1‑150 (MQ=255)
                                                                                aaaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAcc                                                                 <  2:2533436/150‑1 (MQ=255)
                                                                                aaaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAcc                                                                 <  2:3279011/150‑1 (MQ=255)
                                                                                   aaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGc                                                              <  1:3535911/150‑1 (MQ=255)
                                                                                   aaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGc                                                              >  2:3767961/1‑150 (MQ=255)
                                                                                   aaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGc                                                              >  1:3890479/1‑150 (MQ=255)
                                                                                          catcatGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCAt                                                       <  2:3615333/150‑1 (MQ=255)
                                                                                           atcatGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCAtt                                                      >  2:3073657/1‑150 (MQ=255)
                                                                                           atcatGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCAtt                                                      <  2:2969070/150‑1 (MQ=255)
                                                                                             catGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTg                                                    >  2:72613/1‑150 (MQ=255)
                                                                                             catGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTg                                                    >  2:3943466/1‑150 (MQ=255)
                                                                                                gAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGTATCCCATTTGCca                                                 <  2:1793570/150‑1 (MQ=255)
                                                                                                gAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCca                                                 <  2:1465986/150‑1 (MQ=255)
                                                                                                     cTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACaaaa                                            <  2:2390852/150‑1 (MQ=255)
                                                                                                        gATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAatat                                         <  1:1711094/150‑1 (MQ=255)
                                                                                                         atatGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATc                                        >  2:522054/1‑150 (MQ=255)
                                                                                                          tatGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATcc                                       >  1:3049574/1‑150 (MQ=255)
                                                                                                              aaGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgcc                                   >  2:2746920/1‑150 (MQ=255)
                                                                                                              aaGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgcc                                   >  2:1150567/1‑150 (MQ=255)
                                                                                                              aaGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgcc                                   >  1:3262822/1‑150 (MQ=255)
                                                                                                              aaGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgcc                                   >  1:675655/1‑150 (MQ=255)
                                                                                                                gCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGt                                 >  1:377020/1‑150 (MQ=255)
                                                                                                                   tGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGct                              <  2:1582767/150‑1 (MQ=255)
                                                                                                                    gagaGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctc                             >  1:1560433/1‑150 (MQ=255)
                                                                                                                       agTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctcctc                          >  1:1337273/1‑150 (MQ=255)
                                                                                                                       agTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctcctc                          <  1:2179368/150‑1 (MQ=255)
                                                                                                                        gTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCg                         >  2:2614961/1‑150 (MQ=255)
                                                                                                                               cAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAgcgc                  >  1:1080304/1‑150 (MQ=255)
                                                                                                                                aGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCt                 >  2:2557568/1‑150 (MQ=255)
                                                                                                                                aGTGTTTATATAAGGATTATGTCCCCCTCTCTTGAACGGTAACTCTCTATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAACATATCACGCCGTGCTCCTCGGCAGCGCt                 <  1:203856/150‑1 (MQ=255)
                                                                                                                                aGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCt                 <  1:3611061/150‑1 (MQ=255)
                                                                                                                                aGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCt                 >  1:1009168/1‑150 (MQ=255)
                                                                                                                                 gtgtTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCCTTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCtt                <  2:3600194/150‑1 (MQ=255)
                                                                                                                                    tttATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCt             >  1:3184867/1‑150 (MQ=255)
                                                                                                                                        tataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCNCGCCGTGCTCCTCGGCAGCGCTTCCTGTcc         <  1:3851385/150‑1 (MQ=255)
                                                                                                                                        tataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTcc         <  1:2349035/150‑1 (MQ=255)
                                                                                                                                         ataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCa        >  2:1080058/1‑150 (MQ=255)
                                                                                                                                         ataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCa        >  2:1849192/1‑150 (MQ=255)
                                                                                                                                          taAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGGTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCAc       >  2:3913832/1‑150 (MQ=255)
                                                                                                                                           aaGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACg      <  1:791191/150‑1 (MQ=255)
                                                                                                                                           aaGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACg      <  2:1990121/150‑1 (MQ=255)
                                                                                                                                             ggATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGaa    <  1:2757578/150‑1 (MQ=255)
                                                                                                                                             ggATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGaa    >  2:208877/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATC  >  CP009273/321077‑321369

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: