Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,842,793 T→G 24.6% G141G (GGT→GGG gdhA → glutamate dehydrogenase, NADP‑specific

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,842,7930TG24.6% 117.8 / 2.7 70G141G (GGT→GGGgdhAglutamate dehydrogenase, NADP‑specific
Reads supporting (aligned to +/- strand):  ref base T (21/28);  new base G (18/3);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.33e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCGCGTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGG  >  NC_000913/1842558‑1843011
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                            
ttCGCGTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGTTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGt                                                                                                                                                                                                                              <  1:239747/236‑1 (MQ=255)
    cgTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGCGGTAAAGGCGGCAGCGATTTCGATCCGCACGGAAAAACGGAAGGGGAAGGGa                                                                                                                                                                                                                       >  1:1072095/1‑239 (MQ=255)
    cgTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAAAAAACCTTCAAAAATGCCCCTACTCCTCTGCCGATGGGCGGTGGTCAAGGCGGCACCGATTTCGATCCGACAGGAAAAAGCGAAGGTGagggga                                                                                                                                                                                                                       >  1:1043013/1‑234 (MQ=255)
    cgTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAAGCGGTATGCGCTTCCATCCGTCCGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTAATCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGACTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGTGa                                                                                                                                                                                                                       >  1:113858/1‑239 (MQ=255)
       ggttggggTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATg                                                                                                                                                                                                                     <  1:876643/233‑1 (MQ=255)
                  atgatCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGGGGCAGCGATTTCGATCCGAAAGGAAAAAACGAAGGTGAAGAGATGCGGTTTTGCCAg                                                                                                                                                                                                         >  1:779155/1‑239 (MQ=255)
                  atgatCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCCATGGGCGGGGGGAAAGGGGGCAACGATTTCTGTCCCAAAGGAAAAAGGGAGGGGGAAGGGGTGGGGTTTTGccca                                                                                                                                                                                                         >  2:600462/1‑237 (MQ=255)
                   tgatCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCCGTTAACCTTTCCATTCTCAAATTTCTCGGCTTTGAAACAACCTTCCAAAATGCCCTGACCACTCTACCCATGGGCGGTGGGCAAGGCGGCCGCGATCTCTAGCCGCAAGAAAAAAGCGAAGGGGaacagtgggggttttgacacg                                                                                                                                                                                                        >  2:317691/1‑220 (MQ=255)
                   tgatCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGACCGTAAAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCCAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGGAGTGATGCGGTTTTGCCAgg                                                                                                                                                                                                        >  1:860743/1‑239 (MQ=255)
                     actccacaccactaaaGGTAAACCGGGGCTGGCGTGTGCCGTTCAGCTCTGCCATCGGCCCGTACCAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCCCAAATTCCTCGGCTTTGAACAAACTTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGgcg                                                                                                                                                                                                      <  2:959787/224‑1 (MQ=255)
                       cGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAACGGCGGCAGCGATTTCGATCCGAAAGGCAAAAGCGAAGGGAAAGTGTTGCGGTTTTGCCAGgcgag                                                                                                                                                                                                    >  1:406349/1‑237 (MQ=255)
                       cGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGGAAGGGAAAAACGAAGGGGAAGGGGTGCGTTTTTGCCAggg                                                                                                                                                                                                       >  1:467454/1‑235 (MQ=255)
                       cGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCACCGATTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGGGATGCGGTTTTTCCAGgc                                                                                                                                                                                                       >  1:181974/1‑236 (MQ=255)
                       cGCAACCAGATAACGGGCAACCGTGCATGGGTGGTGCATTTCCGCCTTGCCAGCGGCCCGTACAAAGGCGGTTTGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCt                                                                                                                                                                                                   <  2:329027/239‑1 (MQ=255)
                         cAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGGAAGGAAAAAGCGAAGGTGAAGTGATGCGGTTTTGCCAGGCGCtga                                                                                                                                                                                                 >  1:3138/1‑239 (MQ=255)
                          ccccagctgcctgtaacccGTGCATGGCGTGTGCAGTTCCGCTCTGCCATCGGCCCGTACAAAGCCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGCAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCtga                                                                                                                                                                                                 <  2:406349/221‑1 (MQ=255)
                          aaCCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCtga                                                                                                                                                                                                 <  1:1077401/238‑1 (MQ=255)
                                 tACAGGTCAACCGTGCAGGGCTTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTg                                                                                                                                                                                           <  1:621734/237‑1 (MQ=255)
                                    aGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGGGATGGGTTTTTGGCAGGCGCGGCTGACCGACCtt                                                                                                                                                                                     >  1:166026/1‑238 (MQ=255)
                                      gTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCCAATTCCTCGGCTTTTAACCACCCCTCAAAAATGCCCTTACTACCCTGCCGATGGGGGGTGGTAAAGGCAGCAGCGGTTTCGATCCCAAAGGAAAACGGCCGGGTCACGTGAGTCGTTTTTCCCAGgccgcgtagaacaaacagt                                                                                                                                                                                      >  2:621734/1‑221 (MQ=255)
                                      gTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAACGAAGGGGAAGTGATGCGTTTTTGCCAAGCGAGTATGACTGAAATgga                                                                                                                                                                                   >  1:911469/1‑237 (MQ=255)
                                      gTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAAGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAACGTGAAAGGATGGGTTTTTGCAAGGCGTTGATGACTGAACTga                                                                                                                                                                                    >  1:437576/1‑237 (MQ=255)
                                      gTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGGGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAACGAAGGGGAAGGGAGGCGTTTTTGCAAgggcctgtggaacgaacttg                                                                                                                                                                                      >  1:959972/1‑220 (MQ=255)
                                      gTCAACCGTGCATGGCGCGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGGATGCGCTTCCATCCGGCAGTTAACCTTTCCATTCTCAAATTCCTCGGCCTTGAACAAACCTTCAAAAATGACCCGACTAACCTGCCGCTGGGCGGGGGTAAAGGGCGCAGCGATTTCGGTCCGCACGGAAACACCGAAGGGGAAGTAGCTCGTTTTTGTCAGGAGCTGAGGAcgtaccctg                                                                                                                                                                                     >  2:85577/1‑230 (MQ=255)
                                          aCCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGGAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGGGATGCGTGTTTGCCAGGGGCGGGTGACTTAACGGGATCg                                                                                                                                                                                >  1:181333/1‑238 (MQ=255)
                                              tGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTACATTTTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGGAAAAGCGCAGGTGGAGGGATGCGGTTTTGCCAGGCGCTGATGACTGAACTGTAACGCCAcc                                                                                                                                                                           >  1:255870/1‑239 (MQ=255)
                                                     cGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGGGATGCGTTTTTTCCAGGCGCTGAGGACTGAACCGGATGCCCACCTggggcg                                                                                                                                                                    >  1:1061745/1‑236 (MQ=255)
                                                        gtgCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGTGATGGGTTTTTGCCAGGCGCCGATGACTGAACTGTATCGCACCCCGGGCGCGGa                                                                                                                                                                 >  1:944867/1‑239 (MQ=255)
                                                         tgCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGGGATGCGTTTTTGCCAGGCGCTGGTGACTGAACTGTATCGCCACCTGGGCGCGGAt                                                                                                                                                                >  1:973077/1‑239 (MQ=255)
                                                              ttCAGCTCTGCCATCGGCCCGTACAAAGGCGGTCTGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAGATTCCTCGGCTTTGAACAAACCTTCAAAAAGTCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGACTGGATTTCGATCCGAATGGAAAAACCGGATGGTAAGTGAAGCGGTTTTGCCCGGGTCTGTGGACTGGACGGCAACACCACCGGGgccgggacacgac                                                                                                                                                           >  2:431955/1‑228 (MQ=255)
                                                                   ctctGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGGGGTAAAGCCGGCAGCGACTTCGATCCGAAAGGAAAAAGCGCAGGTGAAGGGAGGCGTTTTTGTACTGGGCTGCTGACGAACATGTAAGGCCACCTGGGCGCGAAGACAGACGTTc                                                                                                                                                      >  2:56743/1‑239 (MQ=255)
                                                                      tGCCCTCGGCCCGTACAAAGGCGGTTTGCGCTTCCACCGGTACGTTAACCTTTCCATTCCCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCgg                                                                                                                                                   <  2:982006/239‑1 (MQ=255)
                                                                        ccATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGGGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCa                                                                                                                                                 >  1:257430/1‑239 (MQ=255)
                                                                        ccATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAAACTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGCGATGCGGTTTTGCCAGGCGCTGATGACTGAACTGATTCGCACCCTGGGCCGGGATACCGACGGTCCGGcc                                                                                                                                                 >  1:53692/1‑238 (MQ=255)
                                                                        ccATCGGACCGTACAAAGGCGGGTTGCGCATCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGTCTTTGAACAAACCTTCAAAAATGCCCTGACTAATCAGCCGATGGGCGGAGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGTGATGCGTTTTTTCCGGGCGACAGAGAATTAACGGTATCGACACCTggggggggtaacgcaggtaaccagg                                                                                                                                                 >  2:338498/1‑217 (MQ=255)
                                                                           tCGGCCCTTACAAAGGCTGTATGCGCTTCCATCCGTCAGTTATCCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGt                                                                                                                                              <  2:467454/239‑1 (MQ=255)
                                                                                ctcgacaaggggcgttttggcTTCCAACAGTCCATTAACCTTTCCATGTTCAAATTCCCGGGCTTTCGAAAAACCCTCAAAAATGCCCTGACTACGCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGatat                                                                                                                                         <  2:386490/220‑1 (MQ=255)
                                                                                cccgaCAATGGCGTTTTGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCACGGCTTTGAACAATCCTGCAAAAATGCCCTGACCACCCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGatat                                                                                                                                         <  2:85741/235‑1 (MQ=255)
                                                                                 cGTACAAGGGGGGTATGCGCTTCCATCCGTCAGTTACCCTTTCCATTCTAAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGatat                                                                                                                                         <  1:1067332/238‑1 (MQ=255)
                                                                                           cGGTATGCTCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTtgg                                                                                                                                <  2:966172/237‑1 (MQ=255)
                                                                                                 gcccTTCCATCCGGCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTg                                                                                                                        <  1:259426/236‑1 (MQ=255)
                                                                                                    cTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAGAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGGGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCCCGGTTACCACCGTTCCGGCAGGGGATATCGGGGTTGGTGGTCGTGAAg                                                                                                                     >  1:500249/1‑239 (MQ=255)
                                                                                                     ttCCATCCGTAAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGt                                                                                                                    <  1:861920/239‑1 (MQ=255)
                                                                                                       ccATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGGGGTAACGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGGGATGCGTTTTTGCCAGGCGCTGATGACGGAACTGTATCGCACCCTGGGCGCGGATACCGCAGTTCGGGCAGGGGCTATCgggggttgggggggggaaggcg                                                                                                                  >  2:861920/1‑221 (MQ=255)
                                                                                                                   tttccTTTCCATTCTAAAATTTGTCGGCTTTGACCAAACCTTAAAAAAGGCCCTGTCTACTCGGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTTTCGCCACCTGGGCGCGGATCCCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCg                                                                                                       <  2:467483/235‑1 (MQ=255)
                                                                                                                               tctcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGaaaa                                                                                          <  1:285658/239‑1 (MQ=255)
                                                                                                                                ctcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGGGGGCGTGAAGTAGGCTTTATGGCgggggggttggaaaa                                                                                         >  1:646622/1‑228 (MQ=255)
                                                                                                                                 tcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGACGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGGATCGCCACCTGGGCGCGGATACCGCCGTTCCGGCAGGTGATATCGGGGTTGGGGGGCGGGAAGTCGGCTTTTTGGCgggggggggaagaaaa                                                                                        >  2:642919/1‑227 (MQ=255)
                                                                                                                                 tcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGTTACCGACGTTCCGGCAGGTGATATCGGGGCTGGAGGTCGCGAAGTCGGCTTTATGGCGGGgaggtcggaaaaa                                                                                        >  1:1075669/1‑228 (MQ=255)
                                                                                                                                    tattcctcccgtttgtaaccaacTTCAAAACTGCGCTGTCTCCTCTGCCGAGTGGCGTGTGTAAAGGCGGCAGCTATTTCGTTCCGAAAGGAACACGCGAAGTGGCAGTGTTGCGTTTTTGCCAGGCTCTTATGACTGAAGTGTATCGCCACCTGGGCGTGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAACTCTGCTTTATGGCGGGGATGATGAAAAAGctc                                                                                     <  2:297510/217‑1 (MQ=255)
                                                                                                                                    aaTTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGctc                                                                                     <  1:388481/239‑1 (MQ=255)
                                                                                                                                    aaTTCCCGGGCTTTGAACAAACCCTCAAAAATCCCCCGACTACCAGGGCCATGGGCGGTGGTAAAGGCTGCAGCGATGTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGctc                                                                                     <  2:911469/239‑1 (MQ=255)
                                                                                                                                           cGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGt                                                                                                                    <  2:860116/201‑1 (MQ=255)
                                                                                                                                           cGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGt                                                                                                                    >  1:860116/1‑201 (MQ=255)
                                                                                                                                                 tGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGGGATATCGGGGTTGGTGGACGGGAAGCCGGCTTTATGGCGGGGGGGATGacaaaaaagtcccacaacacc                                                                         >  1:704728/1‑222 (MQ=255)
                                                                                                                                                    aaaacccTCCAAAAATCCCTTGACAACTCTGCCGAGGGGGGTTGGTAAGGGGGGCAGCGTTTTCGATCCAAAAGGAAAAGGCGAGGTTGAAGTGGTGCTTTTTTGCCGGTCGCTGATGACTGAACTGTATCGCCCCCGGGGCGCGGATCCCGACGTCCCGGCAGGTGATATCGGGGTTGTTGTTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCACCAATCCCGCCt                                                                     <  2:113858/234‑1 (MQ=255)
                                                                                                                                                        aCCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCg                                                                                                                                       <  2:328754/169‑1 (MQ=255)
                                                                                                                                                        aCCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCg                                                                                                                                       >  1:328754/1‑169 (MQ=255)
                                                                                                                                                        aCCTTAAAAAAGGCCCTGACTACTTGGCCGATGGGGGGTGGTAAAGGCGGCAGCGATTCCGATCCGAAAGGAAAAAGGGAAGGTAAAGTGATGCTTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGt                                                                 <  1:578378/239‑1 (MQ=255)
                                                                                                                                                             cAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCTGAAAGGAAAAAGCGAAGGTGAAGTGATGCATTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGGGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTagg                                                                                                          >  1:350735/1‑191 (MQ=255)
                                                                                                                                                             cAAAAATCCCCTGCTTCCTGTGCGGGGGGGCGGGTGAAAAGGCGGCGGCGATTTCGATCCGAAGGGAAAAAGCGAAGGGGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGATCTGTACCGCCACGGGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGGCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCa                                                            <  2:854185/239‑1 (MQ=255)
                                                                                                                                                                     cccttgatctcttggcGAGGGGGGGTGGTAAAGGCGGCAGCGTTTTGGATCCGAAAGGAAAAAGCGAGGTTGAAGTGTGGCTTTTTGGCCAGGGGCTGATGACTGAACGTTATGGCCACCTGGGGGCGGATACCGCGTTTCCGGCAGGTGATATCGGGGTTGGTGGTCGGGAAGCCGGCTTTATGGCGGGGATGATGAAAAGGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTaa                                                     <  2:3138/223‑1 (MQ=255)
                                                                                                                                                                     cccggaattctcgGCCGATGGGCGGGGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAg                                                    <  1:613879/230‑1 (MQ=255)
                                                                                                                                                                      ccggACTTCCCTGGCGAGGGGCGGGGGTAAGGGGGGCAGGGATTTGGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGTTGATGACTGAACTGTATGGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGGGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAg                                                    <  1:85577/235‑1 (MQ=255)
                                                                                                                                                                              ctctGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAGAAAAACCTCCCACAACAACACCTCCCGCTTCCCCGGCGAGGGGCTTTCa                                           >  1:643894/1‑239 (MQ=255)
                                                                                                                                                                                      ggggggggTGGTAAAGGGTGCGGGTTTTCCGAGCCGAAAGGAAAAAGCGAAGGTGAAGTTATGCGTTTTTGGCAGGCCGTGATGACTGAATTTTACCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGAGGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTCCACCGGTAAGGGCCTTTCATTTggcgg                                   <  2:181333/237‑1 (MQ=255)
                                                                                                                                                                                       tgggggggtttAAGGGGGGGAGCGTTGTCGACCGGAGGGGAAAAAGAGAAGGGTGAGTGAGGTGTTTTTGCGGGGCGTCGATGACTAAATTGTCTCGCCCCCTGCGCGCGTATCCCGGAGTGCCGGGAGTTGATTCCGGGGGTGTTGTTCGTGAAGCCGGCTTTATGGCGGGGATGATGAAAAAGCTCCCCACCAATGCCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTggcgg                                   <  2:591259/228‑1 (MQ=255)
                                                                                                                                                                                                          cAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGaaa                                                                                           >  1:775231/1‑162 (MQ=255)
                                                                                                                                                                                                          cAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGaaa                                                                                           <  2:775231/162‑1 (MQ=255)
                                                                                                                                                                                                              gATTCCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCCCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTa            <  1:600462/238‑1 (MQ=255)
                                                                                                                                                                                                                          aaaGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAAGATTCCCACAATAACGCCCGCGGCTTCCCCCGTGAGGGCCTTTTCTTTTTCGGGAGTCTTATTTTCCCCGCAGCCACCCGccccgg  >  1:465290/1‑231 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                            
TTCGCGTGGTATGGGTTGATGATCGCAACCAGATACAGGTCAACCGTGCATGGCGTGTGCAGTTCAGCTCTGCCATCGGCCCGTACAAAGGCGGTATGCGCTTCCATCCGTCAGTTAACCTTTCCATTCTCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTGCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCCAACAATACCGCCTGCGTCTTCACCGGTAAGGGCCTTTCATTTGGCGGCAGTCTTATTCGCCCGGAAGCTACCGGCTACGG  >  NC_000913/1842558‑1843011

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