Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,372,891 A→T 10.4% A144A (GCA→GCT epmA → Elongation Factor P Lys34 lysyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,372,8910AT10.4% 520.7 / 40.8 173A144A (GCA→GCTepmAElongation Factor P Lys34 lysyltransferase
Reads supporting (aligned to +/- strand):  ref base A (84/71);  new base T (6/12);  total (90/83)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 181 total aligned reads displayed.

GCCGCAGCTTCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACGCTGCTACAA  >  CP009273/4372746‑4373038
                                                                                                                                                 |                                                                                                                                                   
 ccGCAGCTTCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCTGaa                                                                                                                                                  <  1:431969/148‑1 (MQ=255)
 ccGCAGCTTCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGaa                                                                                                                                                  >  2:1968975/1‑148 (MQ=255)
  cGCAGCTTCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGaaa                                                                                                                                                 <  1:2552778/148‑1 (MQ=255)
      gCTTCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCTGAAAGCCt                                                                                                                                             <  2:103712/148‑1 (MQ=255)
       cTTCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCtt                                                                                                                                            >  1:2255842/1‑148 (MQ=255)
         tCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTc                                                                                                                                          >  2:1057288/1‑148 (MQ=255)
         tCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTc                                                                                                                                          >  1:1730421/1‑148 (MQ=255)
           cGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCTGAAAGCCTTTCtt                                                                                                                                        <  2:1832650/148‑1 (MQ=255)
             tAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTAt                                                                                                                                      >  2:143680/1‑148 (MQ=255)
               aTGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCTGAAAGCCTTTCTTATca                                                                                                                                    <  2:1454140/148‑1 (MQ=255)
                tGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATcaa                                                                                                                                   >  1:2446094/1‑148 (MQ=255)
                tGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATcaa                                                                                                                                   >  1:2446212/1‑148 (MQ=255)
                tGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATcaa                                                                                                                                   >  2:2205875/1‑148 (MQ=255)
                tGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATcaa                                                                                                                                   >  2:3262822/1‑148 (MQ=255)
                  aaGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATcaaca                                                                                                                                 >  2:2912898/1‑148 (MQ=255)
                   agagATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATcaacaa                                                                                                                                >  1:976996/1‑148 (MQ=255)
                    gagaTGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAg                                                                                                                               >  2:1670983/1‑148 (MQ=255)
                    gagaTGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAg                                                                                                                               <  1:1018869/148‑1 (MQ=255)
                    gagaTGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAg                                                                                                                               <  1:2673119/148‑1 (MQ=255)
                     agaTGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGc                                                                                                                              <  1:1881910/148‑1 (MQ=255)
                      gaTGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCt                                                                                                                             >  1:1721190/1‑148 (MQ=255)
                       aTGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCtt                                                                                                                            <  2:2073509/148‑1 (MQ=255)
                         ggggCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTc                                                                                                                          <  2:1298303/148‑1 (MQ=255)
                          gggCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCt                                                                                                                         >  1:1899965/1‑148 (MQ=255)
                          gggCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCt                                                                                                                         <  2:2963560/148‑1 (MQ=255)
                           ggCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCtt                                                                                                                        <  1:2156825/148‑1 (MQ=255)
                           ggCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCtt                                                                                                                        <  1:2831560/148‑1 (MQ=255)
                            gCGTTATCACAACCCTGTGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTg                                                                                                                       <  1:1452071/148‑1 (MQ=255)
                            gCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTg                                                                                                                       <  2:3300144/148‑1 (MQ=255)
                              gTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCg                                                                                                                     <  1:101123/148‑1 (MQ=255)
                               ttATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGt                                                                                                                    >  1:676167/1‑148 (MQ=255)
                                tATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGtt                                                                                                                   >  1:3167632/1‑148 (MQ=255)
                                   cacaACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATc                                                                                                                >  1:2220688/1‑148 (MQ=255)
                                   cacaACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTATTGCGTTATc                                                                                                                >  2:831716/1‑148 (MQ=255)
                                        cccTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGaa                                                                                                           <  2:3250385/148‑1 (MQ=255)
                                        cccTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGaa                                                                                                           <  2:214269/148‑1 (MQ=255)
                                            gAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTg                                                                                                       >  2:1207600/1‑148 (MQ=255)
                                            gAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTCGCGTTATCTGGAAATTg                                                                                                       <  1:2249652/148‑1 (MQ=255)
                                             aGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGa                                                                                                      >  2:803215/1‑148 (MQ=255)
                                                 cACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCg                                                                                                  >  2:1907266/1‑148 (MQ=255)
                                                   cTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGct                                                                                                >  1:1811403/1‑148 (MQ=255)
                                                    tATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGctc                                                                                               >  1:3298108/1‑148 (MQ=255)
                                                     aTGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGctct                                                                                              <  1:1330313/148‑1 (MQ=255)
                                                     aTGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGctct                                                                                              >  1:2183542/1‑148 (MQ=255)
                                                      tGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGctctc                                                                                             >  2:1774878/1‑148 (MQ=255)
                                                       gCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGctctct                                                                                            >  1:1519375/1‑148 (MQ=255)
                                                           gAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGTCGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCg                                                                                        >  1:2158612/1‑148 (MQ=255)
                                                           gAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCg                                                                                        >  2:1881391/1‑148 (MQ=255)
                                                            aGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGa                                                                                       <  1:1949852/148‑1 (MQ=255)
                                                              tGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACa                                                                                     >  1:505658/1‑148 (MQ=255)
                                                              tGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACa                                                                                     >  1:212582/1‑148 (MQ=255)
                                                               ggTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACaa                                                                                    >  2:3020595/1‑148 (MQ=255)
                                                               ggTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACaa                                                                                    >  2:2569285/1‑148 (MQ=255)
                                                                 tATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACaaaa                                                                                  <  1:875796/148‑1 (MQ=255)
                                                                  aTCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAAc                                                                                 >  1:1983769/1‑148 (MQ=255)
                                                                   tCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTTCCGACAAAACg                                                                                >  1:2746537/1‑148 (MQ=255)
                                                                   tCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTTCCGACAAAACg                                                                                >  1:1815273/1‑148 (MQ=255)
                                                                   tCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACg                                                                                >  2:3217864/1‑148 (MQ=255)
                                                                    cGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGc                                                                               >  2:1514352/1‑148 (MQ=255)
                                                                      aCCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCaa                                                                             <  2:3090797/148‑1 (MQ=255)
                                                                      aCCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCaa                                                                             <  2:2977239/148‑1 (MQ=255)
                                                                       ccGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAAc                                                                            >  1:1439395/1‑148 (MQ=255)
                                                                       ccGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAAc                                                                            >  1:1439400/1‑148 (MQ=255)
                                                                         gCACTATTATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGAGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTg                                                                          <  1:2227606/148‑1 (MQ=255)
                                                                           aCTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCg                                                                        <  2:2643351/148‑1 (MQ=255)
                                                                                gATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAATAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAg                                                                   <  2:435598/148‑1 (MQ=255)
                                                                                gATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAATAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAg                                                                   <  2:2099067/148‑1 (MQ=255)
                                                                                 atatGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCTGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGt                                                                  >  2:478133/1‑148 (MQ=255)
                                                                                     gTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCa                                                              <  1:1135766/148‑1 (MQ=255)
                                                                                      tACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAg                                                             <  1:45571/148‑1 (MQ=255)
                                                                                           gTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGaaa                                                        <  2:273197/148‑1 (MQ=255)
                                                                                             tgatgaACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACt                                                      <  2:1701286/148‑1 (MQ=255)
                                                                                                tgaACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGa                                                   >  2:1551665/1‑148 (MQ=255)
                                                                                                 gaACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGAt                                                  >  1:2835410/1‑148 (MQ=255)
                                                                                                  aaCGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGAtt                                                 >  1:238125/1‑148 (MQ=255)
                                                                                                         tGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCaa                                          <  2:1191073/148‑1 (MQ=255)
                                                                                                           gACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATg                                        <  2:676167/148‑1 (MQ=255)
                                                                                                              gATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTg                                     >  2:2385199/1‑148 (MQ=255)
                                                                                                              gATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTg                                     <  2:3167632/148‑1 (MQ=255)
                                                                                                               aTCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGc                                    <  1:3020595/148‑1 (MQ=255)
                                                                                                                tctcTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCt                                   >  1:721004/1‑148 (MQ=255)
                                                                                                                 ctctTACAACAGGTGCTGGACTGCCCGGCAGCTGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTg                                  <  1:2722502/148‑1 (MQ=255)
                                                                                                                 ctctTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTg                                  <  1:482305/148‑1 (MQ=255)
                                                                                                                   ctTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGAt                                <  2:92284/148‑1 (MQ=255)
                                                                                                                    ttACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATa                               <  2:1983769/148‑1 (MQ=255)
                                                                                                                     tACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATAc                              >  1:1786085/1‑148 (MQ=255)
                                                                                                                     tACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATAc                              >  1:999770/1‑148 (MQ=255)
                                                                                                                         acaGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCgaa                          <  2:1040292/148‑1 (MQ=255)
                                                                                                                          caGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCgaag                         >  2:2518287/1‑148 (MQ=255)
                                                                                                                          caGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCgaag                         <  2:567149/148‑1 (MQ=255)
                                                                                                                            ggTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCgaagaa                       <  1:2543290/148‑1 (MQ=255)
                                                                                                                              tGCTGGACTGCCCGGCAGCTGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAagaaga                     <  1:2326214/148‑1 (MQ=255)
                                                                                                                               gCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGAc                    <  1:441035/148‑1 (MQ=255)
                                                                                                                                        cccGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACg           <  1:1881391/148‑1 (MQ=255)
                                                                                                                                        cccGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACg           <  2:364108/148‑1 (MQ=255)
                                                                                                                                             cagcagAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACgctgct      >  1:1114207/1‑148 (MQ=255)
                                                                                                                                               gcagAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACGCTGCTAc    >  1:483356/1‑148 (MQ=255)
                                                                                                                                                cagAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACGCTGCTACa   >  1:1850538/1‑148 (MQ=255)
                                                                                                                                                cagAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACGCTGCTACa   >  1:2213638/1‑148 (MQ=255)
                                                                                                                                                cagAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACGCTGCTACa   >  1:402224/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCGCAGCTTCCGTAATGAAGAGATGGGGCGTTATCACAACCCTGAGTTCACTATGCTGGAGTGGTATCGACCGCACTATGATATGTACCGGTTGATGAACGAGGTGGACGATCTCTTACAACAGGTGCTGGACTGCCCGGCAGCAGAAAGCCTTTCTTATCAACAAGCTTTCTTGCGTTATCTGGAAATTGACCCGCTCTCTGCCGACAAAACGCAACTGCGGGAAGTCGCAGCGAAACTGGATTTGAGCAATGTTGCTGATACCGAAGAAGACCGCGACACGCTGCTACAA  >  CP009273/4372746‑4373038

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: