Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 628,868 T→G 25.4% G106G (GGT→GGG)  entA → 2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913628,8680TG25.4% 110.0 / 2.9 53G106G (GGT→GGG) entA2,3‑dihydro‑2,3‑dihydroxybenzoate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (9/27);  new base G (14/2);  total (23/30)
Fisher's exact test for biased strand distribution p-value = 4.00e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTG  >  NC_000913/628642‑629105
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
aaGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATTGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGTGGGTAACGTCGGCGGGGCGGTTTAccct                                                                                                                                                                                                                                      >  1:8431/1‑237 (MQ=255)
     aCAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGGCCGATCAGCTCAGTAAAGAGGGCTGGCAGCAGACTTTTGCGGGTAACGTCGGCGGGGCGGTTTAACTGT‑Tcc                                                                                                                                                                                                                                 >  2:1151462/1‑238 (MQ=255)
            ttGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGTGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAAcc                                                                                                                                                                                                                        >  1:33815/1‑239 (MQ=255)
                  aaGCGTTCCCTCAGGAGCAATATCCCTTTGCGACCCGAGTGGTGGGTGTTGCCGACGCTGAGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGGAACGGAGGGACTGGGCGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAAc                                                                                                                                                                                                                  <  2:841032/239‑1 (MQ=255)
                   aGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGGGCGACCGATCAGCTCAGTAAAGAGGGCTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGGGCGCGTAACCTGT‑TCCAG‑CAAACCATTAAc                                                                                                                                                                                                                  >  2:212128/1‑238 (MQ=255)
                   aGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCCGCTCAGTAAAGAGGACTGGCAGCAGACTTTTTCGGGTAACGTCTGCGGGGCGTTTTACCTCT‑TCCCG‑CAAACAATGaaac                                                                                                                                                                                                                 >  2:803061/1‑237 (MQ=255)
                   aGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGGGGGGTTAACCTGG‑TCCAG‑CAAAACATGGAcc                                                                                                                                                                                                                 >  2:538480/1‑239 (MQ=255)
                    gCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTTCGGTTAACGTCGGCGGTGCGTTTAAACTGT‑TCCAG‑CCAAACATGAACCa                                                                                                                                                                                                                >  1:989680/1‑239 (MQ=255)
                    gCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCa                                                                                                                                                                                                                >  1:556964/1‑239 (MQ=255)
                              aGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTTCGGTTAACGTCGGCGGGGGGTTTAACCCGT‑TCCCG‑CAAAACATGAACCA‑GTTTTGCCg                                                                                                                                                                                                      >  1:277169/1‑238 (MQ=255)
                                      tatCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCACTGGGGCGACCGATCAGCTTCGTAAAGAGGACTGGCAGCCGACTTTTGCGGTTACAGGCGGGGGGGGGGTTTAGCTGTATCCAG‑CAAACAACACAACA‑TTTTCG‑CGTCCgcgggg                                                                                                                                                                                             >  2:833563/1‑239 (MQ=255)
                                        tCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAgcggggcg                                                                                                                                                                                           <  1:1125617/239‑1 (MQ=255)
                                        tCCCTGTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAgcggggcg                                                                                                                                                                                           <  2:255465/239‑1 (MQ=255)
                                             ttGCGACCGAGGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAACCGGAGCGGCGGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGGGTGATCGATCCGCTCAGTAAAGAGGACTGGCAGCGGACCTTTGCGGTTAACGTCGCCCGGGCGGTTAACCCGT‑TCCCGCCAACCCTTGAACC‑‑GGTTCGCCGCCCCCCgggccgggcg                                                                                                                                                                                      >  2:546850/1‑239 (MQ=255)
                                              tGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGGGGGTTTAAACTGT‑TTCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGgggcggggcgg                                                                                                                                                                                     >  1:779810/1‑238 (MQ=255)
                                                  aCCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGGATTTTACGCATGGGCGCGGCCCCTCAGCTCAGTCAAGAGGACCGGCCGCCGGCTTTTGCGGTTTACGTCGGCGGGGCGTTTTACCCTT‑TCCCG‑CCAACCACGAACCAGGTTTCCCCG‑CCGccgggcgggggcggtgt                                                                                                                                                                                 >  2:64806/1‑232 (MQ=255)
                                                   ccGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGt                                                                                                                                                                                 <  1:607877/238‑1 (MQ=255)
                                                      aaGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAGCGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGGAAAAAGGACTGGCAGCAGACTTTTGCGGTTAACGTCCGCGGGGGGGTTTACCCGT‑TCCAG‑CAAACCATGAACCA‑GTTTTGCCCTCCGcgggggggggggggtgtcgcc                                                                                                                                                                             >  2:1110695/1‑228 (MQ=255)
                                                                   tGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGCATTTTACGCCTGGGCGCGACCGATCAGCCCAGTAAAGAGGACTGGCAGCCCACTTTTTCGGTTAACGTCCGCGGTGCGTTTTACCTCT‑TCCCC‑CCAACCCTGAACCC‑GTTTCGCCGTCAACGGGGCGGGGGGGCTGTCACTGTGGccttcgcc                                                                                                                                                                 >  2:1027235/1‑231 (MQ=255)
                                                                                     cgcgCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGt                                                                                                                                               <  1:1125612/238‑1 (MQ=255)
                                                                                        gcAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTAtt                                                                                                                                            <  2:1189851/238‑1 (MQ=255)
                                                                                        gcAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGGTAACGTCGGCGGGGCGTTTAACCTGT‑TCCAG‑CAAAACATGAACCA‑GGTTTGCCCTCAGCGGGGCGGGGGGGTTGTTCCTGTGGCGTCCGACCCCCCGCACCCGCcccggaatg                                                                                                                                           >  1:39395/1‑234 (MQ=255)
                                                                                          aaGTGTGTCAGCGACTGTTAGCTGAAACGGGGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCATCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGc                                                                                                                                         <  1:546850/239‑1 (MQ=255)
                                                                                             tgtgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGGTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGGTTGTTACTGTGGCGGCCGGCGCCGCGCACACCCCGCGGATTTGCAtt                                                                                                                                      >  1:793689/1‑238 (MQ=255)
                                                                                              gtgtCAGCGACTGTTAGCTGAAGCGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATg                                                                                                                                      <  1:1110695/238‑1 (MQ=255)
                                                                                              gtgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCAt                                                                                                                                       <  1:1151462/237‑1 (MQ=255)
                                                                                              gtgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATg                                                                                                                                      <  1:746339/238‑1 (MQ=255)
                                                                                               tgtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCCGCAGACTTTTGCGGGTAACGTCCGCGGGGCGGTTTACCCTT‑TCCCG‑CAAACCATGACACA‑GTTTTGCCCGCCGCGGGGGGGGGCGGCTTTTCCTTTGGGGCCCCGCGCCCCCCACACGCCGCGTACTgggcgtgg                                                                                                                                    >  2:743143/1‑233 (MQ=255)
                                                                                                gtCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAg                                                                                                                                    <  1:61539/238‑1 (MQ=255)
                                                                                                                                  cgcgCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGc                                                                                                  <  1:1232121/238‑1 (MQ=255)
                                                                                                                                          tCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGtt                                                                                          <  2:556964/238‑1 (MQ=255)
                                                                                                                                                cggcggGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGt                                                                                    <  2:823758/238‑1 (MQ=255)
                                                                                                                                                   cggGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCgg                                                                                 <  1:429152/238‑1 (MQ=255)
                                                                                                                                                         ttttACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGaa                                                                          <  1:128138/239‑1 (MQ=255)
                                                                                                                                                            tACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTg                                                                       <  2:33815/239‑1 (MQ=255)
                                                                                                                                                                 aTGGGCGCGACCGATCAGCTCAGTAAAGCGGACTGGCAGCAGACTTTTGCGGCTAACGTCTGCGGGGCGTTTAACCTGT‑TCCAG‑CAAACCATGGACCA‑GGTTCGCCGTCAGCGGGGCGGGGCGCTTGTCACTGTGGCGGCCGGAGCCGCGCACCCGCCGCCGATTGGGATGAGTGCTTATGGGGCCTGGGAAGGGGCGCCGACAAAGCTGGCTGTGCGCGGCGCGGCGGAACTgggcgg                                                                  >  2:885983/1‑235 (MQ=255)
                                                                                                                                                                 aTGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCggg                                                                  >  1:582694/1‑239 (MQ=255)
                                                                                                                                                                   gggCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCCCCCCCGCGCCCCAGCCCCGCATTGGCCTGAGTGCCTTTGGCGCCTCGCAACGCGCCCTCAAAAGCCTGGCGTTGCGCGTTGGGCCGGGACTGGCCGGtg                                                                >  2:499016/1‑238 (MQ=255)
                                                                                                                                                                           ccGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTg                                                        <  1:99320/239‑1 (MQ=255)
                                                                                                                                                                           ccGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTg                                                        <  1:1144763/239‑1 (MQ=255)
                                                                                                                                                                           ccGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTg                                                        <  1:212128/239‑1 (MQ=255)
                                                                                                                                                                             gATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTgcg                                                      <  1:784597/239‑1 (MQ=255)
                                                                                                                                                                                             agGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGc                                                                               >  1:501306/1‑198 (MQ=255)
                                                                                                                                                                                             agGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGc                                                                               <  2:501306/198‑1 (MQ=255)
                                                                                                                                                                                             agGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGGTGCGGCGGGGGCCGTCAAGTGGgttt                                       >  1:461220/1‑237 (MQ=255)
                                                                                                                                                                                                         agtcTTTTGCGGTTAGCGGCGGCGGTGCGTTTAGCCTGT‑GCCAG‑CAAACCATGACCCA‑GTTGCGCCGTCAGCGGGGCGGGGTGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTcc                             <  2:8431/233‑1 (MQ=255)
                                                                                                                                                                                                            cTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccga                        <  1:538480/238‑1 (MQ=255)
                                                                                                                                                                                                             ttttGCGGTTAACGTCGGCGGGGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccgaca                      <  2:779810/239‑1 (MQ=255)
                                                                                                                                                                                                             ttgtGCGGTTAACGGCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccgaca                      <  2:39395/236‑1 (MQ=255)
                                                                                                                                                                                                             tgttGCGGTTAACGTCGGCGGGGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCcaccgaca                      <  1:885983/237‑1 (MQ=255)
                                                                                                                                                                                                                          gTCGGCGGAGCGTATAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAacg         <  1:245435/239‑1 (MQ=255)
                                                                                                                                                                                                                            cggcggTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAacgca       <  1:64806/239‑1 (MQ=255)
                                                                                                                                                                                                                               cggTGCGTTTAACCTGT‑TCCAG‑AAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAacgcacg     <  1:383808/238‑1 (MQ=255)
                                                                                                                                                                                                                                  tGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCtg  <  1:651070/238‑1 (MQ=255)
                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                
AAGTTACAGGTTTTGATCAAGCGTTCACTCAGGAGCAATATCCCTTTGCGACCGAAGTGATGGATGTTGCCGACGCTGCGCAGGTCGCGCAAGTGTGTCAGCGACTGTTAGCTGAAACGGAGCGACTGGACGCGCTGGTCAATGCGGCGGGAATTTTACGCATGGGCGCGACCGATCAGCTCAGTAAAGAGGACTGGCAGCAGACTTTTGCGGTTAACGTCGGCGGTGCGTTTAACCTGT‑TCCAG‑CAAACCATGAACCA‑GTTTCGCCGTCAGCGGGGCGGGGCGATTGTCACTGTGGCGTCCGACGCCGCGCACACGCCGCGTATTGGCATGAGTGCTTATGGCGCATCGAAAGCGGCGCTGAAAAGCCTGGCGTTGAGCGTCGGGCTGGAACTGGCGGGTAGCGGCGTGCGCTGTAATGTGGTTTCGCCTGGCTCCACCGACACCGATATGCAACGCACGCTG  >  NC_000913/628642‑629105

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