Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,335,651 C→A 42.8% S7* (TCG→TAG)  kdsC → 3‑deoxy‑D‑manno‑octulosonate 8‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,335,6510CA42.8% 248.2 / inf 264S7* (TCG→TAG) kdsC3‑deoxy‑D‑manno‑octulosonate 8‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (68/83);  new base A (61/52);  total (129/135)
Fisher's exact test for biased strand distribution p-value = 1.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 267 total aligned reads displayed.

GAACTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATGTTC  >  CP009273/3335506‑3335800
                                                                                                                                                 |                                                                                                                                                     
gAACTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCtt                                                                                                                                                   >  1:2850583/1‑150 (MQ=255)
  aCTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGc                                                                                                                                                 >  2:1844560/1‑150 (MQ=255)
      aaTGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGAcc                                                                                                                                             <  1:4360699/150‑1 (MQ=255)
       aTGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCt                                                                                                                                            >  1:3546038/1‑150 (MQ=255)
            gTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTAc                                                                                                                                       <  1:128410/150‑1 (MQ=255)
              cccGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACgg                                                                                                                                     >  2:86160/1‑150 (MQ=255)
              cccGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACgg                                                                                                                                     <  1:2286452/150‑1 (MQ=255)
                cGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGcc                                                                                                                                   <  2:2018335/150‑1 (MQ=255)
                 gCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAATAAGGTGAGACGCTTTCGACCTGTTACGGccc                                                                                                                                  >  2:1238945/1‑150 (MQ=255)
                   cATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTg                                                                                                                                >  1:3430754/1‑150 (MQ=255)
                    atatCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGt                                                                                                                               >  1:1787972/1‑150 (MQ=255)
                          ccTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAgcgc                                                                                                                         >  1:454304/1‑150 (MQ=255)
                           cTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGcc                                                                                                                        <  1:2239778/150‑1 (MQ=255)
                              cGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGAc                                                                                                                     >  1:250041/1‑150 (MQ=255)
                                 gATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGtt                                                                                                                  >  2:256337/1‑150 (MQ=255)
                                 gATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGtt                                                                                                                  <  2:3737322/150‑1 (MQ=255)
                                 gATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGtt                                                                                                                  <  1:738043/150‑1 (MQ=255)
                                 gATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGtt                                                                                                                  <  2:3063/150‑1 (MQ=255)
                                   tGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTtat                                                                                                                >  1:3077376/1‑150 (MQ=255)
                                       tGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCa                                                                                                            <  2:3227983/150‑1 (MQ=255)
                                           gATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAg                                                                                                        <  1:2383160/150‑1 (MQ=255)
                                           gATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAg                                                                                                        <  2:2707406/150‑1 (MQ=255)
                                           gATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAg                                                                                                        >  1:3055828/1‑150 (MQ=255)
                                            aTGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGc                                                                                                       >  1:1098465/1‑150 (MQ=255)
                                            aTGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGc                                                                                                       <  2:1195070/150‑1 (MQ=255)
                                            aTGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGc                                                                                                       <  1:1924792/150‑1 (MQ=255)
                                              ggCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCag                                                                                                     >  2:1639032/1‑150 (MQ=255)
                                              ggCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCag                                                                                                     >  2:718129/1‑150 (MQ=255)
                                                    cATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACa                                                                                               >  1:2260432/1‑150 (MQ=255)
                                                     aTTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACAt                                                                                              <  1:2007595/150‑1 (MQ=255)
                                                      tttACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACAtt                                                                                             >  2:2076060/1‑150 (MQ=255)
                                                      tttACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACAtt                                                                                             >  2:290195/1‑150 (MQ=255)
                                                          cTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTc                                                                                         <  2:3318090/150‑1 (MQ=255)
                                                           tCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTct                                                                                        >  1:2422809/1‑150 (MQ=255)
                                                            cGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTctg                                                                                       >  2:1927760/1‑150 (MQ=255)
                                                              gtgtgtTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTctgct                                                                                     <  1:2074643/150‑1 (MQ=255)
                                                               tgtgtTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGTCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTctgctg                                                                                    >  2:1321234/1‑150 (MQ=255)
                                                               tgtgtTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTctgctg                                                                                    >  2:1321670/1‑150 (MQ=255)
                                                                   ttACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATcc                                                                                >  1:4213004/1‑150 (MQ=255)
                                                                    tACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCt                                                                               <  2:2546418/150‑1 (MQ=255)
                                                                    tACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCt                                                                               >  1:2745809/1‑150 (MQ=255)
                                                                      cATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCtcg                                                                             >  1:4476700/1‑150 (MQ=255)
                                                                       atatGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCtcga                                                                            >  2:491381/1‑150 (MQ=255)
                                                                               gATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCtcgatgtcgatg                                                                    <  1:4416841/150‑1 (MQ=255)
                                                                                  ttACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCg                                                                 <  1:1747964/150‑1 (MQ=255)
                                                                                     cTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTAc                                                              >  2:2996028/1‑150 (MQ=255)
                                                                                      tGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACt                                                             >  1:310944/1‑150 (MQ=255)
                                                                                       gCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTg                                                            <  2:230639/150‑1 (MQ=255)
                                                                                             aGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGAt                                                      >  1:1176730/1‑150 (MQ=255)
                                                                                             aGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGAt                                                      >  1:3146334/1‑150 (MQ=255)
                                                                                             aGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGAt                                                      >  1:3146343/1‑150 (MQ=255)
                                                                                              ggCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATg                                                     <  1:2641544/150‑1 (MQ=255)
                                                                                                cGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGc                                                   <  2:1700332/150‑1 (MQ=255)
                                                                                                cGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGc                                                   >  1:4284354/1‑150 (MQ=255)
                                                                                                 gTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGcc                                                  >  1:125878/1‑150 (MQ=255)
                                                                                                  tAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCt                                                 <  1:491381/150‑1 (MQ=255)
                                                                                                   aGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTg                                                <  1:324703/150‑1 (MQ=255)
                                                                                                   aGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTg                                                <  1:2623000/150‑1 (MQ=255)
                                                                                                     tgtAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGAt                                              <  1:3580457/150‑1 (MQ=255)
                                                                                                      gtAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGAtt                                             <  2:518264/150‑1 (MQ=255)
                                                                                                       tAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGAttt                                            <  1:1832374/150‑1 (MQ=255)
                                                                                                        aaaGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTta                                           <  1:1041437/150‑1 (MQ=255)
                                                                                                          aGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTtata                                         <  2:2509042/150‑1 (MQ=255)
                                                                                                          aGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTtata                                         >  2:111621/1‑150 (MQ=255)
                                                                                                          aGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTtata                                         <  1:3158776/150‑1 (MQ=255)
                                                                                                          aGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTtata                                         <  2:4474273/150‑1 (MQ=255)
                                                                                                           gATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTatat                                        >  1:3956861/1‑150 (MQ=255)
                                                                                                             ttCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATgg                                      <  1:812122/150‑1 (MQ=255)
                                                                                                              tCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATggg                                     <  1:2547658/150‑1 (MQ=255)
                                                                                                              tCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATggg                                     <  1:2164608/150‑1 (MQ=255)
                                                                                                                aaGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCa                                   <  2:4220949/150‑1 (MQ=255)
                                                                                                                 aGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCaa                                  >  1:3749336/1‑150 (MQ=255)
                                                                                                                 aGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCaa                                  >  2:2629339/1‑150 (MQ=255)
                                                                                                                  ggATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCaat                                 >  1:2556961/1‑150 (MQ=255)
                                                                                                                   gATAAACAACAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCaata                                <  2:3831498/150‑1 (MQ=255)
                                                                                                                       aacaacAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATgg                            >  2:2841853/1‑150 (MQ=255)
                                                                                                                       aacaacAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATgg                            >  2:2841994/1‑150 (MQ=255)
                                                                                                                       aacaacAATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATgg                            >  2:3666467/1‑150 (MQ=255)
                                                                                                                        acaacaATGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGc                           <  2:3129664/150‑1 (MQ=255)
                                                                                                                          aacaaTGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGa                         >  1:575088/1‑150 (MQ=255)
                                                                                                                          aacaaTGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGa                         >  1:3875163/1‑150 (MQ=255)
                                                                                                                          aacaaTGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGa                         >  1:3874033/1‑150 (MQ=255)
                                                                                                                            caaTGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAag                       <  2:253994/150‑1 (MQ=255)
                                                                                                                             aaTGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCGGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAaga                      <  2:3073014/150‑1 (MQ=255)
                                                                                                                              aTGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAagag                     >  1:896716/1‑150 (MQ=255)
                                                                                                                               tGAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGc                    >  2:2195912/1‑150 (MQ=255)
                                                                                                                                gAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCt                   <  1:2316895/150‑1 (MQ=255)
                                                                                                                                gAGCAAAGCAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCt                   >  1:1495774/1‑150 (MQ=255)
                                                                                                                                     aaGCAGTTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAg              >  2:1544684/1‑150 (MQ=255)
                                                                                                                                      aGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGc             <  2:3292431/150‑1 (MQ=255)
                                                                                                                                      aGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGc             <  2:2547791/150‑1 (MQ=255)
                                                                                                                                      aGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGc             <  2:215453/150‑1 (MQ=255)
                                                                                                                                        cAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGt           >  2:3557431/1‑150 (MQ=255)
                                                                                                                                        cAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGt           >  2:2159820/1‑150 (MQ=255)
                                                                                                                                        cAGGTGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGt           >  2:4403718/1‑150 (MQ=255)
                                                                                                                                           gTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCa        >  2:1592355/1‑150 (MQ=255)
                                                                                                                                            tGCGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCaa       <  1:1691658/150‑1 (MQ=255)
                                                                                                                                              cGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATg     >  1:2545784/1‑150 (MQ=255)
                                                                                                                                              cGTAGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATg     <  2:2936809/150‑1 (MQ=255)
                                                                                                                                                 cGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATGTTc  >  2:2346780/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GAACTTAATGCAGTCCCGCCATATCACCTCCGTGATGGTTGCCGATGGCGACCATTTACTCGGTGTGTTACATATGCATGATTTACTGCGTGCAGGCGTAGTGTAAAGATTCAAGGATAAACAACAATGAGCAAAGCAGGTGCGTCGCTTGCGACCTGTTACGGCCCTGTCAGCGCCGACGTTATAGCAAAAGCAGAGAACATTCGTCTGCTGATCCTCGATGTCGATGGCGTACTGTCAGATGGCCTGATTTATATGGGCAATAATGGCGAAGAGCTGAAAGCGTTCAATGTTC  >  CP009273/3335506‑3335800

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: