Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,521,370 G→C 6.5% T217T (ACC→ACG) ‡ 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,3700GC 6.5% 279.8 / 6.4 87T217T (ACC→ACG) ‡yjhFKpLE2 phage‑like element; putative transporter YjhF
Reads supporting (aligned to +/- strand):  ref base G (39/42);  new base C (4/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAG  >  NC_000913/4521226‑4521515
                                                                                                                                                |                                                                                                                                                 
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:504505/1‑150 (MQ=255)
                                           gTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGa                                                                                                   >  1:270542/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:250713/150‑1 (MQ=255)
                                                                                                                         aTTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATc                     <  2:1098663/150‑1 (MQ=255)
                                                                                                                           ttcgctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAg                   <  2:1612373/150‑1 (MQ=255)
                                                                                                                           ttcgctttcgctGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAg                   <  2:667808/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)
                                                                                                                                                |                                                                                                                                                 
CACGTTGTAAAATAAGCCTGAATCCGCGGCACTTTTGGCATGAGTCATCTGGATGATGGTGACAACCGCCATCAGGATCACCGGCAGCATTGCCGTGAGAAAACTGATCCCAAAGCCTGGCATTTCGCTTTCGCTGAATCGGCGGGTCGCGCCGAGAGAGGCTATATTCCCCTCTTTTTTAAAGGCGTCCGGAATAATGCGCTGGCAGAATTTATTCAATACCGGGCCGCAAAGAATAAAGGTTGGAATGCCGACAATGATACCGTATATCAGTACCAATCCTACATCAG  >  NC_000913/4521226‑4521515

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: ‑