Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 650,331 T→G 12.8% T54P (ACC→CCC)  citE ← citrate lyase, citryl‑ACP lyase (beta) subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913650,3310TG12.8% 229.6 / 2.4 93T54P (ACC→CCC) citEcitrate lyase, citryl‑ACP lyase (beta) subunit
Reads supporting (aligned to +/- strand):  ref base T (15/58);  new base G (16/4);  total (31/62)
Fisher's exact test for biased strand distribution p-value = 1.67e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTA  >  NC_000913/650095‑650560
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        
ggTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGgtgt                                                                                                                                                                                                                                        <  1:54673/239‑1 (MQ=255)
        cGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCa                                                                                                                                                                                                                                <  1:1147189/239‑1 (MQ=255)
        cGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTGTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCCCCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCa                                                                                                                                                                                                                                <  1:1242068/239‑1 (MQ=255)
         gggTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGACCATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCCTTCGGAATCCAGCGCGTTGCCACGCAAAATGGATTCAATATCCCGCTACAGCGGATGTTGCAGCGCGTGGAAAACAATGGGGCGGGGGGGGGTCTTTTACa                                                                                                                                                                                                                                >  2:371115/1‑239 (MQ=255)
         gggTTCACGACCAAAGGCTTTTTCGGTTCGGAGGATCTTTTTTTCAATTCCCAAAACCTCCCGGGGGGTTTCGGTTTTCGGCAGACGCAAAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACCCCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCa                                                                                                                                                                                                                                <  2:64836/238‑1 (MQ=255)
          gttcccccccccccggttttttCGTTTCGGAGGATCTCTTTTTCAACATCCAGAAGATCCTGAGGGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCCCCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGGGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACg                                                                                                                                                                                                                              <  2:688688/222‑1 (MQ=255)
          ggttccccgccccatggttttttGATACCGAGGATTCCTTTTTCAATGTCCAGAACATCCGGAGCGGTTTCGGTTTTCGGCAGACGCACAAGGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACg                                                                                                                                                                                                                              <  2:1189105/222‑1 (MQ=255)
          ggTGCACGACCACGGGCTTTTGCGAATGGCAGGAGTTCTTTTTCAATAGCCAGAACATCCTGAGCGGGATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACg                                                                                                                                                                                                                              <  2:937968/239‑1 (MQ=255)
           gTTCCCGACCACAGTCTTTTGCGATACGCAGGATCTCTTTTCCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATGTGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACg                                                                                                                                                                                                                              <  1:54349/238‑1 (MQ=255)
           gTTCCCGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACg                                                                                                                                                                                                                              <  1:58688/238‑1 (MQ=255)
           gTTCAGCACCACAGGCTTTTTCGAACCGACGGATTTCTTTTTCAATTACCAGAACATCCTGAGCGGTTACGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACg                                                                                                                                                                                                                              <  1:257363/238‑1 (MQ=255)
           gTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACg                                                                                                                                                                                                                              <  1:849004/238‑1 (MQ=255)
             tcccgccccccggtttttttGATTCGCAGGATCTCTTTTTCAATTACCAGAACATCCTGGGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCa                                                                                                                                                                                                                            <  2:383276/224‑1 (MQ=255)
             tcccgcccccaGGCTTTTTCGTTACGCAGGATTTCTTTTTCAATATCCCGACCATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCa                                                                                                                                                                                                                            <  2:583249/229‑1 (MQ=255)
              cACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGCTACACCGGAGGTTGAAGCGCGTGGTAAACCATGCGGC‑GGGGGGGGCCTTTTTCACCCaa                                                                                                                                                                                                                           >  1:9715/1‑238 (MQ=255)
                      caGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTACCACCCCATTCGGAATCCAGCGCGTTGCCACGCACACTGGTTTCAATATCGCGTTTCAGCTGCTGTTGCGGCCCGTGGTAAACCATACCGC‑GGGGGGGGGCTTTTTTAccccaaactacgcg                                                                                                                                                                                                                  >  2:888254/1‑226 (MQ=255)
                      caGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGTCTTTTTCACGCAATGcatcggg                                                                                                                                                                                                                  >  1:1135655/1‑233 (MQ=255)
                         gCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACACACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGCAACGCCCAATGGTTTCAATATCCCGATACAGCGGAGGTGGCCGCGGGGGGTAAACCAGTGGGC‑GGGGGGGGCATTTTTAAccacagggtacgggggg                                                                                                                                                                                                               >  2:531470/1‑223 (MQ=255)
                             tttCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGGGGGGTCTTTTTCCCGCAATCGTTCGGGGTCt                                                                                                                                                                                                              >  1:699074/1‑236 (MQ=255)
                             tttCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTTACGCAATGGTACGGGGTCTTcc                                                                                                                                                                                                           >  1:373204/1‑238 (MQ=255)
                              tttggttcgcatgtgttttttttCAATATCCATAACATCCTGAGCGGTTTCGGTTTTCGGCAGACGCACAACGT‑CCGCACCCCCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGa                                                                                                                                                                                                          <  2:852751/221‑1 (MQ=255)
                                   tACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGTAATCCAACGCGTTGACACGCACAATGGTGTCAATATCGCGATACAGCGGATGTTGCAGCCCGTGGTAAACCAAGCGGC‑GGGGGTGGGTTTTTTCACGCAAAGcgcacaggacatcaagggga                                                                                                                                                                                                     >  2:198381/1‑220 (MQ=255)
                                   ctcGCAGGTCCTCTTTTTCAAGATCCATAACATCCTGGACTGTATCGGTTTTCGGGTGACGCACAACGC‑TCGCAACCCCGCGAACGGCGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCTCGCAATGCTACGGAGTCTTCGAGGTCa                                                                                                                                                                                                     <  2:1234496/237‑1 (MQ=255)
                                     cGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTa                                                                                                                                                                                                                      <  1:472910/220‑1 (MQ=255)
                                          gATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCACTATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGCCTTTTTCACGCAATTCTACGGAGTCTTCGCGGGCAAACATCa                                                                                                                                                                                              >  1:256454/1‑239 (MQ=255)
                                           aTCTCGTTTTCAATATCCAGACCCTCCGGAGCGGTATCGGTTTTCGGCAGACGCGAAACGT‑CCGCACCACCGCCAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACAt                                                                                                                                                                                                <  1:888254/236‑1 (MQ=255)
                                             ctctTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGCCAAACATCAgg                                                                                                                                                                                            >  1:1210416/1‑238 (MQ=255)
                                               ctTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGGCTTTTTTACGCAACGCTACAGAGACTTCgcccggccctacacaaatc                                                                                                                                                                                         >  2:472910/1‑221 (MQ=255)
                                                tttttCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGGCAAACATCAGGGCa                                                                                                                                                                                         >  1:378789/1‑238 (MQ=255)
                                                   gtCAATATCCAGAACATCCTGAGGGGTAACGGGTTTCCGAAGGCTCAAAACGT‑CCGCACCACCGGGAACGGAGGCTACCAGGTCGGTAACACCCCATGCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAg                                                                                                                                                                                     <  1:1326188/238‑1 (MQ=255)
                                                    tCAATCTCCAGAACCTCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCCCCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGc                                                                                                                                                                                    <  1:1327860/239‑1 (MQ=255)
                                                       catcccaaacccaccTGGGGGGTTTCGGTTTTCGGCAGCCCCACAACGG‑CCGACCCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGGGCGTTGACACGCAAAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCgg                                                                                                                                                                                 <  2:1311377/226‑1 (MQ=255)
                                                       atatCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATAAAGCGGATGTTGCAGCGCGGGGTAAACCATGCGGC‑GGGGGGTGTGTTTTTCCCGCAAAGACACGGGGACTTCGGGGGCAAAAAACAAGGACTCAACCgg                                                                                                                                                                                 >  2:1327860/1‑239 (MQ=255)
                                                        tatCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGAAAACCATCCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTCCGGGGGCTTCGAGGTCCAACATCAGGGCTTAAGCCggg                                                                                                                                                                                >  1:820465/1‑239 (MQ=255)
                                                         agccAAAACATCCTGAGCGGTATCGGTTTTCGGGAGACGCACAAAGT‑CCGGACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGt                                                                                                                                                                               <  1:198381/237‑1 (MQ=255)
                                                          tCCAGACCAGCCTGAGCGGTATCGGTTTCCGGCAGACGCACAACGT‑CCGCTACACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCtt                                                                                                                                                                                                             <  1:177733/208‑1 (MQ=255)
                                                               aaCATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGGCAAAACTAAGGGCACCACCCGGGGAGATGa                                                                                                                                                                         >  1:1274482/1‑239 (MQ=255)
                                                                 caccccgggCGGTTTCCGTTTGCGGCCGACGCACAACGT‑CCGCCACACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAg                                                                                                                                                                       <  2:196541/231‑1 (MQ=255)
                                                                    ccTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTTACGCAATGCTACGGGGCCTTCGGGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAgggg                                                                                                                                                                    >  1:105037/1‑237 (MQ=255)
                                                                    ccTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCa                                                                                                                                                                                                                                >  1:928947/1‑179 (MQ=255)
                                                                    ccTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCa                                                                                                                                                                                                                                <  2:928947/179‑1 (MQ=255)
                                                                    ccTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGTCTTTTTCCCGCAATGCTACGGGGCCTTCGAGGGCAAAAACCAGGGCATCACCCGGGGAAATGAAGGAg                                                                                                                                                                    >  2:580547/1‑239 (MQ=255)
                                                                    ccTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGTCTTTTTCACGCAATGCTACGGAGGCTTCGAGGTCAAACATCAGGGCATCAGCCCGGTAGATGAAGGAg                                                                                                                                                                    >  1:589403/1‑239 (MQ=255)
                                                                    ccTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGGCTTTTTTCCGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGGAGATGAAggga                                                                                                                                                                    >  1:295845/1‑237 (MQ=255)
                                                                            tctccGGTTTCCGG‑AGACCCACAACCTCCCGCACC‑CCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTTAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGAc                                                                                                                                                            <  1:1362634/234‑1 (MQ=255)
                                                                               tCGGTTGTCGGCAAACCCACAACGT‑CCGCACCCCCGCGAAGGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCAt                                                                                                                                                         <  1:92983/239‑1 (MQ=255)
                                                                                   ttttcggaccaccccaaccct‑CCCCACCCCCGGGGACGGCGGCTCCCAGGTCGTTAACACCCCCTTCGGAAGCCAGGCGGTTGACCCGCACAATGGTTTCAATATGGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATcgc                                                                                                                                                      <  2:1359237/218‑1 (MQ=255)
                                                                                   ttttCGGCAGACGCACACCGT‑CCGCACCTCCGCGAACGCCGGCTTCCAGGTCGGTAACCCCCCATGCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATcgc                                                                                                                                                      <  2:485250/238‑1 (MQ=255)
                                                                                    tttCGGCAGACGCACAACGT‑CCGCAACACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATcgcg                                                                                                                                                     <  1:766027/238‑1 (MQ=255)
                                                                                       cGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGGCTTTTTTACCCAATGCTACGGGGCCTTCGGGGTCAAACATCCGGGCATCAGCCGGGGAGATGAAGGAGTTGCTGACCCTCCCGGCAt                                                                                                                                                 >  1:163361/1‑239 (MQ=255)
                                                                                             ccGCACCACGT‑CCGCTCCCCCGCAAACGACGGCTCCCAGGTCGTTAACACCCCATCCGGAATCCAGCGCGTTGACTCGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGt                                                                                                                                                                               <  1:299153/202‑1 (MQ=255)
                                                                                              ccgcaaaCGT‑CCGCCCCACCGAGAAAGACGGCTTCCAGGTCGTTAACCCCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCAc                                                                                                                                           <  1:246534/233‑1 (MQ=255)
                                                                                               cgccaaaaggccccccccccGGGAACGGCGGGTTCCAGGGCGTTAACCCCCCATTCGGAATCCAGCGGGTTGACACGCACAATGGTTTCAATACCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCAc                                                                                                                                           <  2:105037/220‑1 (MQ=255)
                                                                                                    aCGT‑CCGCACCACCGCGAACGGCGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑cac                                                                                                                                   <  1:580547/239‑1 (MQ=255)
                                                                                                     acggccggcccaccGCGAACGGCGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑cac                                                                                                                                   <  1:377500/230‑1 (MQ=255)
                                                                                                          cGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGGGGCTTTTTTACGCAATTCTCAGGGGGCTTCGAGGGCGAAACTCAAGGCATCACCCGGGAAGATGGAAGGGATGCCGACCAACGCGGCAATGGGACCAGGCCacaaaa                                                                                                                               >  2:329667/1‑238 (MQ=255)
                                                                                                                    ggAACGGCGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTgcg                                                                                                                    <  1:495803/238‑1 (MQ=255)
                                                                                                                      gtcgcgggttgccAGGTCTTTACCCCCCCTTTCGGAATGCAGCGCGTTGAACCCGATAATGGTGTCAATATCGCGATACAGCGGATGGTGCAGCGCGTGGGAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTGCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTgcggc                                                                                                                  <  2:75260/228‑1 (MQ=255)
                                                                                                                           cgGTTTCCAGGTCGTTTACCCCCCATTCGGAATCCAGCGCGTTGACCCGGACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTg                                                                                                            <  2:378789/239‑1 (MQ=255)
                                                                                                                               ttCCAGGTCGTTAACCCCCCATTCGGAATCCAGCGCGTTGAAACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCTTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGa                                                                                                         <  1:594522/238‑1 (MQ=255)
                                                                                                                                     gTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGGGCGGGtgtttt                                                                                                  >  1:391859/1‑237 (MQ=255)
                                                                                                                                                cccATTCGGAATCCAGCGCGTTGACACGCCAAATGGTTTCAATATCGCGATACAGCTGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCAc                                                                                                                                           <  2:711017/189‑1 (MQ=255)
                                                                                                                                                cccATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCa                                                                                                                                                                                                     <  1:1160941/131‑1 (MQ=255)
                                                                                                                                                cccATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGGCCATCGCGGCATTGGCACCAGG‑CCCAACCAACATGCTGCGGC‑GGGCTCGGGGTTTAcacttttttgg                                                                                        >  1:2189/1‑228 (MQ=255)
                                                                                                                                                cccATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCAc                                                                                                                                           >  1:711017/1‑189 (MQ=255)
                                                                                                                                                cccATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CCCAAAAAACATGCGCCGGC‑ggggggggggttttctgttttgtcac                                                                                       >  1:1360967/1‑216 (MQ=255)
                                                                                                                                                       ggAATCCAGGGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCgg                                                                               <  1:217082/239‑1 (MQ=255)
                                                                                                                                                             cgccgctTTGAC‑CGCAAAGTGGTTTCAATTCCGGGGTACTGGGGATGGTGCAGGGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCAt                                                                        <  2:295845/236‑1 (MQ=255)
                                                                                                                                                                cgcgTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGGGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCCTCGCGGCATTGGCGCCAGG‑CCCAAACAAACTGCTGCGGCGGGGGGCGGGGTTTACgttttttatcgacagggggaaacatt                                                                        >  2:218959/1‑219 (MQ=255)
                                                                                                                                                                     tGACACGCACAATGGTTTCAAGAGCGGGAGACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTGCACGCAATGCTACGGAGTCTTCGAGGTCAAACAGCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAAt                                                                 <  2:163361/239‑1 (MQ=255)
                                                                                                                                                                     gtgacccccaaagggTTTAAATGTCGGGGTACAGCGGGTGTTGCAGGGCGTGGAAAACCATGGGGG‑GGGCGGGGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAAt                                                                 <  2:820465/226‑1 (MQ=255)
                                                                                                                                                                      gACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGGAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCGCCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAAt                                                                 <  1:218959/238‑1 (MQ=255)
                                                                                                                                                                                atgttttcaattctgctgtcgcgCGGTTGTTGAGGGGCGTGGTAAACCATGGGGC‑GGGCGGTGTGTTTTTCCCGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGc                                                      <  2:816538/217‑1 (MQ=255)
                                                                                                                                                                                aTGGTTTCAATTTCGGGATACAGGGGATGTTGCAGCGCGGGGTAAACCATGCGGC‑GGGCGGTGGCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGc                                                      <  2:656111/239‑1 (MQ=255)
                                                                                                                                                                                  ggTTTCATTTGCGGGATAAAGGGGATGTTGGAGCGGTTGGTAACCCAGGCGGC‑GGGCGGGTTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCag                                                    <  2:361522/239‑1 (MQ=255)
                                                                                                                                                                                  ggTTTCAATTGCGCGATACCGGGGTGTTTGCCGCGCGTGGTAAACCCTGCGGC‑GGGCGGTGTTTTTTTCACGCAATGCTACGGAGGCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCag                                                    <  2:1274482/239‑1 (MQ=255)
                                                                                                                                                                                     ttcaagtccgcagtACGGCGGATGTGGCAGCGCTTGGTAAGCCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGc                                                 <  1:314932/226‑1 (MQ=255)
                                                                                                                                                                                                  tACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CCCAAACAACATGCTGCGGC‑GGGGGCGAGGTTTTCGTTGTTTCTGCGA‑AACGGGAAAAATTGGCAATCCTCCCCTTGCAGGGCCcgggaacccccg                                    >  1:114322/1‑230 (MQ=255)
                                                                                                                                                                                                     aGCGGATTTTGCAGGGCGTGGTAAACCATGCGGC‑TGGCGGTTTCCTTTTCACGCAATGCTACGGAGTCTTCGAGGGCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGgcg                                 <  1:544491/239‑1 (MQ=255)
                                                                                                                                                                                                     aGCGGATGTTGCAGGTCGTGGTAAACCATGCGGC‑CGGCGGTGTCTTTTTCCAGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCCGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGgcg                                 <  1:1015032/239‑1 (MQ=255)
                                                                                                                                                                                                       cGGATGTTGCAGCGGTTTGTAAACCATGCGGG‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGcgcg                               <  2:1135655/239‑1 (MQ=255)
                                                                                                                                                                                                         gATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑c                                                                                                                                     <  2:942454/137‑1 (MQ=255)
                                                                                                                                                                                                         gATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑c                                                                                                                                     >  1:942454/1‑137 (MQ=255)
                                                                                                                                                                                                            gTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCCGG‑CACAAACAACATGCCGCGCG‑GGGTGGGGGGCTTTTGTTGTTGTAGCCAGACAGGAAATCATTTGCAGACCCCCCCCGCCAGAGCGGGGagaagcctggcgggggca                           >  2:1019145/1‑222 (MQ=255)
                                                                                                                                                                                                               gggtgcctgtggaaacactgtggc‑gggggGTGTACTTCTCACGCAATGACACGGATGCTTCGAGGGCAAACAGCAGGGCAACAGCGGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGGTTTTCGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCGCCCATGGCAGATCCGTGATACCGCTGGCGCGTGCcagca                       <  2:621475/218‑1 (MQ=255)
                                                                                                                                                                                                               ggcgggtgtgtgtaaCCTTGGGGG‑GGGCGGGTGCTTTTTACCGCAATGCTAGGGAGTCGTCGGGGTCAAACGACAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                       <  2:7844/226‑1 (MQ=255)
                                                                                                                                                                                                               gCAGCGCGTGGTAAACCATGCGGG‑GGGCGGTGTCTTTTTCACGAAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCcagca                       <  1:1260538/239‑1 (MQ=255)
                                                                                                                                                                                                                  gggcgggttAAACCATGGGGC‑GGGGTTGGTGTTTTTCACGCAAGGGTACGGAGTCTTGGAGGTAAAACAACAGGGCATCAGCCGGGGAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAgg                     <  2:52620/230‑1 (MQ=255)
                                                                                                                                                                                                                               taTGCGGG‑GGGCGGTGTCTTTTTCCCGCAATGCTACGGAGGCCTCGAGGGCAAACAGCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGc       <  2:391859/238‑1 (MQ=255)
                                                                                                                                                                                                                                   cggc‑ggGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTc                                                                                                                                                                                                            <  1:117626/41‑1 (MQ=255)
                                                                                                                                                                                                                                   cggc‑ggGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTc                                                                                                                                                                                                            >  2:117626/1‑41 (MQ=255)
                                                                                                                                                                                                                                     gc‑gggcggTGTCTTTTTCACGCAATGCGACGGAATTTTCGAGGTCAAACCACAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCACATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTa  <  1:842859/238‑1 (MQ=255)
                                                                                                                                                                                                                                      gggggcggTGTTTTTTTGCGGCAAGGTTACGGAGTTCACGGTGTCCAACAAGAGGGCAGCACGCGGGTAGAGGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACAGGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTa  <  2:461908/236‑1 (MQ=255)
                                                                                                                                                                                                                                      gggggcggTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCCAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTa  <  1:531470/236‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                        
GGTGCTGCCGGGTTCACGACCACAGGCTTTTTCGATACGCAGGATCTCTTTTTCAATATCCAGAACATCCTGAGCGGTATCGGTTTTCGGCAGACGCACAACGT‑CCGCACCACCGCGAACGACGGCTTCCAGGTCGTTAACACCCCATTCGGAATCCAGCGCGTTGACACGCACAATGGTTTCAATATCGCGATACAGCGGATGTTGCAGCGCGTGGTAAACCATGCGGC‑GGGCGGTGTCTTTTTCACGCAATGCTACGGAGTCTTCGAGGTCAAACATCAGGGCATCAGCCGGGTAGATGAAGGAGTTGCTGACCATCGCGGCATTGGCACCAGG‑CACAAACAACATGCTGCGGC‑GGGTGCGAGTTTTACGTTGTTGCAGCGA‑AGCGGAAATCATTGGCAATCCTCCCATGGCAGAGCCGGGATACCGCTGGCGCGTGCCAGCAGGGCTTCCAGTCGTGCACGTA  >  NC_000913/650095‑650560

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