Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,474,881 G→A 11.8% intergenic (+25/‑103) yfdQ → / → yfdR CPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,474,8810GA11.8% 195.8 / 15.3 68intergenic (+25/‑103)yfdQ/yfdRCPS‑53 (KpLE1) prophage; uncharacterized protein/CPS‑53 (KpLE1) prophage; conserved protein
Reads supporting (aligned to +/- strand):  ref base G (30/29);  new base A (4/4);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-01

GCGAAGACCAAAGACATTATGCCAGTGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGTCATCTTTCTCACTTCT  >  NC_000913/2474614‑2475116
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
gcgaagACCAAAGACATTATGCCAGTGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCAt                                                                                                                                                                                                                                      >  2:68264/1‑275 (MQ=255)
   aagACCAAAGACATTATGCCAGTGGCATTTGAATTTAAATGCGCTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTGGAGAAATTCAAAGAAAGCAAAGGAGAAACCTTTATTGGTCCTTTCACCGCCTGTTTTAANTACTTNNTATGCCCCCGCGGGGGCCTTTa                                                                                                                                                                                                                                   >  3:23244/1‑275 (MQ=255)
   aagACCAAAGACATTATGCCAGTGGCATTTGAATTTAAATGAGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTa                                                                                                                                                                                                                                   >  1:65816/1‑275 (MQ=255)
                ttATGCCAGTGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                                                                                                                                                                       <  1:70407/274‑1 (MQ=255)
                ttATGCCAGTGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCGCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                                                                                                                                                                       <  4:141727/274‑1 (MQ=255)
                         tGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTAt                                                                                                                                                                                                                                  >  4:189393/1‑254 (MQ=255)
                         tGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACg                                                                                                                                                                                                                           >  2:520523/1‑261 (MQ=255)
                         tGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCaata                                                                                                                                                                                                             >  4:70449/1‑275 (MQ=255)
                         tGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCaata                                                                                                                                                                                                             >  2:165292/1‑275 (MQ=255)
                         tGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCaata                                                                                                                                                                                                             >  4:52049/1‑275 (MQ=255)
                             aTTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGaa                                                                                                                                                                                                                              <  2:189394/254‑1 (MQ=255)
                             aTTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                                                                                                                                                                       <  4:520522/261‑1 (MQ=255)
                             aTTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCa                                                                                                                                                                                                         >  2:117137/1‑275 (MQ=255)
                              tttGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTAGGCCTCAGCATTATCACTGGCGATCGTCCTGTATTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCAc                                                                                                                                                                                                        <  4:191269/275‑1 (MQ=255)
                                ggAATTTAAATGCGTTCCGTTTGAAGGGCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCg                                                                                                                                                                                                      <  2:184898/274‑1 (MQ=255)
                                        aaTGCGTTCCTGTTTAAGGTCTGGAATAACGTCCGTTTAAATTACGCCTCAGCATTAGCACTGGCGATCGTCCTGTACGGGTTCTGCGCATTATTCAGCTTGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGCAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCTGGGGCATTTATGAAAACGTAATTAACTCAATAATCACCGGATGGTg                                                                                                                                                                                               <  4:165291/274‑1 (MQ=255)
                                            cGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGc                                                                                                                                                                                          <  4:405733/275‑1 (MQ=255)
                                             gTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCt                                                                                                                                                                                         <  2:70450/275‑1 (MQ=255)
                                                                  gAACGTCCGTTTAAATTACGCCTCATCATTATCACTGGCGATCGTCCGGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAgcg                                                                                                                                                                    <  1:23244/275‑1 (MQ=255)
                                                                  gAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAgcg                                                                                                                                                                    <  4:117136/275‑1 (MQ=255)
                                                                     cGTCCGTTTAAATTACGACTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCACCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAg                                                                                                                                                                                             >  2:382918/1‑247 (MQ=255)
                                                                         cGTTTAAATTTCGCCTCAGCATTATCACTGGCGTTCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACGCAATACTCACCGGATGGTGAGGGCt                                                                                                                                                                                         <  4:382917/247‑1 (MQ=255)
                                                                         cGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCgg                                                                                                                                                                                                     <  1:13500/235‑1 (MQ=255)
                                                                          gTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCa                                                                                                                                                             <  4:68263/274‑1 (MQ=255)
                                                                                      ccTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCa                                                                                                                                                                             >  2:521125/1‑246 (MQ=255)
                                                                                          aGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACt                                                                                                                                                                         <  4:521124/246‑1 (MQ=255)
                                                                                                                         ttCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCgg                                                                                                                                                                                                                                            <  6:582166/148‑1 (MQ=255)
                                                                                                                             gcgcATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAg                                                                                                         >  2:82239/1‑275 (MQ=255)
                                                                                                                                   tattCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGNGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCTGTGCACCNCATAGACGTGGAGAACAAAATGTCANATATTAAAAATTTTTACGGGAAGCAattt                                                                                                   >  3:96695/1‑275 (MQ=255)
                                                                                                                                             gAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGAc                                                                                          <  2:52050/274‑1 (MQ=255)
                                                                                                                                             gAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAAACAAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACa                                                                                         >  4:364622/1‑275 (MQ=255)
                                                                                                                                                    tGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGAt                                                                                                                                                                                                   >  1:50697/1‑162 (MQ=255)
                                                                                                                                                        ggAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTg                                                                                                                                                                                               <  3:50697/162‑1 (MQ=255)
                                                                                                                                                                tGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCa                                                                                                                                                                                                                                       <  7:582593/114‑1 (MQ=255)
                                                                                                                                                                tGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCa                                                                                                                                                                                                                                       >  8:582593/1‑114 (MQ=255)
                                                                                                                                                                   cTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACg                                                                                                                                                                                                                           >  1:98146/1‑123 (MQ=255)
                                                                                                                                                                     aaCGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCg                                                                                                                                                                                                                                             >  7:955844/1‑103 (MQ=255)
                                                                                                                                                                     aaCGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCg                                                                                                                                                                                                                                             <  8:955844/103‑1 (MQ=255)
                                                                                                                                                                       cGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                                                                                                                                                                       <  3:98146/123‑1 (MQ=255)
                                                                                                                                                                         aaTTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTaa                                                             <  1:96695/275‑1 (MQ=255)
                                                                                                                                                                                  aTCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCttt                                                                                                                                                                                   >  5:1094714/1‑148 (MQ=255)
                                                                                                                                                                                    cTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTa                                                                                                                                                                                 <  6:1094714/148‑1 (MQ=255)
                                                                                                                                                                                      gCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTttatta                                                                                              >  4:245256/1‑229 (MQ=255)
                                                                                                                                                                                       cTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAt                                                                                                                                                                                                                                      <  7:38510/92‑1 (MQ=255)
                                                                                                                                                                                       cTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAt                                                                                                                                                                                                                                      >  8:38510/1‑92 (MQ=255)
                                                                                                                                                                                          gttgAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGAc                                                                                          <  2:245257/229‑1 (MQ=255)
                                                                                                                                                                                                 aaTTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCaa                                      >  2:481691/1‑274 (MQ=255)
                                                                                                                                                                                                   ttCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCggggg                                                                                                                                                                                                                                         >  8:109518/1‑77 (MQ=255)
                                                                                                                                                                                                   ttCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCggggg                                                                                                                                                                                                                                         <  7:109518/77‑1 (MQ=255)
                                                                                                                                                                                                        aGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCa                                                                                                                                                                                                         <  7:750011/104‑1 (MQ=255)
                                                                                                                                                                                                        aGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCa                                                                                                                                                                                                         >  8:750011/1‑104 (MQ=255)
                                                                                                                                                                                                             gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                                                                                                                                                                       >  8:555003/1‑85 (MQ=255)
                                                                                                                                                                                                             gCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAAt                                                                                                                                                                                                                       <  7:555003/85‑1 (MQ=255)
                                                                                                                                                                                                                         ccTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtt                                                                                                                                                                                                                                     >  8:241476/1‑59 (MQ=38)
                                                                                                                                                                                                                         ccTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAtt                                                                                                                                                                                                                                     <  7:241476/59‑1 (MQ=38)
                                                                                                                                                                                                                            ttATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCtggg                                                                                                                                                                                                                                          >  6:362482/1‑51 (MQ=255)
                                                                                                                                                                                                                            ttATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCgggg                                                                                                                                                                                                                                          <  5:362482/51‑1 (MQ=255)
                                                                                                                                                                                                                             tATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAg                                                                                                                                                                                             >  2:289017/1‑95 (MQ=255)
                                                                                                                                                                                                                             tATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGTCATCTTtc          >  2:305089/1‑274 (MQ=255)
                                                                                                                                                                                                                                 ggTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCt                                                                                                                                                                                         <  4:289016/95‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGTCATCTTTCTCACTTCt  >  2:36649/1‑275 (MQ=255)
                                                                                                                                                                                                                                       ttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATc                                                                                                                                                                                                          >  6:1126152/1‑72 (MQ=255)
                                                                                                                                                                                                                                       ttCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATc                                                                                                                                                                                                          <  5:1126152/72‑1 (MQ=255)
                                                                                                                                                                                                                                           ccGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAt                                                                                                                                                                                                                                      >  7:816827/1‑40 (MQ=25)
                                                                                                                                                                                                                                           ccGCCTGATTTCATTACTGCAAATGCCCCCGCAGGGGCAt                                                                                                                                                                                                                                      <  8:816827/40‑1 (MQ=25)
                                                                                                                                                                                                                                             gCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATAc                                                                                                                                                   <  8:293433/121‑1 (MQ=255)
                                                                                                                                                                                                                                             gCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATAc                                                                                                                                                   >  7:293433/1‑121 (MQ=255)
                                                                                                                                                                                                                                                   tttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                                                                                                                                                            >  6:266142/1‑106 (MQ=255)
                                                                                                                                                                                                                                                   tttCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                                                                                                                                                            <  5:266142/106‑1 (MQ=255)
                                                                                                                                                                                                                                                        ttACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCt                                                                                                                                                                                     <  5:86921/76‑1 (MQ=255)
                                                                                                                                                                                                                                                        ttACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCt                                                                                                                                                                                     >  6:86921/1‑76 (MQ=255)
                                                                                                                                                                                                                                                         tACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                                                                                                                                                            >  7:311635/1‑100 (MQ=255)
                                                                                                                                                                                                                                                         tACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAg                                                                                                                                                            <  8:311635/100‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACaaa                                                                                                                                       >  6:1003057/1‑107 (MQ=255)
                                                                                                                                                                                                                                                                       cTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACaaa                                                                                                                                       <  5:1003057/107‑1 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                           
GCGAAGACCAAAGACATTATGCCAGTGGCATTTGAATTTAAATGCGTTCCGTTTGAAGGTCTGAAAGAACGTCCGTTTAAATTACGCCTCAGCATTATCACTGGCGATCGTCCTGTACTGGTTCTGCGCATTATTCAGCTGGAAGCGGTGCAGGAAGATATGGCTAACGAATTTCGTGATCTGCTTGTTGAGAAATTCAAAGACAGCAAAGTAGAAACCTTTATTGGTACTTTCACCGCCTGATTTCATTACTGCAAATGCCCCTGCGGGGGCATTTATGGAAACGTAATTAACTCAATAATCACCGGATGGTGAGGGCTTCCTTTTACCCAAACTCAGCGCGGTGCAGCGCATATACGTGGAGAACAAAATGTCATTTATTAAAACTTTTTCCGGGAAGCATTTTTATTATGACAAGATAAATAAAGACGACATCGTGATTAACGATATCGCGGTTTCCCTTTCAAATATCTGCCGCTTTGCCGGTCATCTTTCTCACTTCT  >  NC_000913/2474614‑2475116

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