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

GAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGG  >  NC_000913/2791202‑2791653
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                         
gAGCGGCACGTAGTGTGGGTCCCTTACACGCCGGATTTAATCAATAACCTTTGAAAACAGGATGGAGCGATGAAACTTAACGACAGTAACTTACTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGAcgccgc                                                                                                                                                                                                                       <  1:882477/239‑1 (MQ=255)
                  aTGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGCAAATGGGCGCGGCTGAAACCC‑GCGCCGCTATCAACCCCGCCAACCGCGCCCTGcccc                                                                                                                                                                                                      >  1:858234/1‑237 (MQ=255)
                    gCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGCCGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGCAAATGGGCCCGGATGAAACCC‑GCCCCGCTATCGACCCCCCCAAccccccccccccccccc                                                                                                                                                                                                   >  2:379881/1‑230 (MQ=255)
                      cTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCCCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCCGGCCCTCCCCCCCTg                                                                                                                                                                                                  >  1:792071/1‑238 (MQ=255)
                       ttACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGAGCGCGGATGAAACCC‑GCCCCCCTATCGACGCCCCCAACCGCGCCCCGCCCCCCTGgc                                                                                                                                                                                                >  2:52129/1‑239 (MQ=255)
                       ttACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGAAGCGTGCCGAAAATGGCCGCGGATGAAACCC‑GCGCCGCTATCGACGGCGCCAACCCCGCCCCGCCCGCCGGgc                                                                                                                                                                                                >  2:780985/1‑239 (MQ=255)
                         acacGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGCAACTTATTCCGCCAGCAGGCGTTGATTAACGAGGAATGGCTGGACGCCAACAATGGTGAAGCCATCCACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCCCCGCTATCGACCCCGCCAACCTCGCCCCGCCCCCCTGgggc                                                                                                                                                                                              >  2:642216/1‑239 (MQ=255)
                              ccgcATTTAACCAATAACCTTTGAAAACGGGATGTAGCAAGGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCt                                                                                                                                                                                         <  1:642216/239‑1 (MQ=255)
                               cgcATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGTCGAACGGCGACCAGCTGGGAAGCGTGCCGAAAAAGGCCGCGGATGAAACCC‑GCCCCCCCACCGCCCCCGCCCACGCCGCCCCGcccccccgccccccccc                                                                                                                                                                                         >  2:637933/1‑227 (MQ=255)
                               cgcATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACCCCGCCAACCGCGCCCTGCCCTCCTGGCGCGCGCt                                                                                                                                                                                         >  1:267575/1‑238 (MQ=255)
                               cgcATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGACGCGGCTGAAACCC‑GCGCCGCTATCCACCCCCCCACCCcccccctgcccccctgccccccgcc                                                                                                                                                                                         >  2:87687/1‑229 (MQ=255)
                                gcATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCa                                                                                                                                                                                       <  1:235974/239‑1 (MQ=255)
                                             aaCCTTTGAAAACGGGTTTGGGCGATGAACATTAACGGACGTAACTTTTGCCGCCGGAGGGAGTTGATTAACGGGGAATGGCTGGACGCCAACAAGGGTGAAGCCATGGACGTCCCCAATCCGGCGAACGGCGCCAAGCTGGGTCGCGTCCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAAcg                                                                                                                                                                          <  2:376629/239‑1 (MQ=255)
                                                  ttaaaaaCAGGATGTAGGGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCGCCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGcca                                                                                                                                                                     <  1:204029/236‑1 (MQ=255)
                                                  ttGAAAACAGGATGTAGCGATGAAACTTAACGGCAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGcca                                                                                                                                                                     <  1:379881/239‑1 (MQ=255)
                                                  ttGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCCCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGCCGCGCGCCCACCGCCAAAGAACGCGcc                                                                                                                                                                      >  1:553666/1‑238 (MQ=255)
                                                   tGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCCCCGGGCGCGGGTTCACCGCCAAAGAACGCGccc                                                                                                                                                                     >  1:817230/1‑237 (MQ=255)
                                                               gTACCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCCCCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCaa                                                                                                                                                        <  1:466534/239‑1 (MQ=255)
                                                                   cGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAACATGGGCGCGGATGAAACCC‑GCCCCGCTATCGACGCCGCCCACCGGGCCCTGCCCGCCGGCACCGGCCCCACCGCCAAAGACCGCGcccccctcccccccacctgg                                                                                                                                                    >  2:72006/1‑224 (MQ=255)
                                                                    gATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGA‑TAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCCCCATTCTGCGCAACTGGt                                                                                                                                                   >  1:790975/1‑238 (MQ=255)
                                                                      tGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCCTTCcgcgacaccggtc                                                                                                                                                  >  1:347547/1‑229 (MQ=255)
                                                                         aactttaacgcggTAACTTATTCCGCCGGCGGGCGTTGATTAAGGGGGAAGGGCTGGACTCCAAAAAGGGGGAAGCCATTGACGTCACCAATCCGGCGAACGGCGACAAGCGGGGTAGCGTGCCGAAAATGGGCCCGGATGAAACCC‑GCGCCGTTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCaa                                                                                                                                               <  2:598428/226‑1 (MQ=255)
                                                                                cGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTtgataa                                                                                                                                       >  1:814532/1‑237 (MQ=255)
                                                                                   cAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCa                                                                                                                                                                                 <  1:31953/194‑1 (MQ=255)
                                                                                   aaGTAACTTATTCCGCCAGCGGGGTTGGATTAAGGGGGAAGGGCTGGAGCGCAAAAGTGGGGAGGCATTCGACGTCCCCAATGCGGCGAGCGGGGACAAGGGGGGTGGCGGGCCGAAAATGGGGGCGGATGAACCCC‑GGGCCGCTACCGACGCCGCCAACCGCTCCCCGCCCGCCTGGCGCGCGCTTACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGatgatg                                                                                                                                     <  2:40782/237‑1 (MQ=255)
                                                                                        aCTTATTCCGCCAGCGGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCAt                                                                                                                               <  1:72006/239‑1 (MQ=255)
                                                                                            aTTCCGCCAGAAGGCGTTGATAAACGGGAAATGGCGGGACGCCAAAAATTGGGAAGCCATCGACGGCACCAATCCGGCGAACGGCGACAAGCGGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAgg                                                                                                                           <  2:550581/239‑1 (MQ=255)
                                                                                              tCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACCCCGCCAACCGCGCCCTACCCGCCTGGCCGCCGCCCACCCCCAAAGAACGCCCCACCACTCCCCGAAACCGGTTCACttacttgcaaaggaacaagaac                                                                                                                         >  2:282306/1‑219 (MQ=255)
                                                                                                       cggtgTTGATTAACGGGGAATGGCCGGCCGCCAACAAGGGGGAGGCCAGCGGCGTCACCAATCGGGCGAACGGCGACAAGCGGGGTGGCGTCCCGAAAATGGGCGCGGAGGAAACCC‑GCGCCGGTAGCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                <  2:675592/235‑1 (MQ=255)
                                                                                                       agggagggtttAACGGGGAATGGGTGGACGCCAACAATGGTGAACGCCTGGACGTAACCAACCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAAGGGGCCCGGAGGAACCCC‑GCGCCGCTACCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                <  2:384139/230‑1 (MQ=255)
                                                                                                          cGTTGATTAACGGGGAACGGCGGGCCGCCAACAATGGTGAAGCCATGGACACCCCCAATCCGGCAAACGGCGCCAAGCTGGGTAGCGTGCCGAAAATGGGCGCGAATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                <  2:267575/236‑1 (MQ=255)
                                                                                                               aTTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCAAAAGAACCCGCCACCATTCTCCGCAACTGGTTCAATTTGATGATGGAGCATCGGGACGATTTAGCGCGCCtgat                                                                                                        >  1:134727/1‑239 (MQ=255)
                                                                                                                ttAACGGGGAATGGCTGGACGCCAACAAGGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCtgat                                                                                                        <  1:87687/238‑1 (MQ=255)
                                                                                                                         aaTGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGAccggtc                                                                                                >  1:384139/1‑233 (MQ=255)
                                                                                                                         aaTGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGAGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCCCACCGCCAAACAACGCGCCACCATTCTGCGCACCTGGGTCAATTTGATGACGGAGCATCAGGACGATTTTCACCGCCCGATGACCCTcta                                                                                               >  2:77340/1‑236 (MQ=255)
                                                                                                                          ctggctggcCCCCAAAATTGGGGAAGCCATCGCGGCCCCCAATCCGGTGAACGGCGACAAGCGGGGTAGCGTGCCGAAAAGGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCTCCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGaa                                                                                              <  2:792071/237‑1 (MQ=255)
                                                                                                                          aTGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGaa                                                                                              <  1:780985/238‑1 (MQ=255)
                                                                                                                          aTGGCGGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAAGGGGCGCGGATGAACCCC‑GCGCCGCTATCGACGCCGCCAACCGCCCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGaa                                                                                              <  1:824870/238‑1 (MQ=255)
                                                                                                                          aTGGCGGGACCCAAAAAAGGGGGAAGCCCCCGCAGTCAGCAAGCCGGCAGACGGCGACAAGGGGGGGAGCGGGCCGAAACTGGGCGCGGATGAAACCC‑GCGCCGCTAGCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCAGCATGCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGaa                                                                                              <  2:858234/238‑1 (MQ=255)
                                                                                                                             gCTGGAGCCAAACAAGGTAGAAGCCAGCGACGGCACTACGCCCTGGAACCGCGACAAGTGGGGTGGCGTGACGAAAAGGGGGGGGGAGGAACCCC‑GCGCCGCTACCGACACCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACa                                                                                            <  2:660897/237‑1 (MQ=255)
                                                                                                                             gCTGGACGCCAACAATGGGGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACa                                                                                            <  1:52129/237‑1 (MQ=255)
                                                                                                                                   cGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCACACCGCCAACGAACGCGCCACCATTCTGCGCACCTGGTTCAATTTGATGATGGACCATCAGGACGATTTAGGGCGCCTGATGCCCCTCGAACAGGGTacac                                                                                    >  1:437305/1‑239 (MQ=255)
                                                                                                                                       aacaaTGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTCCCCGCCTGGCCCGCGCTCACCGCCAAAGACCGCGCCACCATTCTGCGCCACTGGTTCAATTTGATGATGGACCATCAGGACCTTTTTGCGCGCCCGATGCCCCCCGAAAAGGGAAAACCAcc                                                                                >  2:804267/1‑238 (MQ=255)
                                                                                                                                        acaaTGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACt                                                                                >  1:80950/1‑238 (MQ=255)
                                                                                                                                                 aaGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCTTCAGGACGATTTAGCGCGCCTGTTGACCCTCGAACAGGGTAAACCACTGGCCGCAGCg                                                                      >  1:665026/1‑239 (MQ=255)
                                                                                                                                                      aTCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCCTCAGGACGGTTTAGCGCGCCGGATGACCCCCGAAAAGGGTAAACCACTGGCCGAAGCGAAAgg                                                                 >  2:117483/1‑239 (MQ=255)
                                                                                                                                                      aTCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACGGGCCGAAGCGAAAg                                                                  >  1:704955/1‑238 (MQ=255)
                                                                                                                                                      aTCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCCCGAACAGGGTAAACCCCCGCCCCAAACGACAgg                                                                 >  2:179785/1‑239 (MQ=255)
                                                                                                                                                        cGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCg                                                               >  1:583138/1‑239 (MQ=255)
                                                                                                                                                                        ggagccgggacaatgtggggatgCGTGCCAAAAGTGGCCGCGGAGAAAACCC‑GCGCCGTTACGGAGCCCGCGACCCGCGCCCTGCCCGCCTGGCGCGCGCCCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  2:80950/217‑1 (MQ=255)
                                                                                                                                                                        gcaaaCGGCGACAAGCGGGGGAGCGTGCCGAAAATGGGCGCAGAGTAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTAACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  1:804267/236‑1 (MQ=255)
                                                                                                                                                                        gCGACCGGGCACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTACCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  1:49593/239‑1 (MQ=255)
                                                                                                                                                                        gCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTACCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATCACCCTCGAACAGGGGAAACCACTCGCCGAAGCGACAGGGGAAAACaaccaagacc                                                >  2:874539/1‑229 (MQ=255)
                                                                                                                                                                        cgcaccgccGACAAGTTGGGTGTGGTCCGAAAAATGGGCCGGGATGAAACCCGGCGCCG‑TACCGCCGCCGCCACCCCCGCCCTGCCCGCCGGGCGCCCCCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  2:841035/231‑1 (MQ=255)
                                                                                                                                                                         ggAACGGCCACAAGCGGGGTAGCGTGCCCAAAATGGGCGCGGATGAAACCC‑GCGCCCCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  1:234719/237‑1 (MQ=255)
                                                                                                                                                                         cGACCGGCGAAACGCTGGGTAGAGGGCCAAAACGGGGCGGGGAGGAAACCC‑GCGCCGTTGTGGACCCAGCCAACCCCGCCCGGCCCGCCTGGCGCGCGCTACCCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  2:437305/238‑1 (MQ=255)
                                                                                                                                                                                cGAAAAGCTGGGGAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGTTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCttt                                        <  1:282306/238‑1 (MQ=255)
                                                                                                                                                                                               cGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                <  1:509211/151‑1 (MQ=255)
                                                                                                                                                                                               cGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAgcg                                                                                                                >  2:509211/1‑151 (MQ=255)
                                                                                                                                                                                                 tGCCGAAAATGGGCCCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCCCCCCCCGCCAAAGAACCGGCCACCCTTCTGCCCAACTGGTTCCATTTGAGTATGGAGCATCAGGACGCTCTAACCCGCCTGATGACCCCCGAACAGGGAAAACTACTGGCCGCAGCGTAAGGCGAAACCAGCTACCACGCCCCCttctttggggggggtgcagaa                      >  2:286157/1‑220 (MQ=255)
                                                                                                                                                                                                                  gATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGCTCAATTTGATGACGGAGCATCAGGACGATTTAGCGCGCCTGATGCCCCTCGAACAGGGTATACCACTGACCGAAGCGGAAGGCGAAATCAGCTACGCCCCCTCCCTTATTTAGTGGTTGGCCAAAAAAAGGAAACGCAATTa     >  2:198123/1‑239 (MQ=255)
                                                                                                                                                                                                                   aTGAACCCC‑GCGCCGCCACGCGCCCCGCCAACCGCCCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTAcgccgc                                               <  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)
                                                                                                                                                                                                                                           |                                                                                                                                                                                                                         
GAGCGGCACGTAGTGTGGATGCCTTACACGCCGCATTTAATCAATAACCTTTGAAAACAGGATGTAGCGATGAAACTTAACGACAGTAACTTATTCCGCCAGCAGGCGTTGATTAACGGGGAATGGCTGGACGCCAACAATGGTGAAGCCATCGACGTCACCAATCCGGCGAACGGCGACAAGCTGGGTAGCGTGCCGAAAATGGGCGCGGATGAAACCC‑GCGCCGCTATCGACGCCGCCAACCGCGCCCTGCCCGCCTGGCGCGCGCTCACCGCCAAAGAACGCGCCACCATTCTGCGCAACTGGTTCAATTTGATGATGGAGCATCAGGACGATTTAGCGCGCCTGATGACCCTCGAACAGGGTAAACCACTGGCCGAAGCGAAAGGCGAAATCAGCTACGCCGCCTCCTTTATTGAGTGGTTTGCCGAAGAAGGCAAACGCATTTATGG  >  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: ‑