Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,521,374 G→C 5.8% A216G (GCG→GGG)  yjhF ← KpLE2 phage‑like element; putative transporter YjhF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,521,3740GC 5.8% 290.7 / 6.6 86A216G (GCG→GGG) yjhFKpLE2 phage‑like element; putative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base G (39/42);  new base C (3/2);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

CACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATA  >  NC_000913/4521226‑4521522
                                                                                                                                                    |                                                                                                                                                    
cACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTcgc                                                                                                                                                     >  1:1090477/1‑150 (MQ=255)
  cGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCgcgc                                                                                                                                                   >  2:1488374/1‑150 (MQ=255)
   gTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGTTGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGcc                                                                                                                                                  <  1:760008/150‑1 (MQ=255)
       tAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCgaga                                                                                                                                              >  1:1257728/1‑150 (MQ=255)
        aaaaTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCgagag                                                                                                                                             <  2:1179704/150‑1 (MQ=255)
          aaTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGagagag                                                                                                                                           >  2:923182/1‑150 (MQ=255)
          aaTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACCCGCCgattcag                                                                                                                                           <  1:120490/150‑6 (MQ=255)
          aaTAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGCGACCCGCCgattcag                                                                                                                                           <  1:17237/150‑6 (MQ=255)
             aaGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCt                                                                                                                                        <  2:1291008/150‑1 (MQ=255)
                   gAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTc                                                                                                                                  <  2:574622/150‑1 (MQ=255)
                      tCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTcccc                                                                                                                               <  1:564635/150‑1 (MQ=255)
                       ccGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAACATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCct                                                                                                                              >  1:1445164/1‑150 (MQ=255)
                        cgcgGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCctc                                                                                                                             >  1:515354/1‑150 (MQ=255)
                        cgcgGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCctc                                                                                                                             <  1:1446067/150‑1 (MQ=255)
                        cgcgGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCGTTCGCTGAATCGGCGGGTCGCGCCGCGCGAGGCTATATTCCCctc                                                                                                                             >  2:1629320/1‑150 (MQ=255)
                         gcgGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCctct                                                                                                                            >  1:404565/1‑150 (MQ=255)
                          cgGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCtt                                                                                                                           <  1:971119/150‑1 (MQ=255)
                             cACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCttttt                                                                                                                        <  1:1470365/150‑1 (MQ=255)
                             cACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCttttt                                                                                                                        <  1:1476159/150‑1 (MQ=255)
                             cACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCttttt                                                                                                                        >  2:335428/1‑150 (MQ=255)
                             cACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCttttt                                                                                                                        <  2:357424/150‑1 (MQ=255)
                               cTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTa                                                                                                                      >  1:1436359/1‑150 (MQ=255)
                               cTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTa                                                                                                                      <  2:440058/150‑1 (MQ=255)
                                  ttGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAg                                                                                                                   <  2:266576/150‑1 (MQ=255)
                                   tGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAgg                                                                                                                  <  2:404565/150‑1 (MQ=255)
                                    ggCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGc                                                                                                                 <  1:898099/150‑1 (MQ=255)
                                    ggCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGc                                                                                                                 <  1:892395/150‑1 (MQ=255)
                                    ggCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGc                                                                                                                 >  1:390871/1‑150 (MQ=255)
                                      cATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGt                                                                                                               >  1:1278284/1‑150 (MQ=255)
                                      cATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGt                                                                                                               >  2:490707/1‑150 (MQ=255)
                                          aGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCgg                                                                                                           <  1:1318001/150‑1 (MQ=255)
                                          aGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCgg                                                                                                           <  1:1312251/150‑1 (MQ=255)
                                           gTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGa                                                                                                          >  1:270542/1‑150 (MQ=255)
                                           gTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGa                                                                                                          >  1:504505/1‑150 (MQ=255)
                                                cTGGATTATGGTGAGAACCGCCATCAGGAACACCGGCAGCATTGCAGTGAGAAAACTGATCCCAAACCCTGGCATTTCGCTTTCGCTGAATCGCCGGGTCGCGCCGAGACAGGCTATAATCCCCTCTATTTTAAACGAGTCCGtaataaa                                                                                                     >  2:223284/1‑149 (MQ=255)
                                                  ggATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATgc                                                                                                   >  1:1612373/1‑150 (MQ=255)
                                                      gatgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTg                                                                                               <  1:788986/150‑1 (MQ=255)
                                                       atgGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTgg                                                                                              <  1:1427571/150‑1 (MQ=255)
                                                           tGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGa                                                                                          >  2:141405/1‑150 (MQ=255)
                                                            gACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGaa                                                                                         >  2:1029682/1‑150 (MQ=255)
                                                             aCAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAAt                                                                                        >  2:338447/1‑150 (MQ=255)
                                                              cAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGGAGAAtt                                                                                       >  2:185173/1‑150 (MQ=255)
                                                               aaCCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAAttt                                                                                      <  1:1306333/150‑1 (MQ=255)
                                                                 ccGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTAt                                                                                    <  2:515354/150‑1 (MQ=255)
                                                                  cGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTAtt                                                                                   <  1:908392/150‑1 (MQ=255)
                                                                   gCCGTCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTc                                                                                  <  1:833926/150‑1 (MQ=255)
                                                                   gCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTc                                                                                  <  1:1447864/150‑1 (MQ=255)
                                                                   gCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTc                                                                                  >  1:1103720/1‑150 (MQ=255)
                                                                     cATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGACTCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCaa                                                                                <  1:1529216/150‑1 (MQ=255)
                                                                          ggATCACCGGCAGCATTGCCGTGAGAAAACTTATCCCAAAGCCTGGCATTTAGCTTTAGCTGAATCGGCGCGTCGCGCCGAGAGAGGCTATATTCCCCTCGTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTAGTCAATACCg                                                                           >  2:433748/1‑150 (MQ=255)
                                                                          ggATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCg                                                                           >  2:10415/1‑150 (MQ=255)
                                                                          ggATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCg                                                                           >  1:189789/1‑150 (MQ=255)
                                                                          ggATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCg                                                                           >  1:1211102/1‑150 (MQ=255)
                                                                           gATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCgg                                                                          <  2:898036/150‑1 (MQ=255)
                                                                                ccGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGc                                                                     >  1:250713/1‑150 (MQ=255)
                                                                                 cGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCa                                                                    >  1:1075985/1‑150 (MQ=255)
                                                                                      gcaTTGCCGTGGGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGaa                                                               >  2:513071/1‑150 (MQ=255)
                                                                                      gcaTTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGaa                                                               >  1:720370/1‑150 (MQ=255)
                                                                                           gCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAg                                                          <  1:1305739/150‑1 (MQ=255)
                                                                                           gCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCGTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAg                                                          >  1:251553/1‑150 (MQ=255)
                                                                                            ccGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAgg                                                         >  1:468982/1‑150 (MQ=255)
                                                                                             cGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGt                                                        >  1:1579489/1‑150 (MQ=255)
                                                                                               tGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTg                                                      <  2:1445164/150‑1 (MQ=255)
                                                                                               tGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTg                                                      <  1:319987/150‑1 (MQ=255)
                                                                                               tGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTg                                                      >  2:319987/1‑150 (MQ=255)
                                                                                                gagaAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTgg                                                     <  1:490707/150‑1 (MQ=255)
                                                                                                gagaAAACTGATCACAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCAGAGAGAGGCTATATTACCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTgg                                                     >  2:568683/1‑150 (MQ=255)
                                                                                                 agaAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGGCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGa                                                    <  2:1519200/150‑1 (MQ=255)
                                                                                                   aaaaCTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAATTAGTATGTAATTATGCGCTGGGAGATTTTATTCAATAGAGTGTCGTTAAGAATAATTGTTGGAAt                                                  <  2:1436359/150‑1 (MQ=255)
                                                                                                       ttGATCCCAAAGCGTGGCATTTCGGTTTCGCTGAATCGGCGGGTCGCGCCGAGAGTGGGTATATTGCGCTCTTTTTTAAAGGGGTCCGGAATAATGGGCTGGCAGAATTTATTGTATTCCGGGCCGCAAAGTATAATGGTTGGAATGCCg                                              <  2:270542/149‑1 (MQ=255)
                                                                                                            cccAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAAGAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAAt                                         >  1:914945/1‑150 (MQ=255)
                                                                                                              cAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGa                                       >  2:1067270/1‑150 (MQ=255)
                                                                                                               aaaGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGAt                                      <  2:1090091/150‑1 (MQ=255)
                                                                                                                       gCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGATTGATGCTATAATCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGtata                              <  1:335428/150‑1 (MQ=255)
                                                                                                                        cATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTatat                             <  2:30507/150‑1 (MQ=255)
                                                                                                                         aTTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATc                            <  2:1098663/150‑1 (MQ=255)
                                                                                                                         aTTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATc                            <  2:250713/150‑1 (MQ=255)
                                                                                                                           ttcgctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAg                          <  2:667808/150‑1 (MQ=255)
                                                                                                                           ttcgctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAg                          <  2:1612373/150‑1 (MQ=255)
                                                                                                                             cgctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTa                        >  2:55171/1‑150 (MQ=255)
                                                                                                                                  tcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGCGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAAt                   >  2:1155226/1‑150 (MQ=255)
                                                                                                                                    gctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATcc                 >  1:499625/1‑150 (MQ=255)
                                                                                                                                     ctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCt                <  1:1155226/150‑1 (MQ=255)
                                                                                                                                      tgtaTCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTa               <  1:923182/147‑1 (MQ=255)
                                                                                                                                       ctctCGGCGCGACCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTAc              >  1:146732/4‑150 (MQ=255)
                                                                                                                                        tctCGGCGCGACCCGCAGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATAACGTATATCAGTACAAATACTACa             >  2:336730/3‑150 (MQ=255)
                                                                                                                                         ctCGGCGCGACCCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAACCCTACAt            >  2:47352/2‑150 (MQ=255)
                                                                                                                                            ggcggGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAg         >  2:868684/1‑150 (MQ=255)
                                                                                                                                                   cgcgCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCata  >  1:714470/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAGCGCCATA  >  NC_000913/4521226‑4521522

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