Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,820,188 A→T 5.6% intergenic (+122/+15) 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,1880AT 5.6% 655.1 / 19.0 195intergenic (+122/+15)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 T (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 = 9.78e-01

Only 100 of 198 total aligned reads displayed.

TCTAAGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGG  >  NC_002947/4820039‑4820337
                                                                                                                                                     |                                                                                                                                                     
 cTAAGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCaa                                                                                                                                                      >  1:3113333/1‑150 (MQ=255)
   aaGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgc                                                                                                                                                    >  2:5284336/1‑150 (MQ=255)
   aaGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgc                                                                                                                                                    <  1:3816806/150‑1 (MQ=255)
   aaGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAgc                                                                                                                                                    >  2:2717110/1‑150 (MQ=255)
           cGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCTTGGCGCCCCttt                                                                                                                                            <  2:893543/150‑3 (MQ=255)
           cGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCTTGGCGCCCCttt                                                                                                                                            <  2:2230216/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)
                 gCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATca                                                                                                                                      >  1:3701956/1‑150 (MQ=255)
                  cGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATcac                                                                                                                                     >  1:2791280/1‑150 (MQ=255)
                   gATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATcaca                                                                                                                                    >  1:907630/1‑150 (MQ=255)
                     ttattaGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGc                                                                                                                                  <  2:2030299/150‑1 (MQ=255)
                        ttaGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGcccc                                                                                                                               >  2:385926/1‑150 (MQ=255)
                             ccTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAg                                                                                                                          >  2:1445137/1‑150 (MQ=255)
                              cTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGt                                                                                                                         >  2:3422464/1‑150 (MQ=255)
                              cTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGt                                                                                                                         >  1:1658788/1‑150 (MQ=255)
                               tAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGtt                                                                                                                        <  2:276788/150‑1 (MQ=255)
                                 gCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTcc                                                                                                                      <  2:1048361/150‑1 (MQ=255)
                                   tttGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGc                                                                                                                    >  1:744547/1‑150 (MQ=255)
                                   tttGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGc                                                                                                                    >  2:2985477/1‑150 (MQ=255)
                                   tttGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGc                                                                                                                    >  1:4629275/1‑150 (MQ=255)
                                           ccTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGt                                                                                                            >  1:3612691/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)
                                               ttCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                        >  2:2488508/1‑150 (MQ=255)
                                               ttCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                        >  2:3349974/1‑150 (MQ=255)
                                               ttCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACt                                                                                                        >  1:1938718/1‑150 (MQ=255)
                                                  aCGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAg                                                                                                     >  1:5044075/1‑150 (MQ=255)
                                                             cAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCgg                                                                                                 >  1:3132629/1‑143 (MQ=255)
                                                             cAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCgg                                                                                                 <  2:3132629/143‑1 (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)
                                                                          cTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAg                                                                             <  1:1183895/150‑1 (MQ=255)
                                                                          cTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAg                                                                             <  2:1570638/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)
                                                                             ccTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTc                                                                                               >  2:2016708/1‑129 (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)
                                                                                     atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt                                                                  <  2:3200597/150‑1 (MQ=255)
                                                                                     atCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCt                                                                  <  1:3365975/150‑1 (MQ=255)
                                                                                      tCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACGCGACCAATGAAAGGGGCGCCAAGTGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTg                                                                 >  2:2669630/1‑150 (MQ=255)
                                                                                      tCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTg                                                                 >  1:391814/1‑150 (MQ=255)
                                                                                       cTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGCCCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                <  2:4626569/150‑1 (MQ=255)
                                                                                       cTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                <  2:1887007/150‑1 (MQ=255)
                                                                                       cTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                >  2:2106196/1‑150 (MQ=255)
                                                                                       cTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                <  2:4137296/150‑1 (MQ=255)
                                                                                        tGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCa                                                               >  1:3748558/1‑150 (MQ=255)
                                                                                        tGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCa                                                               >  1:1007248/1‑150 (MQ=255)
                                                                                         gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                >  2:4411994/1‑148 (MQ=255)
                                                                                         gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGc                                                                <  1:4411994/148‑1 (MQ=255)
                                                                                         gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAg                                                              >  1:302379/1‑150 (MQ=255)
                                                                                         gCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAg                                                              >  2:3121301/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                                                       <  2:4865776/150‑1 (MQ=255)
                                                                                                     gAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATa                                                  >  1:4450698/1‑150 (MQ=255)
                                                                                                        cGACTCAGCCGCAGTGTGCGCAACCCGACCAATGATAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTc                                               <  1:3727673/150‑1 (MQ=255)
                                                                                                         gACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCg                                              <  2:470809/150‑1 (MQ=255)
                                                                                                         gACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCg                                              <  1:2350187/150‑1 (MQ=255)
                                                                                                               gCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                        <  1:353956/150‑1 (MQ=255)
                                                                                                                   cAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCt                                    <  2:1151311/150‑1 (MQ=255)
                                                                                                                      tgtgCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCa                                 <  2:788235/150‑1 (MQ=255)
                                                                                                                        tgCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCAcc                               <  1:2717110/150‑1 (MQ=255)
                                                                                                                          cgcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGt                             >  1:4174357/1‑150 (MQ=255)
                                                                                                                           gcAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGtg                            <  2:3574267/150‑1 (MQ=255)
                                                                                                                                ccGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCa                       >  1:2517294/1‑150 (MQ=255)
                                                                                                                                   aCCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTc                    <  2:2524922/150‑1 (MQ=255)
                                                                                                                                     gatgggcaGGGGCGCTTGGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGa                  >  1:3778649/8‑150 (MQ=255)
                                                                                                                                     cAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGa                  <  1:4909721/150‑1 (MQ=255)
                                                                                                                                     cAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGa                  >  2:3541108/1‑150 (MQ=255)
                                                                                                                                      aaTGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGaa                 >  2:4041115/1‑150 (MQ=255)
                                                                                                                                        tGATAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAAc                                        <  2:1806500/125‑1 (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:3951635/3‑150 (MQ=255)
                                                                                                                                          aaaGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATc             <  1:3736497/150‑1 (MQ=255)
                                                                                                                                           caGGGGCGCTTGGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGccc                                                           >  2:3879304/2‑103 (MQ=255)
                                                                                                                                             ggggCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGcc          <  2:2491852/150‑1 (MQ=255)
                                                                                                                                               ggCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCg                   >  1:5155415/1‑139 (MQ=255)
                                                                                                                                                gcgcCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGc       >  1:1529222/1‑150 (MQ=255)
                                                                                                                                                gcgcCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGc       >  1:429648/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)
                                                                                                                                                     |                                                                                                                                                     
TCTAAGGGAACCGTTTGGCGATTATTAGCCCTAGCTTTGAGGGCCTCTTCACGAAAATAGCCAGCAACGAGCGACTGCCTGATTGATCTGCAAATCCGGATGAACGACTCAGCCGCAGTGTGCGCAACCCGACCAATGAAAGGGGCGCCAAGCGCCCCTGCCCATCACAGCCCCTCCAGTTCCGCCATCAGGTCACTGAGCCGGTCGACCTTGTCATCATCCAGCGCACTGCCCTGCAGCCCCTGTACATACTCGGCCAACGCCTCCACCGTGCTGCACTCGAACATCGCCCGCAGCGG  >  NC_002947/4820039‑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: