Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,820,189 A→G 5.6% intergenic (+123/+14) PP_4242 → / ← pvdL conserved protein of unknown function with Rhs repeat‑associated core domain/non‑ribosomal peptide synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,820,1890AG 5.6% 620.9 / 15.9 195intergenic (+123/+14)PP_4242/pvdLconserved protein of unknown function with Rhs repeat‑associated core domain/non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base A (93/91);  new base G (5/6);  total (98/97)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

Only 100 of 197 total aligned reads displayed.

CTAAGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGG  >  NC_002947/4820040‑4820337
                                                                                                                                                     |                                                                                                                                                    
cTAAGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCaa                                                                                                                                                      >  1:3113333/1‑150 (MQ=255)
  aaGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgc                                                                                                                                                    >  2:5284336/1‑150 (MQ=255)
          cttttGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGc                                                                                                                                            <  1:3671031/148‑1 (MQ=255)
          cGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCTTGGCGCCCCttt                                                                                                                                            <  2:893543/150‑3 (MQ=255)
          cGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGc                                                                                                                                            <  2:2897219/150‑1 (MQ=255)
           gTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGcc                                                                                                                                           >  1:3053479/1‑150 (MQ=255)
           gTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGAGGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGcc                                                                                                                                           <  1:4635017/150‑1 (MQ=255)
            tttGGCGATTATTAGCCCTAGATTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGccc                                                                                                                                          <  1:122425/150‑1 (MQ=255)
              tGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCAt                                                                                                                                        <  2:5044116/150‑1 (MQ=255)
              tGGCGATTATTAGCCCTAGCTTTGACGGCCTCTTCACGAAAATAGCCACCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCAt                                                                                                                                        <  2:4392084/150‑1 (MQ=255)
               ggCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCTTGGCGCCCCtttcattg                                                                                                                                       <  1:614232/150‑8 (MQ=255)
                 cGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATcac                                                                                                                                     >  1:2791280/1‑150 (MQ=255)
                       ttaGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGcccc                                                                                                                               >  2:385926/1‑150 (MQ=255)
                         aGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTc                                                                                                                             <  2:2111001/150‑1 (MQ=255)
                         aGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTc                                                                                                                             <  1:3739193/150‑1 (MQ=255)
                           cccTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCa                                                                                                                           <  2:1728477/150‑1 (MQ=255)
                           cccTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCa                                                                                                                           <  2:2017690/150‑1 (MQ=255)
                             cTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGt                                                                                                                         >  2:3422464/1‑150 (MQ=255)
                             cTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGt                                                                                                                         >  1:1658788/1‑150 (MQ=255)
                              tAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGtt                                                                                                                        <  2:276788/150‑1 (MQ=255)
                                  tttGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGc                                                                                                                    >  2:2985477/1‑150 (MQ=255)
                                  tttGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGc                                                                                                                    >  1:4629275/1‑150 (MQ=255)
                                  tttGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGc                                                                                                                    >  1:744547/1‑150 (MQ=255)
                                            tctTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAAACCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGATCCGCCATCAGGTCa                                                                                                          >  2:3718543/1‑150 (MQ=255)
                                             ctTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCAc                                                                                                         <  2:5057294/150‑1 (MQ=255)
                                             ctTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCAc                                                                                                         >  1:5077970/1‑150 (MQ=255)
                                             ctTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCAc                                                                                                         <  2:2674747/150‑1 (MQ=255)
                                              ttCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                        >  2:3349974/1‑150 (MQ=255)
                                              ttCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                        >  1:1938718/1‑150 (MQ=255)
                                                 aCGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAg                                                                                                     >  1:5044075/1‑150 (MQ=255)
                                                       aTAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTc                                                                                               <  1:3763462/150‑1 (MQ=255)
                                                       aTAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTc                                                                                               <  1:3712249/150‑1 (MQ=255)
                                                           ccAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGAcc                                                                                           <  1:4293984/150‑1 (MQ=255)
                                                            cAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCgg                                                                                                 <  2:3132629/143‑1 (MQ=255)
                                                            cAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCgg                                                                                                 >  1:3132629/1‑143 (MQ=255)
                                                             aGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCCGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCtt                                                                                         >  2:4956292/1‑150 (MQ=255)
                                                                aaCGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtc                                                                                      <  1:2276953/150‑1 (MQ=255)
                                                                  cGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtcat                                                                                    <  1:82626/150‑1 (MQ=255)
                                                                   gAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCCCTGAACGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGtcatc                                                                                   <  2:4549491/150‑1 (MQ=255)
                                                                       gACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATcc                                                                               >  2:2616747/1‑150 (MQ=255)
                                                                         cTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAg                                                                             <  1:1183895/150‑1 (MQ=255)
                                                                         cTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAg                                                                             <  1:1454792/150‑1 (MQ=255)
                                                                          tGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAgc                                                                            >  1:5382902/1‑150 (MQ=255)
                                                                            ccTTATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                        <  1:11925/120‑1 (MQ=255)
                                                                            ccTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                        >  2:11925/1‑120 (MQ=255)
                                                                             cTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCa                                                                         <  1:4652250/150‑1 (MQ=255)
                                                                              tgattgatGTGCAAATCCGGATGAACGACTCTGCCGCAGTGTGCGCTACGCGACCAATGAAAGGGGCGGCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCAGTGAGCCGGTCGAACTTGTCATCATGCAGCGCAc                                                                        >  2:443988/1‑150 (MQ=255)
                                                                                attgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGccc                                                                                                                                                >  1:5349854/1‑76 (MQ=255)
                                                                                attgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGccc                                                                                                                                                <  2:5349854/76‑1 (MQ=255)
                                                                                attgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGt                                                                                                            <  2:307378/112‑1 (MQ=255)
                                                                                attgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTg                                                                      >  2:2493514/1‑150 (MQ=255)
                                                                                attgatCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTg                                                                      >  2:1523330/1‑150 (MQ=255)
                                                                                    atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt                                                                  <  1:3365975/150‑1 (MQ=255)
                                                                                    atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt                                                                  <  2:3200597/150‑1 (MQ=255)
                                                                                      cTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGCCCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                <  2:4626569/150‑1 (MQ=255)
                                                                                      cTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                <  2:4137296/150‑1 (MQ=255)
                                                                                       tGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCa                                                               >  1:1007248/1‑150 (MQ=255)
                                                                                        gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                <  1:4411994/148‑1 (MQ=255)
                                                                                        gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                >  2:4411994/1‑148 (MQ=255)
                                                                                        gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAg                                                              >  1:302379/1‑150 (MQ=255)
                                                                                        gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAg                                                              >  2:3121301/1‑150 (MQ=255)
                                                                                        gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAg                                                              >  2:1821279/1‑150 (MQ=255)
                                                                                          aaaTCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGcc                                                            <  1:598334/150‑1 (MQ=255)
                                                                                            aTCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGcccc                                                          <  1:4958844/150‑1 (MQ=255)
                                                                                             tCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCAAGCGCACTGCCCTGCAGCCCCt                                                         >  1:3234209/1‑150 (MQ=255)
                                                                                               cGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACCGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGt                                                       <  2:5104897/150‑1 (MQ=255)
                                                                                               cGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGt                                                       <  1:4032946/150‑1 (MQ=255)
                                                                                               cGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGt                                                       <  2:4865776/150‑1 (MQ=255)
                                                                                                    gAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATa                                                  >  1:4450698/1‑150 (MQ=255)
                                                                                                        gACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCg                                              <  2:470809/150‑1 (MQ=255)
                                                                                                        gACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCg                                              <  1:2350187/150‑1 (MQ=255)
                                                                                                          cTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGc                                            <  2:1944091/150‑1 (MQ=255)
                                                                                                              gCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                        <  1:353956/150‑1 (MQ=255)
                                                                                                                     tgtgCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCa                                 <  2:788235/150‑1 (MQ=255)
                                                                                                                       tgCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCAcc                               <  1:2717110/150‑1 (MQ=255)
                                                                                                                         cgcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGt                             >  2:2732270/1‑150 (MQ=255)
                                                                                                                         cgcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGt                             >  2:1094646/1‑150 (MQ=255)
                                                                                                                          gcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGtg                            <  2:3574267/150‑1 (MQ=255)
                                                                                                                          gcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGtg                            <  1:4350109/150‑1 (MQ=255)
                                                                                                                              cccGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGtgctgc                        >  1:373402/1‑150 (MQ=255)
                                                                                                                               ccGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCa                       >  1:2517294/1‑150 (MQ=255)
                                                                                                                                 gACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACt                     >  2:2044962/1‑150 (MQ=255)
                                                                                                                                    cAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGa                  >  2:3541108/1‑150 (MQ=255)
                                                                                                                                    cAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGa                  <  1:4909721/150‑1 (MQ=255)
                                                                                                                                    cAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGa                  >  1:3779342/1‑150 (MQ=255)
                                                                                                                                    cAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGa                  >  2:2315376/1‑150 (MQ=255)
                                                                                                                                     aaTGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGaa                 >  2:4041115/1‑150 (MQ=255)
                                                                                                                                       tGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                        >  1:1806500/1‑125 (MQ=255)
                                                                                                                                       tGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACa               >  1:2867521/1‑150 (MQ=255)
                                                                                                                                        gAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACAt              >  2:5398636/1‑150 (MQ=255)
                                                                                                                                         gcaGGGGCGCTTGGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATc             >  1:3951519/3‑150 (MQ=255)
                                                                                                                                          caGGGGCGCTTGGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGccc                                                           >  2:3879304/2‑103 (MQ=255)
                                                                                                                                          caGGGGCGCTTGGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGccc                                                           <  1:3879304/102‑1 (MQ=255)
                                                                                                                                            ggggCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGcc          <  2:2491673/150‑1 (MQ=255)
                                                                                                                                              ggCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCg                   >  1:5155415/1‑139 (MQ=255)
                                                                                                                                               gcgcCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGc       >  1:929172/1‑150 (MQ=255)
                                                                                                                                               gcgcCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGc       >  1:1529222/1‑150 (MQ=255)
                                                                                                                                                 gcCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGtgctgc                        >  1:4693214/1‑131 (MQ=255)
                                                                                                                                                   cAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCg   <  2:3612691/150‑1 (MQ=255)
                                                                                                                                                    aaGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCgg  >  1:1878870/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CTAAGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGG  >  NC_002947/4820040‑4820337

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: