Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,489,016 C→T E350K (GAA→AAA)  yjgR ← DUF853 domain‑containing protein YjgR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,489,0160CT97.0% 205.1 / ‑2.6 67E350K (GAA→AAA) yjgRDUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (30/34);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

GCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGCATCGTTAAACAGCAGATGTGC  >  NC_000913/4488770‑4489262
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
gctgctCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGt                                                                                                                                                                                                                                                     <  1:3053/250‑1 (MQ=255)
          aCTGGCCTGAAAGCCTTTTTGCAGCATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATCCACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGAACGCTCCCCCACAGAAGGGCTTCCTTTCGCCTCCAGCAACGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATAGCCTTTTTGGTATCAAACGCCg                                                                                                                                                                                                                                          >  1:27335/1‑251 (MQ=255)
           cTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGAATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCgg                                                                                                                                                                                                                                         >  1:622017/1‑251 (MQ=255)
              cccggAAACCCTTTTTGCAACGTCTTCTAGGCGGATTCCCGGACCATCACACCCTCATTGTTGCCATACCCCGGAGAGTGATTAATCAAGCCAATACGCTCATCCTCCGTCACCGGCCCCATCCGCGAACAATGCGCGATCACCATCGCACGCTCCACCACAGAACGGCTTCCTTTCGCACCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGAtt                                                                                                                                                                                                                                      <  2:565003/247‑1 (MQ=255)
                               cAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACAACGGAGAGAGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGAGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGGAAAGAGATCAACGACTCGCCGGTGCCCAGTTCCTGAATCGACTTTTTGGTATCAAACGCCGGAGTAGCCCGCATGGTTTgcg                                                                                                                                                                                                                     >  1:180808/1‑251 (MQ=255)
                                    tctcATAGGCGGATTCCCGGTACACCTCATCCTCATATTTGCCATACACCGGAGCGTGATTAACCAAGCCATTACGCTCATCTTACGTCACCGGCCCCATCCGCGCACCAAGCGCGATCCCCATCGCACGCCCCACCACAGAAGGGCATCCTTTCGCCTCCTGAAACGAGCTCAACGCCTCGCCGGCGCCCAGTTCCTGAATCGACTTTTTGGTATCAAACGCCGGATTGGCCCGCacgggatgccccact                                                                                                                                                                                                                >  2:275769/1‑237 (MQ=255)
                                         gagccggctcccaggcccccctcatcctcatATTTCCCATACACCGGATAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCTATCAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCTCATCCAGAAAAGAGATCAACGCCTCGCCTGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCTCATGGTTTGCGCCGCTGCCtt                                                                                                                                                                                                           <  2:156755/232‑1 (MQ=255)
                                           ggCGGATTCCCTTCCCACCTCAGCCTCATATTTTCCATACGCCTGTGTACTATTAATCAAGCCACTACGCTCATCTCCCGTCCCCGGCCCCATCCGCGAACAAGTCGCGATCACCATCGCACGCTCCACCACAGAATGGCTTCCTTTCGCATCCAGAAAATAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCTCTGCCTTCa                                                                                                                                                                                                         <  2:1134/251‑1 (MQ=255)
                                           ggCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTc                                                                                                                                                                                                          >  1:315775/1‑250 (MQ=255)
                                              ggATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGccc                                                                                                                                                                                                                                 >  2:203200/1‑224 (MQ=255)
                                              ggATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGccc                                                                                                                                                                                                                                 <  1:203200/224‑1 (MQ=255)
                                                    ccGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCACCGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCCCTGCTTTCt                                                                                                                                                                                                >  2:470108/1‑251 (MQ=255)
                                                     cgggccACCTTATCCTCACATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTGCACTGCTTTCTg                                                                                                                                                                                               <  2:339211/247‑1 (MQ=255)
                                                         ccACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCATTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCt                                                                                                                                                                                           >  1:384464/1‑251 (MQ=255)
                                                          cACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCtt                                                                                                                                                                                          <  1:85185/251‑1 (MQ=255)
                                                             ctcatcctcatATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTg                                                                                                                                                                                       >  1:566537/1‑251 (MQ=255)
                                                              tcatcctcagATTCGCCCTCCACCGTAAAGTGATTAATCATGCCATTGCGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCGCAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGTTGCCCCGTTCCTGAATCGCCGTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTTCGCCGCTGCCTTCACGGCTTTCTGATCttgt                                                                                                                                                                                        <  2:661872/249‑3 (MQ=255)
                                                               cctcctcatATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTg                                                                                                                                                                                       <  1:194464/247‑1 (MQ=255)
                                                                           tGCCATACCCCGGACATTGATGACTCAAGCCCTTACGCTCATCTTCCGTCCCCGGCCCCATCCGCGAACAAGGCGCGACTCCCATCGCACGCGCCCCCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCTAACGCCGGATTGGGCCGCATGGTTTGCGCCGCGGCCTTCACTGCTTTCTGATCTTTTTGCTTAAAAGCCCg                                                                                                                                                                          <  2:119420/250‑1 (MQ=255)
                                                                           tGCCATACACCGGGGATTTGTTAAACAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTTCGCCGCt                                                                                                                                                                                                                <  2:244470/212‑1 (MQ=255)
                                                                           tGCCATACACCGGAGAGTGATTAATCACGCCATTAAGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGCAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCt                                                                                                                                                                                                                >  1:244470/1‑212 (MQ=255)
                                                                            gCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAAc                                                                                                                                                                                                                                              >  1:9104/1‑181 (MQ=255)
                                                                            cccATACACCGGAGAGCGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAAc                                                                                                                                                                                                                                              <  2:9104/180‑1 (MQ=255)
                                                                               aTACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAg                                                                                                                                                                     <  1:466600/251‑1 (MQ=255)
                                                                                 accccGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCa                                                                                                                                                                   <  1:249927/248‑1 (MQ=255)
                                                                                   acCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCAGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATg                                                                                                                                                                 >  1:559040/1‑251 (MQ=255)
                                                                                     cGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAACAGCCCGCAAAGCATGtt                                                                                                                                                               >  1:539971/1‑251 (MQ=255)
                                                                                            gTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAcgc                                                                                                                                                         >  1:32650/1‑251 (MQ=255)
                                                                                              gATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAcgcg                                                                                                                                                        >  2:224380/1‑250 (MQ=255)
                                                                                                  aaTCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTAc                                                                                                                                                   >  1:153221/1‑251 (MQ=255)
                                                                                                   aTCAAGCCATTACGCTCATCGTCCTTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCAGCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTAcc                                                                                                                                                  <  1:472301/251‑1 (MQ=255)
                                                                                                           tttccGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCACACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCAATGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGATCTGcc                                                                                                                                          >  2:362733/5‑251 (MQ=255)
                                                                                                           aTTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAcg                                                                                                                                                          <  2:144833/235‑1 (MQ=255)
                                                                                                           aTTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAAcg                                                                                                                                                          >  1:144833/1‑235 (MQ=255)
                                                                                                               cGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTg                                                                                                                                                                                       >  1:450133/1‑202 (MQ=255)
                                                                                                                      cTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGCGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTCGCGCCGCTGCCTTCACTGCTTTCTGCTCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCTCGCACAtt                                                                                                                                >  2:135176/1‑250 (MQ=255)
                                                                                                                                 ggCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGCAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGGGCTGCCCGAGAACATTCTCCGGAATATc                                                                                                                    >  2:272773/1‑251 (MQ=255)
                                                                                                                                       agccGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACgg                                                                                                              <  2:32650/249‑1 (MQ=255)
                                                                                                                                       aTCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCt                                                                                                                                                                                                            <  2:524756/157‑1 (MQ=255)
                                                                                                                                       aTCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCt                                                                                                                                                                                                            >  1:524756/1‑157 (MQ=255)
                                                                                                                                         ccGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGt                                                                                                            <  1:272773/251‑1 (MQ=255)
                                                                                                                                           gcgAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGttt                                                                                                          <  2:229490/251‑1 (MQ=255)
                                                                                                                                             gAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATATTTTGGCGTAAAAGCCCGCAAAGCATGTTGACCGCGATTACCGAGCTGCCCGAGCACATTCTCCGGACTATCAGCCGGGTTTTg                                                                                                        >  2:206517/1‑251 (MQ=255)
                                                                                                                                                 aaGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGa                                                                                                     >  1:236919/1‑250 (MQ=255)
                                                                                                                                                             cccATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGAc                                                                                         <  1:224380/250‑1 (MQ=255)
                                                                                                                                                                      cGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTGTGCGCCGCTGCCTTCACTGCTGTCTGAGCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATGACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGcc                                                                                <  1:615705/251‑1 (MQ=255)
                                                                                                                                                                        cTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATcc                                                                                                                            <  1:486195/205‑1 (MQ=255)
                                                                                                                                                                          ccaccaCAGAAGGGTGTCCTTTCGCATCCATAAAAGAGATCAACGCCTCGCCGGTGCCCATTTCCTGAATCGCCTTTTTGGTATCAAACTCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTCCACTGCTTGCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTGTTGCGAAACTAACCAGACGCCCTCGccgt                                                                              <  2:908/249‑3 (MQ=255)
                                                                                                                                                                               acaGAAGGGCTTCCTGTCGCAGTCAGAAAAGAGATCAACGCCGCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATAACAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCg                                                                       <  2:286185/251‑1 (MQ=255)
                                                                                                                                                                                    aGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATaa                                                                  <  1:54453/251‑1 (MQ=255)
                                                                                                                                                                                     gggCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCAGGTGCCAAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTAGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATa                                                                   >  1:244082/1‑249 (MQ=255)
                                                                                                                                                                                      ggCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAACGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGCATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGAGAATa                                                                   >  2:28681/1‑248 (MQ=255)
                                                                                                                                                                                                 gCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTg                                                     <  1:679435/251‑1 (MQ=255)
                                                                                                                                                                                                       aGAAAAGAGATCAACGCCTCGCCGGAGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAAt                                               >  1:219413/1‑251 (MQ=255)
                                                                                                                                                                                                           aaGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCtt                                            >  1:480581/1‑250 (MQ=255)
                                                                                                                                                                                                           aaGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTa                                           <  1:28681/251‑1 (MQ=255)
                                                                                                                                                                                                               gaTCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATcc                                                                                                                            >  1:76808/1‑166 (MQ=255)
                                                                                                                                                                                                               gaTCAACGCCTCGCCGGTGCCCAGTTCCGTAATCGCCTTTGTGGTATCAAACGCCGGAGTGTCCCGCACGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCGTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCca                                       <  2:206983/251‑1 (MQ=255)
                                                                                                                                                                                                               gaTCAACGCCTCGCCGGTGCCCAGGTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTGTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGATCACATTATcc                                                                                                                            <  2:76808/166‑1 (MQ=255)
                                                                                                                                                                                                                  cAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATATTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCcagca                                    <  1:206517/251‑1 (MQ=255)
                                                                                                                                                                                                                     cGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTa                                 <  1:470108/251‑1 (MQ=255)
                                                                                                                                                                                                                        cTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCAGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGATTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGGTTTGCGAAACGAACCAGACGCCCACGCATTTTGAGCGAATAAGCCTTATCACCTGCCCAATCTTATCCAGCAGCACCt                              >  2:333954/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ccAGTTCCTGAATCTGCGTTTTGACATCAACCGCCGGATTGGCCCGCATGGTGTGCGCCGCTGCCTTCACTGCTTTCGTATCGTTGGGCGGAAAAGCCCGCAAATCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGATACGATCCAGACGCCCACGCCGTTTGCGCGAATAAGCCTTATCCCCTGCTCAATCTTATCCAGCAGTCACTGCggtgtagcg                   <  2:39181/251‑8 (MQ=255)
                                                                                                                                                                                                                                   ccAGTTCCTGAATCGCCGTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCGGCTCAATCTTATCCAGCAGTACCTGCGGGGCATCg                   <  1:601789/251‑1 (MQ=255)
                                                                                                                                                                                                                                            gAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGCATCGTTAAAcag           <  1:625275/250‑1 (MQ=255)
                                                                                                                                                                                                                                                  gCCTTTTTGGTATCAAACGCCGGAGTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAAGCTTAGCCAGCAGTACCTGCGGGGGATCGTTAAACAGCAGAtg     <  1:179651/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     ttttTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAATACCTGCGGGGCATCGTTAAACAGCAGAtgggc  >  2:411788/1‑248 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
GCTGCTCGGTACTGGCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGCATCGTTAAACAGCAGATGTGC  >  NC_000913/4488770‑4489262

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: