Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,373,522 C→T 5.0% intergenic (+85/‑134) epmA → / → yjeM Elongation Factor P Lys34 lysyltransferase/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,373,5220CT 5.0% 1203.9 / 34.8 319intergenic (+85/‑134)epmA/yjeMElongation Factor P Lys34 lysyltransferase/putative transporter
Reads supporting (aligned to +/- strand):  ref base C (175/128);  new base T (5/11);  total (180/139)
Fisher's exact test for biased strand distribution p-value = 4.19e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

Only 100 of 325 total aligned reads displayed.

TGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAA  >  CP009273/4373373‑4373671
                                                                                                                                                     |                                                                                                                                                     
tGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGc                                                                                                                                                       <  2:4852638/150‑1 (MQ=255)
tGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGc                                                                                                                                                       >  2:3652604/1‑150 (MQ=255)
   tGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCcac                                                                                                                                                    >  1:784322/1‑150 (MQ=255)
     ttGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACaa                                                                                                                                                  <  2:328949/150‑1 (MQ=255)
      tGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAAt                                                                                                                                                 >  1:3696101/1‑150 (MQ=255)
      tGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAAt                                                                                                                                                 <  1:2468655/150‑1 (MQ=255)
         cgcTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTa                                                                                                                                              <  1:1605395/150‑1 (MQ=255)
            tGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGa                                                                                                                                           >  1:394716/1‑150 (MQ=255)
             gggCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGAc                                                                                                                                          >  1:2781747/1‑150 (MQ=255)
             gggCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGAc                                                                                                                                          <  2:3013464/150‑1 (MQ=255)
                    gagaCACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGatat                                                                                                                                   <  1:362781/150‑1 (MQ=255)
                      gaCACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATAttt                                                                                                                                 <  2:5109551/150‑1 (MQ=255)
                        cacTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGa                                                                                                                               >  1:5098242/1‑150 (MQ=255)
                          cTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGata                                                                                                                             <  2:4812706/150‑1 (MQ=255)
                            ggCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATatc                                                                                                                           >  1:4389066/1‑150 (MQ=255)
                              cTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATatcat                                                                                                                         >  1:77476/1‑150 (MQ=255)
                                 aaGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCa                                                                                                                      <  1:226387/150‑1 (MQ=255)
                                 aaGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCa                                                                                                                      <  1:1701243/150‑1 (MQ=255)
                                  aGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAg                                                                                                                     >  2:506351/1‑150 (MQ=255)
                                   gTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAgg                                                                                                                    <  2:1208749/150‑1 (MQ=255)
                                   gTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAgg                                                                                                                    <  1:3540279/150‑1 (MQ=255)
                                    tCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGt                                                                                                                   <  1:4932140/150‑1 (MQ=255)
                                     cATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTa                                                                                                                  <  2:1863534/150‑1 (MQ=255)
                                      aTCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTAt                                                                                                                 <  2:821460/150‑1 (MQ=255)
                                      aTCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTAt                                                                                                                 <  2:4270438/150‑1 (MQ=255)
                                       tCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATc                                                                                                                <  2:1400506/150‑1 (MQ=255)
                                         gCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCaa                                                                                                              >  1:4149792/1‑150 (MQ=255)
                                             ttAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCt                                                                                                          <  1:3083467/150‑1 (MQ=255)
                                             ttAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCt                                                                                                          >  2:2130960/1‑150 (MQ=255)
                                               aGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGt                                                                                                        <  2:2384006/150‑1 (MQ=255)
                                                 cGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTAt                                                                                                      >  2:160010/1‑150 (MQ=255)
                                                   ttGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTg                                                                                                    <  1:3575244/150‑1 (MQ=255)
                                                   gtGACCGGGCATAATTCGGAAAAATGACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTg                                                                                                    <  2:4880679/149‑1 (MQ=255)
                                                     gACCGGGCATATTTCTGACAAATTACTGATTTAACATCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCAACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGt                                                                                                  <  2:4089171/150‑1 (MQ=255)
                                                       ccGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGttt                                                                                                <  1:3260076/150‑1 (MQ=255)
                                                        cGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACACTTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGtttt                                                                                               >  2:4984879/1‑150 (MQ=255)
                                                           gCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTAc                                                                                            <  2:2751785/150‑1 (MQ=255)
                                                           gCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTAc                                                                                            >  1:115083/1‑150 (MQ=255)
                                                             aTAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGc                                                                                          >  1:593315/1‑150 (MQ=255)
                                                               aaTTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGctg                                                                                        >  1:4313743/1‑150 (MQ=255)
                                                                aTTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGctgc                                                                                       >  1:4486566/1‑150 (MQ=255)
                                                                  tCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGctgcct                                                                                     <  1:3386634/150‑1 (MQ=255)
                                                                   cTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTAATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGctgcctg                                                                                    >  1:2435798/1‑150 (MQ=255)
                                                                      aaaaaTTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTc                                                                                 <  1:2948991/150‑1 (MQ=255)
                                                                        aaaTTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCaa                                                                               >  1:1165994/1‑150 (MQ=255)
                                                                          aTTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGTCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATg                                                                             <  1:4790934/150‑1 (MQ=255)
                                                                          aTTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATg                                                                             <  1:2203008/150‑1 (MQ=255)
                                                                          aTTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATg                                                                             <  2:4391869/150‑1 (MQ=255)
                                                                            tACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGtt                                                                           >  1:1028852/1‑150 (MQ=255)
                                                                            tACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGtt                                                                           >  2:581155/1‑150 (MQ=255)
                                                                             aCTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTg                                                                          >  2:4942875/1‑150 (MQ=255)
                                                                                    tAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCgcagtg                                                                   >  1:553792/1‑150 (MQ=255)
                                                                                      aCAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCgcagtgca                                                                 >  1:1278502/1‑150 (MQ=255)
                                                                                      aCAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCgcagtgca                                                                 <  2:3718097/150‑1 (MQ=255)
                                                                                      aCAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCgcagtgca                                                                 >  1:3278242/1‑150 (MQ=255)
                                                                                      aCAGAGAAGAATGGCGTGATAAGCGACATTCTGTATCATATGTTTTCATTTTTCCGTTAAAGGCCACTATTTAAGACTGGATATTTGCTATCATCAAGGTATGAATTATGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCgcattgca                                                                 <  2:1145801/150‑1 (MQ=255)
                                                                                       cAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCgcagtgcag                                                                <  2:2308250/150‑1 (MQ=255)
                                                                                            aagaaTGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAg                                                           >  1:2603637/1‑150 (MQ=255)
                                                                                            aagaaTGGCGTGATAAGCGCCATTATGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTAGTGGTTTTTACGCTGCATGCTCAATGTATCGCAGTCCAGTTCAg                                                           >  1:4692940/1‑150 (MQ=255)
                                                                                              gaaTGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGTCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTg                                                         <  2:1650658/150‑1 (MQ=255)
                                                                                               aaTGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTgg                                                        >  1:3301566/1‑150 (MQ=255)
                                                                                                   gCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCAc                                                    <  1:1004514/150‑1 (MQ=255)
                                                                                                    cGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACg                                                   >  2:278535/1‑150 (MQ=255)
                                                                                                     gTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACgt                                                  >  2:3574620/1‑150 (MQ=255)
                                                                                                     gTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTAAGCTGCCTGCTCAATGTTGAGCAGTGCAGTTCAGTGGGCACgt                                                  >  2:4205195/1‑150 (MQ=255)
                                                                                                        aTAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGtt                                               >  2:2341268/1‑150 (MQ=255)
                                                                                                        aTAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGtt                                               >  2:3470767/1‑150 (MQ=255)
                                                                                                           aGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTtata                                            >  1:4203024/1‑150 (MQ=255)
                                                                                                                ccttCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGAGTTCGACAcgc                                       >  1:446130/3‑150 (MQ=255)
                                                                                                                cATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACAcgc                                       <  1:3149246/150‑1 (MQ=255)
                                                                                                                 aTTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACAcgcg                                      >  1:516803/1‑150 (MQ=255)
                                                                                                                 aTTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACAcgcg                                      >  1:2986160/1‑150 (MQ=255)
                                                                                                                   tCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCt                                    >  1:1476512/1‑150 (MQ=255)
                                                                                                                     tGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGa                                  >  2:2208870/1‑150 (MQ=255)
                                                                                                                       tAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGaaa                                >  1:1990020/1‑150 (MQ=255)
                                                                                                                        aGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAAt                               >  2:3896134/1‑150 (MQ=255)
                                                                                                                          cATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGa                             >  2:4708335/1‑150 (MQ=255)
                                                                                                                          cATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGa                             <  2:1427437/150‑1 (MQ=255)
                                                                                                                            tatGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAg                           <  1:2289505/150‑1 (MQ=255)
                                                                                                                              tGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGa                         <  2:2250692/150‑1 (MQ=255)
                                                                                                                               gTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGAt                        >  2:3003895/1‑150 (MQ=255)
                                                                                                                               gTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGAt                        <  2:1030869/150‑1 (MQ=255)
                                                                                                                                ttttCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACTCGCTGAAATGAAGGATg                       >  2:13731/1‑150 (MQ=255)
                                                                                                                                 tttCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATgg                      >  2:1747235/1‑150 (MQ=255)
                                                                                                                                    cATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGttt                   >  2:2527364/1‑150 (MQ=255)
                                                                                                                                    cATTTGTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGAGATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGttt                   <  2:8801/150‑1 (MQ=255)
                                                                                                                                      tttttCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCa                 >  1:1576100/1‑150 (MQ=255)
                                                                                                                                       ttttCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCAt                <  2:3912644/150‑1 (MQ=255)
                                                                                                                                         ttCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATga              <  2:4054408/150‑2 (MQ=255)
                                                                                                                                            cGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTc           >  2:307575/1‑150 (MQ=255)
                                                                                                                                             gTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTca          >  1:720904/1‑150 (MQ=255)
                                                                                                                                             gTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTca          >  1:4499262/1‑150 (MQ=255)
                                                                                                                                             gTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTca          >  1:3809684/1‑150 (MQ=255)
                                                                                                                                              ttAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTcac         <  1:4441639/150‑1 (MQ=255)
                                                                                                                                               tAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTcaca        >  1:1490786/1‑150 (MQ=255)
                                                                                                                                                aaaGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCacac       <  2:1540582/150‑1 (MQ=255)
                                                                                                                                                aaaGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCacac       >  1:5001530/1‑150 (MQ=255)
                                                                                                                                                 aaGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACg      >  2:650409/1‑150 (MQ=255)
                                                                                                                                                  aGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGa     >  1:1331176/1‑150 (MQ=255)
                                                                                                                                                     ccACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATaa  <  1:2339045/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATGTTGGCGCTGGGCGCGGAGACACTGGCTGAAGTCATCGCCTTTAGCGTTGACCGGGCATAATTCTGAAAAATTACTGAATTAACAGCGAAGAATGGCGTGATAAGCGCCATTCTGTAGCATATGTTTTCATTTTTCCGTTAAAGGCCACAATTTAAGACTGGATATTTGATATCATCCAGGTATCAATTCTGTATTGGTTTTTACGCTGCCTGCTCAATGTTGCGCAGTGCAGTTCAGTGGGCACGTGTTATACACGCGCTGAAATGAAGGATGGTTTCATGCCTCACACGATAA  >  CP009273/4373373‑4373671

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: