Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,271,025 T→A 13.0% I373N (ATC→AAC)  yjcE → putative cation/proton antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,271,0250TA13.0% 629.5 / 71.0 200I373N (ATC→AAC) yjcEputative cation/proton antiporter
Reads supporting (aligned to +/- strand):  ref base T (92/82);  new base A (15/11);  total (107/93)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 205 total aligned reads displayed.

ATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACA  >  CP009273/4270876‑4271174
                                                                                                                                                     |                                                                                                                                                     
atatGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAt                                                                                                                                                       >  1:2764831/1‑150 (MQ=255)
atatGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAt                                                                                                                                                       >  2:4057486/1‑150 (MQ=255)
  atGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCa                                                                                                                                                     >  2:554929/1‑150 (MQ=255)
   tGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCAc                                                                                                                                                    <  2:4286840/150‑1 (MQ=255)
   tGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCAc                                                                                                                                                    >  2:2255014/1‑150 (MQ=255)
     cggcgCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCActc                                                                                                                                                  >  1:1809523/1‑150 (MQ=255)
           tGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCg                                                                                                                                            >  2:3864609/1‑150 (MQ=255)
           tGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCg                                                                                                                                            >  2:4069935/1‑150 (MQ=255)
              tGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGgtg                                                                                                                                         >  2:2527351/1‑150 (MQ=255)
                 tGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGc                                                                                                                                      <  2:3924058/150‑1 (MQ=255)
                 tGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGc                                                                                                                                      >  2:4047413/1‑150 (MQ=255)
                   gTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctc                                                                                                                                    >  1:1972368/1‑150 (MQ=255)
                   gTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctc                                                                                                                                    <  1:4039106/150‑1 (MQ=255)
                   gTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctc                                                                                                                                    <  2:209948/150‑1 (MQ=255)
                    tCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctct                                                                                                                                   <  1:2389092/150‑1 (MQ=255)
                    tCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctct                                                                                                                                   <  2:4814567/150‑1 (MQ=255)
                    tCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctct                                                                                                                                   <  2:648881/150‑1 (MQ=255)
                     ccGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctctc                                                                                                                                  <  1:2144150/150‑1 (MQ=255)
                     ccGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctctc                                                                                                                                  <  2:1266048/150‑1 (MQ=255)
                      cGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGctctct                                                                                                                                 <  1:4110312/150‑1 (MQ=255)
                       gTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTa                                                                                                                                <  2:1554647/150‑1 (MQ=255)
                           cGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATccc                                                                                                                            >  2:1596683/1‑150 (MQ=255)
                           cGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATccc                                                                                                                            >  1:380097/1‑150 (MQ=255)
                            ggctggctGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCg                                                                                                                           <  2:1174146/150‑1 (MQ=255)
                                  ctgTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCTGTGTGCTCTCTATCCCGCTGctc                                                                                                                     <  2:4169393/150‑1 (MQ=255)
                                   tgtgGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGctct                                                                                                                    <  2:3246177/150‑1 (MQ=255)
                                    gtgGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCtt                                                                                                                   <  2:2418512/150‑1 (MQ=255)
                                       gACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGcc                                                                                                                >  2:1961235/1‑150 (MQ=255)
                                          gATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGa                                                                                                             >  1:1261077/1‑150 (MQ=255)
                                            tGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATg                                                                                                           >  2:4483533/1‑150 (MQ=255)
                                             gAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATgg                                                                                                          >  1:848384/1‑150 (MQ=255)
                                             gAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATgg                                                                                                          >  1:1028699/1‑150 (MQ=255)
                                              aaaaaGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGt                                                                                                         <  1:3534659/150‑1 (MQ=255)
                                                 aaGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAAc                                                                                                      <  1:4628520/150‑1 (MQ=255)
                                                 aaGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAAc                                                                                                      <  1:3036469/150‑1 (MQ=255)
                                                  aGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACg                                                                                                     >  2:4431840/1‑150 (MQ=255)
                                                  aGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACg                                                                                                     >  1:5141013/1‑150 (MQ=255)
                                                      cAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCtt                                                                                                 <  2:4550054/150‑1 (MQ=255)
                                                       aGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTc                                                                                                >  1:3108075/1‑150 (MQ=255)
                                                          aaCCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCg                                                                                             >  2:5135974/1‑150 (MQ=255)
                                                                ttCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGgcgcgc                                                                                       >  1:4555108/1‑150 (MQ=255)
                                                                    tGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATg                                                                                   >  2:2585471/1‑150 (MQ=255)
                                                                       agaagaagCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGc                                                                                >  2:4831905/1‑150 (MQ=255)
                                                                       agaagaagCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGc                                                                                >  1:2659344/1‑150 (MQ=255)
                                                                         aagaagCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTg                                                                              >  1:2810216/1‑150 (MQ=255)
                                                                          agaagCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTgg                                                                             >  2:131322/1‑150 (MQ=255)
                                                                          agaagCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTgg                                                                             >  1:4330367/1‑150 (MQ=255)
                                                                           gaagCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGgt                                                                            >  2:4538615/1‑150 (MQ=255)
                                                                             agCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGgtgt                                                                          >  2:2156444/1‑150 (MQ=255)
                                                                               ccGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTc                                                                        <  1:2721174/150‑1 (MQ=255)
                                                                               ccGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCAGGGGTGCGTGGTGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCACGTCGCGCTATGAGCTGGTGTTc                                                                        >  1:581696/1‑150 (MQ=255)
                                                                                cGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTcc                                                                       >  1:4646700/1‑150 (MQ=255)
                                                                                cGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTcc                                                                       >  1:2048362/1‑150 (MQ=255)
                                                                                 gATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCt                                                                      >  1:2658646/1‑150 (MQ=255)
                                                                                 gATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCt                                                                      <  2:2762315/150‑1 (MQ=255)
                                                                                   tGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTgg                                                                    <  1:389306/150‑1 (MQ=255)
                                                                                    ggAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTggc                                                                   >  1:4435580/1‑150 (MQ=255)
                                                                                    ggAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTggc                                                                   <  2:1573942/150‑1 (MQ=255)
                                                                                     gAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTggcg                                                                  <  2:705569/150‑1 (MQ=255)
                                                                                     gAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTggcg                                                                  >  2:2791734/1‑150 (MQ=255)
                                                                                      aGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTggcgg                                                                 >  1:4837110/1‑150 (MQ=255)
                                                                                       gTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTggcggc                                                                >  2:1650292/1‑150 (MQ=255)
                                                                                       gTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTggcggc                                                                >  1:939159/1‑150 (MQ=255)
                                                                                       gTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTggcggc                                                                <  1:2053724/150‑1 (MQ=255)
                                                                                         ttGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTg                                                              <  2:3747037/150‑1 (MQ=255)
                                                                                          tGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTgg                                                             >  2:332494/1‑150 (MQ=255)
                                                                                           ggTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGc                                                            <  2:3184660/150‑1 (MQ=255)
                                                                                            gTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCg                                                           <  2:2774302/150‑1 (MQ=255)
                                                                                              tCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTc                                                         <  1:2527351/150‑1 (MQ=255)
                                                                                               cGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCa                                                        >  2:2468531/1‑150 (MQ=255)
                                                                                                     ccACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCAttctct                                                  >  2:1939394/1‑150 (MQ=255)
                                                                                                     ccACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCAttctct                                                  <  2:149786/150‑1 (MQ=255)
                                                                                                      cacaCGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCAttctctt                                                 >  1:420737/1‑150 (MQ=255)
                                                                                                      cacaCGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCAttctctt                                                 >  1:661504/1‑150 (MQ=255)
                                                                                                           gagaAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGc                                            >  2:805933/1‑150 (MQ=255)
                                                                                                            agaAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCt                                           >  1:2973255/1‑150 (MQ=255)
                                                                                                               aaTCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGtt                                        >  1:1089040/1‑150 (MQ=255)
                                                                                                                 tCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTg                                      <  2:4744943/150‑1 (MQ=255)
                                                                                                                    tGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCg                                   <  2:2087516/150‑1 (MQ=255)
                                                                                                                     gATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCgg                                  <  1:4710749/150‑1 (MQ=255)
                                                                                                                       tCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCg                                >  2:2668812/1‑150 (MQ=255)
                                                                                                                             cTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGa                          <  2:2283434/150‑1 (MQ=255)
                                                                                                                             cTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGa                          <  2:4866969/150‑1 (MQ=255)
                                                                                                                              tttCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGAt                         <  1:2385028/150‑1 (MQ=255)
                                                                                                                                tCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGt                       <  2:881947/150‑1 (MQ=255)
                                                                                                                                tCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGt                       <  1:1596683/150‑1 (MQ=255)
                                                                                                                                     ggggTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTCTGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCt                  >  1:4005622/1‑150 (MQ=255)
                                                                                                                                     ggggTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCt                  >  1:993525/1‑150 (MQ=255)
                                                                                                                                     ggggTGCGTGGGGCGATCACTCTGGCCGGGGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCGGGTTGTCGGCGTGGTGATGTTGCCt                  >  1:1164393/1‑150 (MQ=255)
                                                                                                                                     ggggTGCGTGGGGCGAACACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCt                  >  1:3452450/1‑150 (MQ=255)
                                                                                                                                        gTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTAtt               >  2:304904/1‑150 (MQ=255)
                                                                                                                                         tGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTc              <  2:2624058/150‑1 (MQ=255)
                                                                                                                                          gCGTGGGGCGATCACTCTGGCCGCTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCt             <  2:4634281/150‑1 (MQ=255)
                                                                                                                                           cGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTg            <  1:784391/150‑1 (MQ=255)
                                                                                                                                            gTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGc           >  2:3811815/1‑150 (MQ=255)
                                                                                                                                                 gCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTacaa      >  1:1387218/1‑150 (MQ=255)
                                                                                                                                                   gATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGGGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTacaaca    >  1:4437330/1‑150 (MQ=255)
                                                                                                                                                   gATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTacaaca    >  1:466720/1‑150 (MQ=255)
                                                                                                                                                   gATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTacaaca    >  1:470248/1‑150 (MQ=255)
                                                                                                                                                    aTCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAacac   <  1:1637247/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATATGCGGCGCTGATGCTGGTCCGTTTCGGCTGGCTGTGGACGATGAAAAAGTTCAGCAACCGCTTCCTGAAGAAGAAGCCGATGGAGTTTGGTTCGTGGACCACACGAGAAATCCTGATCGCGTCTTTCGCCGGGGTGCGTGGGGCGATCACTCTGGCCGGTGTGCTCTCTATCCCGCTGCTCTTGCCGGATGGTAACGTCTTCCCGGCGCGCTATGAGCTGGTGTTCCTGGCGGCTGGCGTCATTCTCTTCTCGCTGTTTGTCGGCGTGGTGATGTTGCCTATTCTGCTACAACACA  >  CP009273/4270876‑4271174

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: