Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 321,223 C→A 94.1% 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,2230CA94.1% 689.1 / 28.5 212intergenic (‑148/+65)ykgH/betAputative inner membrane protein/choline dehydrogenase, a flavoprotein
Reads supporting (aligned to +/- strand):  ref base C (10/3);  new base A (107/92);  total (117/95)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 213 total aligned reads displayed.

AATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATCATA  >  CP009273/321076‑321372
                                                                                                                                                   |                                                                                                                                                     
aataTATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAgg                                                                                                                                                     <  1:2404825/150‑1 (MQ=255)
 atatatCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGt                                                                                                                                                    >  2:106151/1‑150 (MQ=255)
  tatatCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTc                                                                                                                                                   >  1:1522711/1‑150 (MQ=255)
     atCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTcccc                                                                                                                                                >  1:1961163/1‑150 (MQ=255)
       ccAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCct                                                                                                                                              <  2:1732782/150‑1 (MQ=255)
       ccAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCct                                                                                                                                              >  2:3210009/1‑150 (MQ=255)
          gagaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctc                                                                                                                                           <  1:853789/150‑1 (MQ=255)
           agaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctct                                                                                                                                          >  1:1147935/1‑150 (MQ=255)
           agaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctct                                                                                                                                          >  2:3166540/1‑150 (MQ=255)
           agaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCctctct                                                                                                                                          <  1:1734942/150‑1 (MQ=255)
            gaATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTg                                                                                                                                         <  1:1573868/150‑1 (MQ=255)
             aaTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTgg                                                                                                                                        >  1:188272/1‑150 (MQ=255)
             aaTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTgg                                                                                                                                        <  1:982826/150‑1 (MQ=255)
             aaTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTgg                                                                                                                                        >  2:738221/1‑150 (MQ=255)
              aTAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGa                                                                                                                                       <  2:3600099/150‑1 (MQ=255)
                 aaaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACg                                                                                                                                    <  2:2285029/150‑1 (MQ=255)
                 aaaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACg                                                                                                                                    >  2:3752675/1‑150 (MQ=255)
                 aaaTCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACg                                                                                                                                    >  2:4331064/1‑150 (MQ=255)
                       gTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAAct                                                                                                                              <  1:900843/150‑1 (MQ=255)
                       gTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAAct                                                                                                                              <  1:4157777/150‑1 (MQ=255)
                       gTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAAct                                                                                                                              >  1:2587547/1‑150 (MQ=255)
                       gTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAAct                                                                                                                              <  1:900857/150‑1 (MQ=255)
                       gTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAAct                                                                                                                              <  1:1968621/150‑1 (MQ=255)
                        tCGCAGATAATTTTGTATTTATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAATGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAActc                                                                                                                             <  2:1029267/150‑1 (MQ=255)
                         cGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAActct                                                                                                                            <  2:3424389/150‑1 (MQ=255)
                          gCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACtctc                                                                                                                           <  2:2941499/150‑1 (MQ=255)
                            aGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCaa                                                                                                                         <  1:1075856/150‑1 (MQ=255)
                            aGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCaa                                                                                                                         >  2:307110/1‑150 (MQ=255)
                             gATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAAt                                                                                                                        >  2:2244907/1‑150 (MQ=255)
                                 aGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGa                                                                                                                    >  2:4259048/1‑150 (MQ=255)
                                    ttGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCg                                                                                                                 <  2:3240241/150‑1 (MQ=255)
                                      gTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGt                                                                                                               <  1:1666169/150‑1 (MQ=255)
                                         ttAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTcc                                                                                                            >  1:2029856/1‑150 (MQ=255)
                                          tAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCt                                                                                                           <  1:2928744/150‑1 (MQ=255)
                                           aaTAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTg                                                                                                          >  1:2683231/1‑150 (MQ=255)
                                           aaTAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTg                                                                                                          >  2:949736/1‑150 (MQ=255)
                                               gTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGtt                                                                                                      >  1:2783098/1‑150 (MQ=255)
                                                 cTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAg                                                                                                    <  2:1538288/150‑1 (MQ=255)
                                                 cTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAg                                                                                                    >  1:1183560/1‑150 (MQ=255)
                                                     aTCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTca                                                                                                >  1:945306/1‑150 (MQ=255)
                                                      tCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTcac                                                                                               <  2:1632818/150‑1 (MQ=255)
                                                      tCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTcac                                                                                               >  1:400052/1‑150 (MQ=255)
                                                       cAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTcaca                                                                                              >  1:1904233/1‑150 (MQ=255)
                                                        aGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACAt                                                                                             <  1:2095788/150‑1 (MQ=255)
                                                         ggAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATc                                                                                            >  1:2837173/1‑150 (MQ=255)
                                                                 cGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCa                                                                                    <  1:1267661/150‑1 (MQ=255)
                                                                   ggTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAtt                                                                                  >  2:1533888/1‑150 (MQ=255)
                                                                    gTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAttt                                                                                 >  2:2598116/1‑150 (MQ=255)
                                                                       aaaTATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCGCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAtttttt                                                                              >  2:289686/1‑150 (MQ=255)
                                                                       aaaTATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAtttttt                                                                              >  1:1777941/1‑150 (MQ=255)
                                                                       aaaTATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCAtttttt                                                                              <  1:24793/150‑1 (MQ=255)
                                                                             cAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGctct                                                                        <  2:2214446/150‑1 (MQ=255)
                                                                             cAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTATTACGctct                                                                        >  2:2860564/1‑150 (MQ=255)
                                                                                cAAAAAGCCCGCATCATTAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAc                                                                     <  1:2545276/150‑1 (MQ=255)
                                                                                cAAAAAGCCCGCATCATTAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAc                                                                     <  1:2544615/150‑1 (MQ=255)
                                                                                cAAAAAGCCCGCATCATGAATACTAGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAc                                                                     >  2:3880559/1‑150 (MQ=255)
                                                                                 aaaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCAcc                                                                    <  2:764145/150‑1 (MQ=255)
                                                                                  aaaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCg                                                                   >  1:2066725/1‑150 (MQ=255)
                                                                                   aaaGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCgg                                                                  >  2:2265410/1‑150 (MQ=255)
                                                                                      gCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCAt                                                               <  2:1137780/150‑1 (MQ=255)
                                                                                      gCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCAt                                                               <  2:4310622/150‑1 (MQ=255)
                                                                                      gCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCAt                                                               <  1:1945878/150‑1 (MQ=255)
                                                                                             tcatGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCAttt                                                        >  1:3059307/1‑150 (MQ=255)
                                                                                             tcatGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCAttt                                                        <  1:3166540/150‑1 (MQ=255)
                                                                                               atGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGc                                                      <  1:4004679/150‑1 (MQ=255)
                                                                                                tGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGcc                                                     >  2:2370096/1‑150 (MQ=255)
                                                                                                 gAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCAGTTGCca                                                    >  2:2665803/1‑150 (MQ=255)
                                                                                                 gAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCca                                                    <  2:1401125/150‑1 (MQ=255)
                                                                                                       tGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAat                                              >  1:2089757/1‑150 (MQ=255)
                                                                                                        ggATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAata                                             <  2:505604/150‑1 (MQ=255)
                                                                                                        ggATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAata                                             >  1:734825/1‑150 (MQ=255)
                                                                                                         gATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAatat                                            <  2:2070135/150‑1 (MQ=255)
                                                                                                          atatGAAGCATGAGAGTTAACTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATATTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATc                                           >  2:78564/1‑150 (MQ=255)
                                                                                                            atGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATccc                                         <  1:3972328/150‑1 (MQ=255)
                                                                                                               aaGCATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCccgcc                                      >  1:4083384/1‑150 (MQ=255)
                                                                                                                  cATGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTg                                   >  1:2184486/1‑150 (MQ=255)
                                                                                                                   aTGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGc                                  >  2:907711/1‑150 (MQ=255)
                                                                                                                    tGAGAGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGct                                 <  1:2010643/150‑1 (MQ=255)
                                                                                                                     gagaGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctc                                <  1:318098/150‑1 (MQ=255)
                                                                                                                     gagaGTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctc                                >  2:4231296/1‑150 (MQ=255)
                                                                                                                       gagTTACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGctcct                              <  2:2233983/150‑1 (MQ=255)
                                                                                                                           tACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGc                          >  1:1374638/1‑150 (MQ=255)
                                                                                                                           tACCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGc                          >  2:3954540/1‑150 (MQ=255)
                                                                                                                            aCCTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCa                         >  2:4273452/1‑150 (MQ=255)
                                                                                                                              cTCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAgc                       <  2:1480887/150‑1 (MQ=255)
                                                                                                                               tCAGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAgcg                      >  2:1171996/1‑150 (MQ=255)
                                                                                                                                 aGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCt                    <  1:3308208/150‑1 (MQ=255)
                                                                                                                                 aGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCt                    <  1:1629208/150‑1 (MQ=255)
                                                                                                                                 aGTGTTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCt                    >  1:3228568/1‑150 (MQ=255)
                                                                                                                                   tgtTTATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTc                  >  1:589768/1‑150 (MQ=255)
                                                                                                                                     tttATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCt                <  1:3557493/150‑1 (MQ=255)
                                                                                                                                     tttATATAAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCt                >  2:4047383/1‑150 (MQ=255)
                                                                                                                                       tatataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGt              >  2:2620377/1‑150 (MQ=255)
                                                                                                                                        atataAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTc             >  1:215346/1‑150 (MQ=255)
                                                                                                                                          ataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCa           <  1:4331064/150‑1 (MQ=255)
                                                                                                                                          ataAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCa           <  2:3361857/150‑1 (MQ=255)
                                                                                                                                           taAGGATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCAc          <  2:3605479/150‑1 (MQ=255)
                                                                                                                                              ggATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGaa       <  1:3348870/150‑1 (MQ=255)
                                                                                                                                               gATTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAAt      <  2:352960/150‑1 (MQ=255)
                                                                                                                                                aTTAGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATc     <  1:3824861/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AATATATCCAGAGAATAAAATCTGTCGCAGATAAGGTTGTATTAATAGTCTGTATCAGGAATGTTCGGGTTAAATATCAGCAAAAAGCCCGCATCATGAATACTGGATATGAAGCATGAGAGTTACCTCAGTGTTTATATAAGGATTCGGTCCCCCTCTCTGGAACGGTAACTCTCAATCTGATCGGTTCCTGCGTTAGTTCACATCACGACTCATTTTTTCGCTCTCACCGGCATCCCATTTGCCACAAAATATCCCGCCGTGCTCCTCGGCAGCGCTTCCTGTCCACGAATCATA  >  CP009273/321076‑321372

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: