Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,325,644 C→T S33N (AGC→AAC)  yjdN ← PF06983 family protein YjdN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,6440CT81.0% 79.1 / 12.9 44S33N (AGC→AAC) yjdNPF06983 family protein YjdN
Reads supporting (aligned to +/- strand):  ref base C (5/3);  new base T (19/15);  total (26/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

CATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGC  >  NC_000913/4325400‑4325880
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
cATTTCGATTTTTCCGTTAGCGGCAAGCTTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACCTTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGAttt                                                                                                                                                                                                                                        >  1:81139/1‑251 (MQ=255)
    ccGATTTTTCCGTTAGCGGCAAGAGTGACAAACCAGAGTTTTCCGTCTTCGACCTGTGGCGAATCGAGCCCCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCAGCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCTATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATGTtata                                                                                                                                                                                                                                    <  2:81139/250‑1 (MQ=255)
         ttttCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGACGTTGATTTTATAGAGc                                                                                                                                                                                                                                >  1:28272/1‑250 (MQ=255)
         ttttCCGTTAGCGGAAAGATTGTCAAACCAGCGTTTTCCTTCTTCGCCCTGTTGCGAATCGCACACCAGCGTAAAGCCGGAGGATCTGGCTTTTCCTAACGGCATGGCATCGATCATCATGATGT‑CGCTTCCTGCAATGCGAACATTGGCGAGGGCGAGGGCGGTATCGGGGACTTGCATTCCTGAAGGGCAGTTCACGACGCTGACCTGAGCTGCTTTTGGCATTTCGCCGAAGGTGGTTATATagggc                                                                                                                                                                                                                                >  2:221749/1‑247 (MQ=255)
               gTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTACTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCcg                                                                                                                                                                                                                         >  1:334992/1‑251 (MQ=255)
                       cAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAac                                                                                                                                                                                                                 >  1:120591/1‑251 (MQ=255)
                           aTTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGACATTGCATTCCGGAAGGGCAGTTCTCGGCGCGGGCCTGCGCTGCTTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCcgggatcggtag                                                                                                                                                                                                               >  2:62985/1‑239 (MQ=255)
                           aTTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCcg                                                                                                                                                                                                                         <  1:62985/239‑1 (MQ=255)
                                cAAACCAGCTTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTg                                                                                                                                                                                                         <  1:258806/250‑1 (MQ=255)
                                  aaCCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGat                                                                                                                                                                                                       >  1:306457/1‑250 (MQ=255)
                                    ccatcGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATTGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTTCATTCCGGAAGGTCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGGCGAAGTTGATTTTATAAAGCAGTTCCGCGCCTAACGTACGTTGataat                                                                                                                                                                                                    <  1:85673/247‑1 (MQ=255)
                                             tCCTTCTTCGACCTGTTGCGAATCGAGACCCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATATAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATc                                                                                                                                                                                            <  2:288837/250‑1 (MQ=255)
                                                      gacatGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGCCGGCATGGCATAGCTCATCATGATGT‑CGCTTCCGGCAATGAGCACATTGGCGTGGGCGATGGCGGTCTCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGCTTTTGGCATTTCGCCGAAGTTGATTTTATAGAACAGTTCCGCGCCCAACGTACTGTGATAATCGGCAATCGCGTCGGacc                                                                                                                                                                                  >  2:182723/5‑249 (MQ=255)
                                                      gACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAAc                                                                                                                                                                                  >  1:108964/1‑251 (MQ=255)
                                                      gACCTGATGCGATTCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGATTACGGCAATGCGCACATTGGCGTGGGCGATGGCGGAATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGCTTTTGGCATTTCGCCGAAGATGAATTTATAGAGCAGTTCCGCGCCCACCGTACGTTGATAATAGGCAATCGCGTAGGAAc                                                                                                                                                                                  >  1:57683/1‑251 (MQ=255)
                                                                 aaTCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGc                                                                                                                                                                        >  1:279685/1‑250 (MQ=255)
                                                                        cACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGGAGGAACAGTTACCAGCAAAagggg                                                                                                                                                                >  2:313622/1‑248 (MQ=255)
                                                                             gCGTAAAGCCGGAGTAGCGTGCTTTTCCGGACGGAGTGGCATCGCTCAGCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGTAAATTGCATTCCTGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAGGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAg                                                                                                                                                           <  2:376157/251‑1 (MQ=255)
                                                                                tAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAggg                                                                                                                                                         >  1:257356/1‑250 (MQ=255)
                                                                                             gCTGGCTTTTCCTTACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACAGTGGTGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGTATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGGTGATAATAGGCAATCGCGTCGGAACAGTTACCGTCAAAAGAGAGGTAGGGACTTAACGGCatga                                                                                                                                           <  2:378627/251‑1 (MQ=255)
                                                                                                       cctggcgtgatggcatCGCTCAGCATGA‑GTGCGCGTCCTGCAAGGCGCACATTGGCGTGGGCGATGTCGGTATCGGGAAATTGCATTCCGTAAGTGCATTTCTTGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGTAACATTTACCGGCAAAAGAGAGGTAGGGACTTAACGTCATGATGGTGACCt                                                                                                                                  <  2:82909/241‑1 (MQ=255)
                                                                                                               catggcatCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGTGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGTTAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAaga                                                                                                                         <  1:317142/251‑1 (MQ=255)
                                                                                                                    catCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGc                                                                                                                    <  1:352864/251‑1 (MQ=255)
                                                                                                                       cGCTCATCATGATGG‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAg                                                                                                                 <  1:188527/251‑1 (MQ=255)
                                                                                                                              catGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAgtg                                                                                                          <  1:234060/251‑1 (MQ=255)
                                                                                                                                              ggCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGACGTTGATTTCATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGCACAGTTACCGGCAAAAGAGAGGGAGGGACTTAACGGCATGATGGTGACCTCCGTTAAGAGAAGCCAGGTAAGTGTAGTTCGGAATTtat                                                                                           >  2:116092/1‑251 (MQ=255)
                                                                                                                                              ggCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCatg                                                                                                                                            >  1:153853/1‑202 (MQ=255)
                                                                                                                                                cAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTGCGCCGAAGTTGATTTTATAGAGCAGTTCCTCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATAc                                                                                         <  1:182723/251‑1 (MQ=255)
                                                                                                                                                 aaTGCGCACATTGGCGTGGGCGATGGCGGGATCGGGAAATTGCATTCCGGAAGGGCAGTTCGCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACCGTTACCGGCAAAGGAGAGGTAGGGACTTAACGGCATGATGGGGTCCTCAGTTACGAGAAGCCAGTTAAGAGACGCTCGAAATTTATAc                                                                                         >  2:120988/1‑250 (MQ=255)
                                                                                                                                                  aTGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATTGTGACCTCAGTTAAGAGAAGCCATTTAAGTGTAGTTCGAAATGTATACAg                                                                                       <  1:331507/251‑1 (MQ=255)
                                                                                                                                                     cgcACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGCCTTAACGGCATGATGGGGACCTCCGTTAAGAGAAGCCAGTTAAGTGTAGTTCGGAATTTATACAGAt                                                                                     >  2:332028/1‑250 (MQ=255)
                                                                                                                                                      gcACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAA‑GAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATg                                                                                    >  1:80304/1‑249 (MQ=255)
                                                                                                                                                            tGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGGGGCCTCAGTTAAGAGAAACCAGTTAAGTGTTGTTCGTAATTTATACAGATgcaaggcg                                                                             >  2:143325/1‑244 (MQ=255)
                                                                                                                                                                        tGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCCGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACATATGAGAGGCAGCCTGATAAGAc                                                                 >  2:172069/1‑251 (MQ=255)
                                                                                                                                                                         ggCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCTAGTTAAGTGTAGTTCGAAATGTATACAGATGAGAGACGGCCTGATAAGAc                                                                 <  1:306791/250‑1 (MQ=255)
                                                                                                                                                                                tCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGGGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCATGATAAGAagcgcgcgc                                                         >  2:341718/1‑248 (MQ=255)
                                                                                                                                                                                     aaaTTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGc                                                    >  1:280314/1‑251 (MQ=255)
                                                                                                                                                                                               ccGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAg                                           <  1:230297/250‑1 (MQ=255)
                                                                                                                                                                                               ccGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTAATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATACTAGGCAATCGCGTCGGCACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGt                                          >  1:20202/1‑251 (MQ=255)
                                                                                                                                                                                                                         cTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTa                 >  1:117335/1‑250 (MQ=255)
                                                                                                                                                                                                                           gcgcTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcggcg                   <  2:33539/246‑1 (MQ=255)
                                                                                                                                                                                                                           gcgcTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATgcggcg                   >  1:33539/1‑246 (MQ=255)
                                                                                                                                                                                                                           gcgcTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAAc              <  1:172069/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gCATTTCGCCGAAGCTGATCTTATAGAGCAGTTCCGCGACCAACGTACGATGATAATAGGCAATCGCGTCGGAACAGCTACCGGCAACAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGACAGCGACGCATCAGGCAGTCGGACCTATTGCCGGTCGCGGCGTAAAAGCCttttacggc  >  2:155229/1‑244 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                            
CATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGT‑CGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGC  >  NC_000913/4325400‑4325880

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