Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,566,967 C→T E94K (GAA→AAA)  yjiN ← DUF445 domain‑containing protein YjiN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,566,9670CT98.6% 231.3 / NA 69E94K (GAA→AAA) yjiNDUF445 domain‑containing protein YjiN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (34/34);  minor base G (0/1);  total (34/35)
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.00e+00

CTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTCCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAG  >  NC_000913/4566723‑4567204
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
ctctCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGAATCAAGAAATTTTTTCTGCAc                                                                                                                                                                                                                                         >  1:322116/1‑251 (MQ=255)
  ctcCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACg                                                                                                                                                                                                                                        >  2:226732/1‑250 (MQ=255)
       gCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCACCCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTgg                                                                                                                                                                                                                                  >  2:463056/1‑251 (MQ=255)
         aTCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGcc                                                                                                                                                                                                                                >  1:97337/1‑251 (MQ=255)
              cGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCg                                                                                                                                                                                                                               >  1:417846/1‑247 (MQ=255)
                    ggTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTc                                                                                                                                                                                                                       <  2:119470/249‑1 (MQ=255)
                      tGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGcc                                                                                                                                                                                                                    >  1:659529/1‑250 (MQ=255)
                          ggAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATcc                                                                                                                                                                                                               <  1:243804/251‑1 (MQ=255)
                          ggAAAGATCGACCTTATCAATCGCCCGATGGAACGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCATCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCAGAGACTTTCGCCTATcc                                                                                                                                                                                                               <  2:23780/251‑1 (MQ=255)
                            aaaGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCgg                                                                                                                                                                                                             <  1:633497/251‑1 (MQ=255)
                              agagcgtccTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCAGCATCTGTAAGTGCGAGAACCCCGCTCATGGTCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGTATTGAGTATCAAGAAATTTTTTCTGCACGAACTTGCCGTGATTTTCGGCAATCCGGTc                                                                                                                                                                                                           <  2:51060/244‑1 (MQ=255)
                                           cAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCggg                                                                                                                                                                                              >  2:18311/1‑251 (MQ=255)
                                            aaTCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTCAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCggg                                                                                                                                                                                              >  1:92008/1‑250 (MQ=255)
                                            aaTCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGa                                                                                                                                                                                             <  1:196342/251‑1 (MQ=255)
                                               cGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTTCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATa                                                                                                                                                                                           <  1:375196/250‑1 (MQ=255)
                                                gCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGGTCGGGGCGTCGAATCAATGCCACCAGGGATTGAGGATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCGTTATTACGCGGGATAAt                                                                                                                                                                                         >  2:493268/1‑251 (MQ=255)
                                                   cGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGc                                                                                                                                                                                      >  1:636510/1‑251 (MQ=255)
                                                    gATGTACCGCACGCTTAAGCAGACGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCCACAGATGCTGACCTACGCTGCGCGCGTTTTCTGGCTGGGTGAACCAGTTACACATCAGCAACGGCGGGTCGTGGCGTCGAATGACTGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTTCACGAAATGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGcc                                                                                                                                                                                     >  2:419603/1‑251 (MQ=255)
                                                      gggACCGCGCGCGTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCACTCATGATCTGCGACAGATGCTGACCAACGCGGCGGGCGTTTTCTTGCTTGCTAAACCAGTTGCCAAGCAGCAACGCCGGTTCGGGGCGTCGGATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTTGTCGAGATTTTCTCCAATCCGGTCTTTATTACGCGGGATAATCTCCGt                                                                                                                                                                                   <  2:79798/250‑1 (MQ=255)
                                                           cgcgcgCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGc                                                                                                                                                                              <  2:659529/251‑1 (MQ=255)
                                                             cgcgcTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCCCCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCg                                                                                                                                                                             >  1:383476/1‑250 (MQ=255)
                                                                gcTTAAGCACGCGCTGAATAAGCGCATCATCGGTAAGTTAGATAAAACCGCTCATGATCTGCAACAGATGCTGACGAACGAAGCGGGCGGTTTTTTACTGGCTAAACCAGTTGCCAATCAGCAACACCGGTTCGTGTCGTTCAATCAATGCCACCAGGGGTTGAGAATCAAGAAATTTTTTCTGCACGAACTGGCCGTGGTTTTCTCCAATCCGGTCTTTCTTACGCGGGATAATTGCCgaattgcggga                                                                                                                                                                          >  2:314622/1‑240 (MQ=255)
                                                                 cTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGaa                                                                                                                                                                         >  2:547023/1‑250 (MQ=255)
                                                                 cTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGGCTTTATTACGCGGGATAATCGCAGTATGGCGAGaa                                                                                                                                                                         >  2:341000/1‑250 (MQ=255)
                                                                    aaGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGa                                                                                                                                                                     <  2:97337/251‑1 (MQ=255)
                                                                        aGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCgg                                                                                                                                                                 >  2:447812/1‑251 (MQ=255)
                                                                          gcgcTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGGCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGa                                                                                                                                                                >  1:374370/1‑250 (MQ=255)
                                                                            gcTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATc                                                                                                                                                             >  1:566995/1‑251 (MQ=255)
                                                                            gcTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATc                                                                                                                                                             <  1:390198/251‑1 (MQ=255)
                                                                                   aCGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCAcgc                                                                                                                                                      <  1:45040/251‑1 (MQ=255)
                                                                                   aCGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAACAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCAcgc                                                                                                                                                      >  1:4644/1‑251 (MQ=255)
                                                                                      cgcATCATCGGTAAGTTCGAGCAAACCGCTCATGACCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTCAACCAGTCGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGCAATGATCGGAATCGGCCCgcggc                                                                                                                                                   >  2:82145/1‑251 (MQ=255)
                                                                                           catcGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTg                                                                                                                                                                                                                                   <  2:511280/166‑1 (MQ=255)
                                                                                           catcGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTg                                                                                                                                                                                                                                   >  1:511280/1‑166 (MQ=255)
                                                                                           catcGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGaaa                                                                                                                                               <  2:407044/250‑1 (MQ=255)
                                                                                                    gTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATc                                                                                                                                                             >  1:217882/1‑227 (MQ=255)
                                                                                                    gTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATc                                                                                                                                                             <  2:217882/227‑1 (MQ=255)
                                                                                                           aaaaCCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCa                                                                                                                               >  1:188305/1‑250 (MQ=255)
                                                                                                                  cTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTc                                                                                                                       >  1:290389/1‑251 (MQ=255)
                                                                                                                                  gATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGAcc                                                                                                       <  2:559107/251‑1 (MQ=255)
                                                                                                                                    ggCTGAACAACGCGGCGGGCGTTTGCGGGCTGGCTAACACAGGTGCCAATCAGCAACGCCGGTGCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTGCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCAt                                                                                                     <  2:188351/250‑1 (MQ=255)
                                                                                                                                        gACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATcgc                                                                                                  <  2:87149/250‑1 (MQ=255)
                                                                                                                                                     ggcgTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCg                                                                                    >  2:382643/1‑251 (MQ=255)
                                                                                                                                                       cgTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGACCTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGcc                                                                                  <  1:358086/251‑1 (MQ=255)
                                                                                                                                                          tttCTGGCTGGCTAAACCAGTTGCCAATCAGCAGCGCCTGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCtt                                                                                <  2:430022/250‑1 (MQ=255)
                                                                                                                                                                 ctggctAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCg                                                                                    >  1:652966/1‑239 (MQ=255)
                                                                                                                                                                 ctggctAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCg                                                                                    <  2:652966/239‑1 (MQ=255)
                                                                                                                                                                  tggctAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCt                                                                       >  1:352262/1‑251 (MQ=255)
                                                                                                                                                                   ggctAAACCAGTTGCCAATCTGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACGGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAAGGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCt                                                                       <  1:602988/250‑1 (MQ=255)
                                                                                                                                                                    gctAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTc                                                                      >  2:177262/1‑250 (MQ=255)
                                                                                                                                                                      tAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGATTAAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAGCTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGACTCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCtnccc                                                                   >  2:607122/1‑247 (MQ=255)
                                                                                                                                                                      tAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCggg                                                                                                                                                                                              <  2:146214/128‑1 (MQ=255)
                                                                                                                                                                      tAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCggg                                                                                                                                                                                              >  1:146214/1‑128 (MQ=255)
                                                                                                                                                                       aaaCCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTGCTGCACGAACTGGCCGAGATTTGCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCAc                                                                    <  2:141819/249‑1 (MQ=255)
                                                                                                                                                                               ggCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCGCGAACTTGCCGAGATTTTCGCCAATCCGGTTTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAACCGGCTCGCGGCGAAACAGCGCCACCACCGCAACCCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCAAAAAAtt                                                           <  2:500250/249‑1 (MQ=255)
                                                                                                                                                                                      cAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTAGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAACATTGggcggc                                                    >  1:604003/1‑250 (MQ=255)
                                                                                                                                                                                               cGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTATGCACGAACTGGCCGAGATTTTCGCAAATCAGGTATTTATTAAGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGAGCCACCACCGCAAACCAGTCAGACAGCGCGCAGACCATCGACGCTTCAGCAATCGCCTTCACGCCGATCACACAACAATTGGGCGGCAAAAACAGCg                                          >  1:518684/1‑251 (MQ=255)
                                                                                                                                                                                                    cGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACg                                     <  1:29284/251‑1 (MQ=255)
                                                                                                                                                                                                    cGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACg                                     <  1:226732/251‑1 (MQ=255)
                                                                                                                                                                                                     gTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAGCAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAACAACAGCGTAACg                                     >  1:584729/1‑250 (MQ=255)
                                                                                                                                                                                                                 cAATGCCACCAGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTagcag                        <  1:18311/251‑1 (MQ=255)
                                                                                                                                                                                                                      ccaccaGGGATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGTGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAgcggcg                   <  1:253251/251‑1 (MQ=255)
                                                                                                                                                                                                                         ccaGGGATTGAGTATCAATAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACATCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATa                <  2:292071/251‑1 (MQ=255)
                                                                                                                                                                                                                          cggggATTTAGTATCAAGAAATTTTTTCCGCACAAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATaa               <  2:566995/249‑1 (MQ=255)
                                                                                                                                                                                                                           agtgATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAAGCGGCACGCGGCGAAACAGCGCCACCCCCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAg              <  2:544916/248‑1 (MQ=255)
                                                                                                                                                                                                                              gATTGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAg           <  2:4644/251‑1 (MQ=255)
                                                                                                                                                                                                                                 tGAGTATCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTa         >  1:144969/1‑250 (MQ=255)
                                                                                                                                                                                                                                      aTCAAGAAATTTTTTCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAagag    <  1:463056/250‑1 (MQ=255)
                                                                                                                                                                                                                                        cAAGAAATTTTTTCTGCACGAACGGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGAAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCTCCAGCGCGCCGACCATCGCCGCTGCAGCAATCGCCTTCACGCCGCTCACCCAAAAATGGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAagagag  <  1:314622/250‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                             
CTCTCCAGCATCAACGCACTGGTGCCGGAAAGATCGACCTTATCAATCGCCCGATGGACCGCGCGCTTAAGCAGGCGCTGAATACGCGCATCATCGGTAAGTTCGAGAAAACCGCTCATGATCTGCAACAGATGCTGACCAACGCGGCGGGCGTTTTCTGGCTGGCTAAACCAGTTGCCAATCAGCAACGCCGGTTCGTGGCGTCGAATCAATGCCACCAGGGATTGAGTATCAAGAAATTTTTCCTGCACGAACTGGCCGAGATTTTCGCCAATCCGGTCTTTATTACGCGGGATAATCGCCGTATGGCGAGAAATGATCGGAATCGGCACGCGGCGAAACAGCGCCACCACCGCAAACCAGTCCGCCAGCGCGCCGACCATCGCCGCTTCAGCAATCGCCTTCACGCCGCTCACCCAAAAATTGGGCGGCAAAAACAGCGTAACGACAAAGGTAGCAGCGGCGATAAGCAGTAAAGAGAG  >  NC_000913/4566723‑4567204

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