Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,791,436 G→C 13.0% A56P (GCC→CCC)  gabD → succinate‑semialdehyde dehydrogenase I, NADP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,791,4360GC13.0% 122.7 / 2.4 65A56P (GCC→CCC) gabDsuccinate‑semialdehyde dehydrogenase I, NADP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (26/28);  new base C (7/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 3.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGG  >  NC_000913/2791202‑2791653
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                         
gAGCGGCACGTAGTGTGGGTCCCTTACACGCCGGATTTAATCAATAACCTTTGAAAACAGGATGGAGCGATGAAACTTAACGACAGTAACTTACTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGAcgccgc                                                                                                                                                                                                                       <  1:882477/239‑1 (MQ=255)
                  aTGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGCAAATGGGCGCGGCTGAAACCCGCGCCGCTATCAACCCCGCCAACCGCGCCCTGcccc                                                                                                                                                                                                      >  1:858234/1‑237 (MQ=255)
                    gCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGCCGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGCAAATGGGCCCGGATGAAACCCGCCCCGCTATCGACCCCCCCAAccccccccccccccccc                                                                                                                                                                                                   >  2:379881/1‑230 (MQ=255)
                      cTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCCCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCCGGCCCTCCCCCCCTg                                                                                                                                                                                                  >  1:792071/1‑238 (MQ=255)
                       ttACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGAGCGCGGATGAAACCCGCCCCCCTATCGACGCCCCCAACCGCGCCCCGCCCCCCTGgc                                                                                                                                                                                                >  2:52129/1‑239 (MQ=255)
                       ttACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGAAGCGTGCCGAAAATGGCCGCGGATGAAACCCGCGCCGCTATCGACGGCGCCAACCCCGCCCCGCCCGCCGGgc                                                                                                                                                                                                >  2:780985/1‑239 (MQ=255)
                         acacGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGCAACTTATTCCGCCAGCAGGCGTTGATTAACGAGGAATGGCTGGACGCCAACAATGGTGAAGCCATCCACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCCCCGCTATCGACCCCGCCAACCTCGCCCCGCCCCCCTGgggc                                                                                                                                                                                              >  2:642216/1‑239 (MQ=255)
                              ccgcATTTAACCAATAACCTTTGAAAACGGGATGTAGCAAGGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCt                                                                                                                                                                                         <  1:642216/239‑1 (MQ=255)
                               cgcATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGTCGAACGGCGACCAGCTGGGAAGCGTGCCGAAAAAGGCCGCGGATGAAACCCGCCCCCCCACCGCCCCCGCCCACGCCGCCCCGcccccccgccccccccc                                                                                                                                                                                         >  2:637933/1‑227 (MQ=255)
                               cgcATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACCCCGCCAACCGCGCCCTGCCCTCCTGGCGCGCGCt                                                                                                                                                                                         >  1:267575/1‑238 (MQ=255)
                               cgcATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGACGCGGCTGAAACCCGCGCCGCTATCCACCCCCCCACCCcccccctgcccccctgccccccgcc                                                                                                                                                                                         >  2:87687/1‑229 (MQ=255)
                                gcATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCa                                                                                                                                                                                       <  1:235974/239‑1 (MQ=255)
                                             aaCCTTTGAAAACGGGTTTGGGCGATGAACATTAACGGACGTAACTTTTGCCGCCGGAGGGAGTTGATTAACGGGGAATGGCTGGACGCCAACAAGGGTGAAGCCATGGACGTCCCCAATCCGGCGAACGGCGCCAAGCTGGGTCGCGTCCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAAcg                                                                                                                                                                          <  2:376629/239‑1 (MQ=255)
                                                  ttaaaaaCAGGATGTAGGGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCGCCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGcca                                                                                                                                                                     <  1:204029/236‑1 (MQ=255)
                                                  ttGAAAACAGGATGTAGCGATGAAACTTAACGGCAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGcca                                                                                                                                                                     <  1:379881/239‑1 (MQ=255)
                                                  ttGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCCCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGCCGCGCGCCCACCGCCAAAGAACGCGcc                                                                                                                                                                      >  1:553666/1‑238 (MQ=255)
                                                   tGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCCCCGGGCGCGGGTTCACCGCCAAAGAACGCGccc                                                                                                                                                                     >  1:817230/1‑237 (MQ=255)
                                                               gTACCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCCCCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCaa                                                                                                                                                        <  1:466534/239‑1 (MQ=255)
                                                                   cGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAACATGGGCGCGGATGAAACCCGCCCCGCTATCGACGCCGCCCACCGGGCCCTGCCCGCCGGCACCGGCCCCACCGCCAAAGACCGCGcccccctcccccccacctgg                                                                                                                                                    >  2:72006/1‑224 (MQ=255)
                                                                    gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCCCCATTCTGCGCAACTGGt                                                                                                                                                   >  1:790975/1‑238 (MQ=255)
                                                                      tGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCCTTCcgcgacaccggtc                                                                                                                                                  >  1:347547/1‑229 (MQ=255)
                                                                         aactttaacgcggTAACTTATTCCGCCGGCGGGCGTTGATTAAGGGGGAAGGGCTGGACTCCAAAAAGGGGGAAGCCATTGACGTCACCAATCCGGCGAACGGCGACAAGCGGGGTAGCGTGCCGAAAATGGGCCCGGATGAAACCCGCGCCGTTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCaa                                                                                                                                               <  2:598428/226‑1 (MQ=255)
                                                                                cGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTtgataa                                                                                                                                       >  1:814532/1‑237 (MQ=255)
                                                                                   cAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCa                                                                                                                                                                                 <  1:31953/194‑1 (MQ=255)
                                                                                   aaGTAACTTATTCCGCCAGCGGGGTTGGATTAAGGGGGAAGGGCTGGAGCGCAAAAGTGGGGAGGCATTCGACGTCCCCAATGCGGCGAGCGGGGACAAGGGGGGTGGCGGGCCGAAAATGGGGGCGGATGAACCCCGGGCCGCTACCGACGCCGCCAACCGCTCCCCGCCCGCCTGGCGCGCGCTTACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGatgatg                                                                                                                                     <  2:40782/237‑1 (MQ=255)
                                                                                        aCTTATTCCGCCAGCGGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCAt                                                                                                                               <  1:72006/239‑1 (MQ=255)
                                                                                            aTTCCGCCAGAAGGCGTTGATAAACGGGAAATGGCGGGACGCCAAAAATTGGGAAGCCATCGACGGCACCAATCCGGCGAACGGCGACAAGCGGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAgg                                                                                                                           <  2:550581/239‑1 (MQ=255)
                                                                                              tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACCCCGCCAACCGCGCCCTACCCGCCTGGCCGCCGCCCACCCCCAAAGAACGCCCCACCACTCCCCGAAACCGGTTCACttacttgcaaaggaacaagaac                                                                                                                         >  2:282306/1‑219 (MQ=255)
                                                                                                       cggtgTTGATTAACGGGGAATGGCCGGCCGCCAACAAGGGGGAGGCCAGCGGCGTCACCAATCGGGCGAACGGCGACAAGCGGGGTGGCGTCCCGAAAATGGGCGCGGAGGAAACCCGCGCCGGTAGCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                <  2:675592/235‑1 (MQ=255)
                                                                                                       agggagggtttAACGGGGAATGGGTGGACGCCAACAATGGTGAACGCCTGGACGTAACCAACCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAAGGGGCCCGGAGGAACCCCGCGCCGCTACCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                <  2:384139/230‑1 (MQ=255)
                                                                                                          cGTTGATTAACGGGGAACGGCGGGCCGCCAACAATGGTGAAGCCATGGACACCCCCAATCCGGCAAACGGCGCCAAGCTGGGTAGCGTGCCGAAAATGGGCGCGAATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                <  2:267575/236‑1 (MQ=255)
                                                                                                               aTTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCAAAAGAACCCGCCACCATTCTCCGCAACTGGTTCAATTTGATGATGGAGCATCGGGACGATTTAGCGCGCCtgat                                                                                                        >  1:134727/1‑239 (MQ=255)
                                                                                                                ttAACGGGGAATGGCTGGACGCCAACAAGGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCtgat                                                                                                        <  1:87687/238‑1 (MQ=255)
                                                                                                                         aaTGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGAccggtc                                                                                                >  1:384139/1‑233 (MQ=255)
                                                                                                                         aaTGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGAGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCCCACCGCCAAACAACGCGCCACCATTCTGCGCACCTGGGTCAATTTGATGACGGAGCATCAGGACGATTTTCACCGCCCGATGACCCTcta                                                                                               >  2:77340/1‑236 (MQ=255)
                                                                                                                          ctggctggcCCCCAAAATTGGGGAAGCCATCGCGGCCCCCAATCCGGTGAACGGCGACAAGCGGGGTAGCGTGCCGAAAAGGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCTCCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGaa                                                                                              <  2:792071/237‑1 (MQ=255)
                                                                                                                          aTGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGaa                                                                                              <  1:780985/238‑1 (MQ=255)
                                                                                                                          aTGGCGGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAAGGGGCGCGGATGAACCCCGCGCCGCTATCGACGCCGCCAACCGCCCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGaa                                                                                              <  1:824870/238‑1 (MQ=255)
                                                                                                                          aTGGCGGGACCCAAAAAAGGGGGAAGCCCCCGCAGTCAGCAAGCCGGCAGACGGCGACAAGGGGGGGAGCGGGCCGAAACTGGGCGCGGATGAAACCCGCGCCGCTAGCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCAGCATGCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGaa                                                                                              <  2:858234/238‑1 (MQ=255)
                                                                                                                             gCTGGAGCCAAACAAGGTAGAAGCCAGCGACGGCACTACGCCCTGGAACCGCGACAAGTGGGGTGGCGTGACGAAAAGGGGGGGGGAGGAACCCCGCGCCGCTACCGACACCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACa                                                                                            <  2:660897/237‑1 (MQ=255)
                                                                                                                             gCTGGACGCCAACAATGGGGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACa                                                                                            <  1:52129/237‑1 (MQ=255)
                                                                                                                                   cGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCACACCGCCAACGAACGCGCCACCATTCTGCGCACCTGGTTCAATTTGATGATGGACCATCAGGACGATTTAGGGCGCCTGATGCCCCTCGAACAGGGTacac                                                                                    >  1:437305/1‑239 (MQ=255)
                                                                                                                                       aacaaTGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTCCCCGCCTGGCCCGCGCTCACCGCCAAAGACCGCGCCACCATTCTGCGCCACTGGTTCAATTTGATGATGGACCATCAGGACCTTTTTGCGCGCCCGATGCCCCCCGAAAAGGGAAAACCAcc                                                                                >  2:804267/1‑238 (MQ=255)
                                                                                                                                        acaaTGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACt                                                                                >  1:80950/1‑238 (MQ=255)
                                                                                                                                                 aaGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCTTCAGGACGATTTAGCGCGCCTGTTGACCCTCGAACAGGGTAAACCACTGGCCGCAGCg                                                                      >  1:665026/1‑239 (MQ=255)
                                                                                                                                                      aTCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCCTCAGGACGGTTTAGCGCGCCGGATGACCCCCGAAAAGGGTAAACCACTGGCCGAAGCGAAAgg                                                                 >  2:117483/1‑239 (MQ=255)
                                                                                                                                                      aTCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACGGGCCGAAGCGAAAg                                                                  >  1:704955/1‑238 (MQ=255)
                                                                                                                                                      aTCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCCCGAACAGGGTAAACCCCCGCCCCAAACGACAgg                                                                 >  2:179785/1‑239 (MQ=255)
                                                                                                                                                        cGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCg                                                               >  1:583138/1‑239 (MQ=255)
                                                                                                                                                                        ggagccgggacaatgtggggatgCGTGCCAAAAGTGGCCGCGGAGAAAACCCGCGCCGTTACGGAGCCCGCGACCCGCGCCCTGCCCGCCTGGCGCGCGCCCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  2:80950/217‑1 (MQ=255)
                                                                                                                                                                        gcaaaCGGCGACAAGCGGGGGAGCGTGCCGAAAATGGGCGCAGAGTAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTAACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  1:804267/236‑1 (MQ=255)
                                                                                                                                                                        gCGACCGGGCACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTACCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  1:49593/239‑1 (MQ=255)
                                                                                                                                                                        gCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTACCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATCACCCTCGAACAGGGGAAACCACTCGCCGAAGCGACAGGGGAAAACaaccaagacc                                                >  2:874539/1‑229 (MQ=255)
                                                                                                                                                                        cgcaccgccGACAAGTTGGGTGTGGTCCGAAAAATGGGCCGGGATGAAACCCGGCGCCGTACCGCCGCCGCCACCCCCGCCCTGCCCGCCGGGCGCCCCCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  2:841035/231‑1 (MQ=255)
                                                                                                                                                                         ggAACGGCCACAAGCGGGGTAGCGTGCCCAAAATGGGCGCGGATGAAACCCGCGCCCCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  1:234719/237‑1 (MQ=255)
                                                                                                                                                                         cGACCGGCGAAACGCTGGGTAGAGGGCCAAAACGGGGCGGGGAGGAAACCCGCGCCGTTGTGGACCCAGCCAACCCCGCCCGGCCCGCCTGGCGCGCGCTACCCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  2:437305/238‑1 (MQ=255)
                                                                                                                                                                                cGAAAAGCTGGGGAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGTTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCttt                                        <  1:282306/238‑1 (MQ=255)
                                                                                                                                                                                               cGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                <  1:509211/151‑1 (MQ=255)
                                                                                                                                                                                               cGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                >  2:509211/1‑151 (MQ=255)
                                                                                                                                                                                                 tGCCGAAAATGGGCCCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCCCCCCCCGCCAAAGAACCGGCCACCCTTCTGCCCAACTGGTTCCATTTGAGTATGGAGCATCAGGACGCTCTAACCCGCCTGATGACCCCCGAACAGGGAAAACTACTGGCCGCAGCGTAAGGCGAAACCAGCTACCACGCCCCCttctttggggggggtgcagaa                      >  2:286157/1‑220 (MQ=255)
                                                                                                                                                                                                                  gATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGCTCAATTTGATGACGGAGCATCAGGACGATTTAGCGCGCCTGATGCCCCTCGAACAGGGTATACCACTGACCGAAGCGGAAGGCGAAATCAGCTACGCCCCCTCCCTTATTTAGTGGTTGGCCAAAAAAAGGAAACGCAATTa     >  2:198123/1‑239 (MQ=255)
                                                                                                                                                                                                                   aTGAACCCCGCGCCGCCACGCGCCCCGCCAACCGCCCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  2:485835/196‑1 (MQ=255)
                                                                                                                                                                                                                        accccggcccgttttCGACGCCCCCACCCGCGCCCTGCCCGCCTGGCGCCCGCTCACCGCCAAAGAACGCGCCACCAGTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTGTATTGAGTGGTTTGCCGAAGAAGGCAATCGCATTTATgg  <  1:278148/223‑1 (MQ=255)
                                                                                                                                                                                                                              cgccgcTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCt                                                                                                                                                                                         <  2:57212/48‑1 (MQ=255)
                                                                                                                                                                                                                              cgccgcTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCt                                                                                                                                                                                         >  1:57212/1‑48 (MQ=255)
                                                                                                                                                                                                                              cgccgcTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTAtt                                     >  1:485835/1‑196 (MQ=255)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                         
GAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCCGCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGG  >  NC_000913/2791202‑2791653

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