Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,590,430 T→C intergenic (‑429/+424) ydeT ← / ← hipA fimbrial usher domain‑containing protein YdeT/serine/threonine‑protein kinase toxin HipA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,590,4300TC97.8% 161.1 / ‑4.0 46intergenic (‑429/+424)ydeT/hipAfimbrial usher domain‑containing protein YdeT/serine/threonine‑protein kinase toxin HipA
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (25/20);  total (25/21)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CAGAATATTATTGCGAACAAGGGAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCGAATTGCTCTAAACAGATAAGTGCAGGGTTTATTAAAGCTGTCCGACAGACTTTCTTAACCTGGTAATATCTCTCTTGGGTGACTCTCCAAACATCCGCGAATATTCCCGACTGAAATGGGATAGCTTTCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAATTATGA  >  NC_000913/1590190‑1590673
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
cctgagcgataccgtaggcgatggatagtttgctggcaacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGAc                                                                                                                                                                                                                                           >  1:249596‑M2/211‑251 (MQ=255)
      cgataccgtaggcgatggctagtttgctggcaacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCAt                                                                                                                                                                                                                                      >  2:409972‑M2/205‑250 (MQ=255)
         taccgtaggcgatggatagtttgctggcaacaccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaaccgataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTACCGACTGAAATGGGTTAGCTTCCGTAACCGACAGCATAGgc                                                                                                                                                                                                                                  >  1:216876‑M2/202‑251 (MQ=255)
                                       caccggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGCCTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGttcattc                                                                                                                                                                                                      >  1:162918‑M2/172‑249 (MQ=255)
                                        accggttaccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgattgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGc                                                                                                                                                                                                   >  1:173492‑M2/171‑251 (MQ=255)
                                               accagaatgcgcttaccggaaagaaaacccatagctttaatccttcttgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTc                                                                                                                                                                                            >  2:26293‑M2/164‑251 (MQ=255)
                                               accagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgagtcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagtggccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTc                                                                                                                                                                                            <  1:388306‑M2/88‑1 (MQ=255)
                                                ccagaatgcgcttaccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaagcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAAGATTCCCGACTGAAATGGGATAGCTGCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCAGTCATCATCTGTCg                                                                                                                                                                                           <  2:249596‑M2/89‑1 (MQ=255)
                                                             accggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCAtt                                                                                                                                                                               <  1:319418‑M2/101‑1 (MQ=255)
                                                              ccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGttcat                                                                                                                                                                                                        <  2:84087‑M2/76‑1 (MQ=255)
                                                              ccggaaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGttcat                                                                                                                                                                                                        >  1:84087‑M2/149‑224 (MQ=255)
                                                                  aaagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACg                                                                                                                                                                         >  1:68112‑M2/145‑251 (MQ=255)
                                                                   aagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGCCTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCGGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGt                                                                                                                                                                        <  1:140637‑M2/108‑1 (MQ=255)
                                                                   aagaaaacccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGt                                                                                                                                                                        <  1:306200‑M2/108‑1 (MQ=255)
                                                                          cccatagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCAtt                                                                                                                                                                 <  1:237345‑M2/115‑1 (MQ=255)
                                                                            catagctttaatccttattgttgatgcttgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTc                                                                                                                                                               <  1:370814‑M2/117‑1 (MQ=255)
                                                                                     aatccttattgttgatgctggttgtgccttaaaatcaggcgattcgtcgttttaggaaacagtacgaacagagaaacggttattataagcaacctggctgtgagtagttatagttgccagtgccgCGTTTATTCCCGACTGAAATGGGATAGCTTCCGTAGCCGACAGTATAGGCTGCTGTGGTGACATCGTAATGTTCATTCATCATGTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCg                                                                                                                                                       <  1:199015‑M2/125‑1 (MQ=255)
                                                                                                        tgttgtgcctgaaaatcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTACATTCAGCATCTGTCGCCTGACTTCATGGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGc                                                                                                                                     <  1:178737‑M2/143‑1 (MQ=255)
                                                                                                             tgccggaaaatcaggcgattcttcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGGGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGGCGTCAGTGAGACGTAGCCATTGCTGATACCGGAGTGGACGCATTCCCGCTTGATc                                                                                                                               <  2:30583‑M2/149‑1 (MQ=255)
                                                                                                                ctgaaactcaggcgattcgtcgttttagtaaacagtacgaacagataaacggttcttataatcaacctggctgtgagtagctatagttgccaggaccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCAATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTc                                                                                                                           >  2:349812‑M2/99‑251 (MQ=255)
                                                                                                                        caggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATa                                                                                                                   <  1:409972‑M2/161‑1 (MQ=255)
                                                                                                                          ggcgattcgtcgttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGTGGACATCGTAATGTTAATTCAGCATCTGTCGCCTGACTTCATTGGGACGTAGCCATTTCTGATACGGGAGGGAACCCATTCCCGCTTGATCCGCCcttaatta                                                                                                                   >  1:160883‑M2/89‑246 (MQ=255)
                                                                                                                                      ttttagtaaacagtacgaacagataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCct                                                                                                     >  2:3055‑M2/77‑251 (MQ=255)
                                                                                                                                             aaactgtacgaacagataaacggttattataagcaccctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGgtgt                                                                                              <  2:255036‑M2/182‑1 (MQ=255)
                                                                                                                                                      gaacagaaaaacggttattataatgaaactggctgtgagtcgctatagttaccaggtcagCGACTATTCCCGACTGAAATGGGATAGCTTCCGTAACCGAAAGCATAGGCTGCTGAGGTGACTTCGTAATGTTCATTCAGAATCTGACGCCTGACTTCATTGAGACGTAGACATTTCTGATACCGGAGAGGACTCATTACAGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTtcgtctg                                                                                     >  2:53924‑M2/61‑248 (MQ=255)
                                                                                                                                                          agataaacggttattataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAgt                                                                                 >  2:103078‑M2/57‑251 (MQ=255)
                                                                                                                                                                    ttattataatcaacctggctgggagtagatatagttgccaggtccgCGAATCTTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCTCCTGACTTCATTGAGACTTCGCCATTTCTTATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGTTGTTTTCGTATGCAGTTTCTAtgctc                                                                       >  1:165963‑M2/47‑249 (MQ=255)
                                                                                                                                                                         ataatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCgacagac                                                                  >  2:70628‑M2/42‑251 (MQ=255)
                                                                                                                                                                           aatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGacac                                                                <  1:216447‑M2/212‑1 (MQ=255)
                                                                                                                                                                           aatcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGacac                                                                >  2:131285‑M2/40‑251 (MQ=255)
                                                                                                                                                                            atcaacctggctgtgagtagctatagttgccaggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACt                                                               >  1:49182‑M2/39‑251 (MQ=255)
                                                                                                                                                                                              agctatagttgccaggtccgAGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATgggg                                              <  1:118855‑M2/230‑1 (MQ=255)
                                                                                                                                                                                                          caggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGCCACTGCGCTCTTAGCATGGGGCAGGGAATAACCa                                 <  1:364231‑M2/243‑1 (MQ=255)
                                                                                                                                                                                                          caggtccgCGAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAAcc                                  >  2:223853‑M2/9‑250 (MQ=255)
                                                                                                                                                                                                                   gAATATTCCCGACTGAAATGGGATAGCTTTCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACTCTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTtc                        <  1:289502/251‑1 (MQ=255)
                                                                                                                                                                                                                   gAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTtc                        >  1:241192/1‑251 (MQ=255)
                                                                                                                                                                                                                   gAATATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTtc                        >  1:29622/1‑251 (MQ=255)
                                                                                                                                                                                                                    aaTATTCCCGACTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTtct                       >  2:325814/1‑251 (MQ=255)
                                                                                                                                                                                                                              aCTGAAATGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTGAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAg             <  1:41984/251‑1 (MQ=255)
                                                                                                                                                                                                                                   aaTGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCa         >  1:222964/1‑250 (MQ=255)
                                                                                                                                                                                                                                    aTGGGATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAAt       >  1:348019/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ggATAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAATTAt    <  1:173880/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aTAGCTTCCTTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTGTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCCCAGGCCAATTATGa  <  2:331625/251‑1 (MQ=255)
                                                                                                                                                                                                                                         aTAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAg                                           >  2:280346/1‑210 (MQ=255)
                                                                                                                                                                                                                                         aTAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAg                                           <  1:280346/210‑1 (MQ=255)
                                                                                                                                                                                                                                          tAGCTTCCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCACTTATGa  >  1:61057/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                   
CAGAATATTATTGCGAACAAGGGAAGTGTCAGGAAAGGTCAAGAGCTACCAATGGTTAATATGGCTATAAATCATTAGCTATGGTTAATGTTTGCTTGCTGAAACTAACACCGAATTGCTCTAAACAGATAAGTGCAGGGTTTATTAAAGCTGTCCGACAGACTTTCTTAACCTGGTAATATCTCTCTTGGGTGACTCTCCAAACATCCGCGAATATTCCCGACTGAAATGGGATAGCTTTCGTAACCGACAGCATAGGCTGCTGTGGTGACATCGTAATGTTCATTCAGCATCTGTCGCCTGACTTCATTGAGACGTAGCCATTTCTGATACCGGAGTGGACTCATTCCCGCTTGATCCGTCCTTAAATAAGTAACAGCACCCTCTGGTGTTTTCGTATGCAGTGTCTATGCGACAGACACTGCGCTCTTAGCATGGGGCAGGGAATAACCAGGCTATTTCTCCATGCACAGGCCAATTATGA  >  NC_000913/1590190‑1590673

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