Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,403,735 T→C 7.3% L85P (CTC→CCC)  acuI → putative acryloyl‑CoA reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,403,7350TC 7.3% 223.4 / 2.2 78L85P (CTC→CCC) acuIputative acryloyl‑CoA reductase
Reads supporting (aligned to +/- strand):  ref base T (32/39);  new base C (3/4);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGACGGGGAGATTGTCGTGACGGGTGCCAGTGGTGGCGTCGGCAGTACCGCCGTG  >  NC_000913/3403506‑3403973
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
cagGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGGGTTACTCCctggctg                                                                                                                                                                                                                                       >  1:501199/1‑239 (MQ=255)
         tgcaggtgtgtataaaggaagtctctGGAGGAAAGTCGCTCGCCGGAGGGCGATGTCACGGTCGATGTTCCCTGGGCGGTCCAGGACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGt                                                                                                                                                                                                                                 <  2:238634/214‑1 (MQ=255)
          aaaCTCTCCCCTCAGTACCGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATTTCACGGTCGATGTTCACTGCGCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCCTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGt                                                                                                                                                                                                                             <  1:906875/239‑1 (MQ=255)
             cgctcGCATCCGTACAGACTCTGGGCGAAAGTCGCCTGCCGGATGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGaa                                                                                                                                                                                                                          <  1:283272/237‑1 (MQ=255)
             cctccgtaTCAGTACAGGCTTGGGACGAAATTCGCCTGCTGGAGGGAGATGGCACGGTCGATGTTCAATGGTCGAGCCTGAACTATAAA‑GATGCGCTGGTGATTCCCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGaa                                                                                                                                                                                                                          <  2:448002/232‑1 (MQ=255)
                  gCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAAcc                                                                                                                                                                                                                      >  1:853841/1‑238 (MQ=255)
                  ccctccgtccAGACTTCGGAAGAAAGGCGCCTGCCGGGGGGCGATGTCACGGTCTATGTTGACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCa                                                                                                                                                                                                                     <  2:772423/230‑1 (MQ=255)
                        ttttgggttttgtgggAAAGTCGCCTGGCGGAGGGCGGTTTTACGGTCGGTGGTCACTGGGCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTgggg                                                                                                                                                                                                               <  2:347170/224‑1 (MQ=255)
                         tACAAACTCTGGAAGGAAGTCCGCTGCCGGGGGGCGATGGCACGGTCGATGTTCACGGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTGCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCGCCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTgggg                                                                                                                                                                                                               <  2:501199/238‑1 (MQ=255)
                            cgACTTTGGAAGAAAGTCGCCTGCGGGGGGGGGCTGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAACCATCCGTAATTTGCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTggg                                                                                                                                                                                                           <  2:369426/238‑1 (MQ=255)
                              aCTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTCCCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATTCCGGTCAGGAGGTGTTACTCACTGCCTGGGGGGTTGGGGAAAACCACTGGGGTgg                                                                                                                                                                                                            >  1:122924/1‑236 (MQ=255)
                              aCTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGCACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCCGGCGGCGTTACTCACTGGCTGGGGCGTTGGTGAAAACCAATTGGGTgg                                                                                                                                                                                                            >  1:886016/1‑236 (MQ=255)
                              aCTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGGGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTgg                                                                                                                                                                                                            >  1:750979/1‑236 (MQ=255)
                               ccttGGACGAAAGTGGCCTGCCGGAGGGCGATGGCACGGTCGATTTTCACTGGTCGGGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAAGTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTggg                                                                                                                                                                                                           <  2:853841/233‑1 (MQ=255)
                                tctGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGaa                                                                                                                                                                                                                          <  1:846817/220‑1 (MQ=255)
                                  tGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGGGTTACCCACTGGCTGGGGCGTTGGTGAAAACCCCTGGGGGGGGCggggg                                                                                                                                                                                                     >  2:625968/1‑237 (MQ=255)
                                  cggACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGGAAACCACTGGGGGGGGCTGGCg                                                                                                                                                                                                     >  1:903197/2‑239 (MQ=255)
                                   ggACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCGGGTGGTGTTACCCACGGGATGGGGCGGTGTGGAAAACCACTGGGGGGGGATGGCgg                                                                                                                                                                                                    >  1:379046/1‑239 (MQ=255)
                                   ggACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGGGGGGTTACTCACTGGCTGGGGCGTTGGGGAAAACCACGggggggggggggggg                                                                                                                                                                                                    >  2:718957/1‑232 (MQ=255)
                                   ggACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTc                                                                                                                                                                                                                                               <  2:94842/196‑1 (MQ=255)
                                   ggACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTc                                                                                                                                                                                                                                               >  1:94842/1‑196 (MQ=255)
                                   ggACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGGGTTACTCCCTGGCTGGGGGGGTGGGGAAAACCCCTGGGGGGGGGcggcgg                                                                                                                                                                                                    >  2:549312/1‑239 (MQ=255)
                                          aGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTTCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATCCCGGTAGGGGGGTGTTACTCACTGGCTGGGGCGGTGGTGAAACACCCTgggggtgtctctttaccacatc                                                                                                                                                                                             >  2:846817/1‑221 (MQ=255)
                                          aGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCCGGAGGTGTTACTCACTGACTGGGGCGTTGGCGAAAACCACTGGGGTGGGCTGCCCGAGCCGgc                                                                                                                                                                                             >  1:84147/1‑239 (MQ=255)
                                          aGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCAGGGGCGTTGTTGACAAc                                                                                                                                                                                                                       >  1:238634/1‑213 (MQ=255)
                                                    cGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGaa                                                                                                                                                                                    <  1:56617/238‑1 (MQ=255)
                                                             aTGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGACCCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGCGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGGGActgg                                                                                                                                                                          >  1:206704/1‑239 (MQ=255)
                                                                         aTGTTCACTGGCCGAGCCTTAACTATAAA‑GATGCGCTGGCGATTGCCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCAt                                                                                                                                                               <  1:411532/238‑1 (MQ=255)
                                                                              cACGGGGCGAGCCTGGACTATAAAGGGTACG‑TGGCGAGTCCCGGTAGGGGAAAAATCATCCGTAATTTTCCGATGATTCCTTGGATCGATGTTGGCGGAACTG‑TACGCTCCAGCGAAGATCCGCGTTTTCATGCCGGTCGGGAGGGGTTACTCACTGGCTGGGGCGTTGGGGAAACCCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCa                                                                                                                                                         <  2:404622/239‑1 (MQ=255)
                                                                               cgtggtcggacctgatatgaaaa‑TATGCGCGGGCGTTTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGTGTTTGCCGGAACTT‑GACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCa                                                                                                                                                         <  2:886016/218‑1 (MQ=255)
                                                                                    tCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCa                                                                                                                                                                                                                     >  2:622296/1‑173 (MQ=255)
                                                                                    tCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCa                                                                                                                                                                                                                     <  1:622296/173‑1 (MQ=255)
                                                                                    tCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGGGACGGCGGGGGGGGAAGGGGACGGGCTGGGTGCCATCCCCCAAggggg                                                                                                                                                   >  2:176759/1‑237 (MQ=255)
                                                                                        cgcTGAACTTAAAA‑GGTGCGCTGGCAATTACCGGTAAGGGAAAAATCCTCCGTAATTTTCCGATGATTCCTGGGATCGATTGTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGAc                                                                                                                                               <  1:23333/237‑1 (MQ=255)
                                                                                                 ttaAA‑GATGCGCTTGGGATTCCCGGTAAGGG‑AAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAg                                                                                                                                     <  1:710953/238‑1 (MQ=255)
                                                                                                  taAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCActggctg                                                                                                                                                                                                                                       >  2:439038/1‑141 (MQ=255)
                                                                                                  taAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCActggctg                                                                                                                                                                                                                                       <  1:439038/141‑1 (MQ=255)
                                                                                                  taAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGGGAGTGAAAGGTGACTGGCCGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAg                                                                                                                                     >  1:400709/1‑239 (MQ=255)
                                                                                                  taAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAg                                                                                                                                     <  1:608410/239‑1 (MQ=255)
                                                                                                     aggttgtgctgtggtttcccGGTAAGGGAAAAATACTCCGTAATTTCCCGATGATGCCTGGGATCGATTTTGCCGGACCTG‑GACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCa                                                                                                                                   <  2:122924/221‑1 (MQ=255)
                                                                                                        aTGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGGGAAAGGTGACTGGCTGGTTGACATGCCGCAAGGGCTGGACGAGGGCAAAGCAATg                                                                                                                                >  1:139755/1‑239 (MQ=255)
                                                                                                         agcactggggtgtACCGGTAAGGGAAAAATCACCCGTAATTTTCCGATGATTCCGGGGATCGATTTTGCCGGAACTG‑GACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTGACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATg                                                                                                                                <  1:568399/226‑1 (MQ=255)
                                                                                                           cgcTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGAtt                                                                                                                             >  1:654604/1‑239 (MQ=255)
                                                                                                           cgcTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGGGAAAGGGGACTGGCTGGTTGCCCTGCCGCAAGGGCTGGAGCCGCGGAAAGCAAGGAtt                                                                                                                             >  1:47753/1‑239 (MQ=255)
                                                                                                              tGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTAt                                                                                                                           <  1:549312/238‑1 (MQ=255)
                                                                                                                 gggttACCGGTAAGGGAAAAATTCTCCGTAATTTTCCGTTTATTCCTGGGATCGATTTTGCCGGAACGT‑TACGCCCCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCgg                                                                                                                        <  2:287868/235‑1 (MQ=255)
                                                                                                                 ctaTTACCGGTAAGGGAAAAATCATCCGTAAGTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGAACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCgg                                                                                                                        <  1:176759/236‑1 (MQ=255)
                                                                                                                 ccaTTACCGGTAAGGGAAAAATCATCCGTAATTTCCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCgg                                                                                                                        <  1:291490/236‑1 (MQ=255)
                                                                                                                       ccGGAAAGGGAAAAATCCCCCGCTATTTCCCGATTATTCCTGGGGTCGTTTTGGCCGGACCTG‑TACGCCCCGGCGAAGATCCGCGTTTTCATGCCGGTCGGGAGGTTTTACTCACTGGCTGGGGCGTTGTTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGcc                                                                                                                 <  2:654604/239‑1 (MQ=255)
                                                                                                                          gTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGAAGCGGGTAAGGCAATGATAATCGGTAATGCCGGt                                                                                                              >  1:681441/1‑239 (MQ=255)
                                                                                                                                gAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCCCCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTa                                                                                                          <  1:625968/237‑1 (MQ=255)
                                                                                                                                aaaaaaTCATCCGTAATTTTCCGATGATTCCTGGTATCGATTTTGCCGGAACTG‑TTCGGACCAGCGAAGATCCGCGTTTTCAGGCCGGTCAGGAGGTGTTACCCACTGGCTGGGGCGTTGGTGAAACCCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTa                                                                                                          <  1:915018/236‑1 (MQ=255)
                                                                                                                                   aaaTCCTCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGTTCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGcc                                                                                                     <  1:718957/239‑1 (MQ=255)
                                                                                                                                   aaaTCATCCTTATTTTTCCGTTTTTTCCTGGGTTCGTTTTTGCCGAGACGG‑TCCGCCCCGGGGAAGTTCCGGGTTTTCATGCCGGTCAGGAGTTGTTACTCCCTGGTGGGGGCGTTGGTGAAAACCCTTGGGGGGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCAGGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGcc                                                                                                     <  2:47753/239‑1 (MQ=255)
                                                                                                                                    aaTCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCa                                                                                                                                                         <  2:786731/186‑1 (MQ=255)
                                                                                                                                    aaTCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCa                                                                                                                                                         >  1:786731/1‑186 (MQ=255)
                                                                                                                                         tcCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGGGGGAAAGGGAACCGGCGGGTTGCCATGCCCCCAGGGGTGGGCCCGCGCAAAAGAATGATTTTTGGTCATGCCGGGTTTTACGCCCTgggg                                                                                               >  2:896904/1‑236 (MQ=255)
                                                                                                                                                  ttcccgtttgttcctgggctggctttGGCGGGACTGTTACGCACCAGCGAAGATCCGGTTTTTAATGCCGGTCAGGGGGTGTTACTCACTGGCTGGGGGGTTGGGGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGGGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGAAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGAt                                                                                       <  2:400709/216‑1 (MQ=255)
                                                                                                                                                            tCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCCTGCTGTTCGGGATGGCGCTGGAAg                                                                            >  1:15064/1‑239 (MQ=255)
                                                                                                                                                                              ggggagtG‑TCCGCACAGGGGAATATCCGGGTTTTCAGGCCGGTTGGGGGGTGTTACCCCCGGGGTGGGGCGTGGGTGAAACCCCCGGGGGTGGGCTGGCGGGGCAGGCGCGGGTGAAAGGTGCCGGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCg                                                           <  2:522080/232‑1 (MQ=255)
                                                                                                                                                                              cGGAACTG‑TCCGCCCCAGGGAAGTTCCGCTTTTTCATGCCGGTCAGGAGGTGTTACTCGCTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGc                                                          <  1:135795/239‑1 (MQ=255)
                                                                                                                                                                                       tACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGGGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGACggg                                                  >  1:522080/1‑239 (MQ=255)
                                                                                                                                                                                        aCGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACCCACTGGCTGGGGCGTTGGTGAAGACCACTGGGGTGGGCTgg                                                                                                                                                                                                       >  2:351026/1‑89 (MQ=255)
                                                                                                                                                                                        aCGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACCCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTgg                                                                                                                                                                                                       <  1:351026/89‑1 (MQ=255)
                                                                                                                                                                                            ccaagggAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTCCTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGACGGGgaga                                              <  1:531293/232‑1 (MQ=255)
                                                                                                                                                                                                      ttCCGGGTTTGCAGGCCGGTCAGGGGGTTTTACCCCCTGGCGGGGGGGTTGGGGAAACCCGATGGGGTGGGCTGGCGGAGCAGGCGCGAGTTAAAGGTGCCTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGACGGGGAGATTGTCGTGACg                                   <  1:135527/238‑1 (MQ=255)
                                                                                                                                                                                                           cgTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGACGGGGGGATTGTCGTGACGGGGGc                              >  1:395625/1‑239 (MQ=255)
                                                                                                                                                                                                                                 tgtTACTCACTGGCTGGGGCGTTGGTGAAAACCACTgtgg                                                                                                                                                                                                               <  1:847321/40‑1 (MQ=255)
                                                                                                                                                                                                                                 tgtTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTa                                                                                                          <  1:705221/141‑1 (MQ=255)
                                                                                                                                                                                                                                  gtTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGAGGGGGAGATTGTCGTGACGGGTGCCAGGGGTGGCGTCGGCCGTAcc        >  1:882867/1‑238 (MQ=255)
                                                                                                                                                                                                                                   ttACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGAc                                                                                                                                               <  2:452368/102‑1 (MQ=255)
                                                                                                                                                                                                                                   ttACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGAc                                                                                                                                               >  1:452368/1‑102 (MQ=255)
                                                                                                                                                                                                                                   ttACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGAtt                                                                                                                             <  2:656077/120‑1 (MQ=255)
                                                                                                                                                                                                                                   ttACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGAtt                                                                                                                             >  1:656077/1‑120 (MQ=255)
                                                                                                                                                                                                                                   ttACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGcc                                                                                                                 <  1:848035/132‑1 (MQ=255)
                                                                                                                                                                                                                                   ttACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGAAATGGTTATCGGTACTGcc                                                                                                                 >  2:848035/1‑132 (MQ=255)
                                                                                                                                                                                                                                     aCTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGACGGGGAGATTGTCCGGACGGGTGCCCGGGGGGGGGTCGTCAGGAccggcg    >  1:687793/1‑239 (MQ=255)
                                                                                                                                                                                                                                     aCTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGGGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGGAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGGGATGGCGCTGGAAGATGCCGGGGTTCGCCCGCCGGCCGGGGGGGTTGGCGGGACGGGGGCCAGggggggggcgggcgggacgcccg    >  2:813132/1‑222 (MQ=255)
                                                                                                                                                                                                                                       tCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGACGGGGGGATTGTCGTGACGGGTGCCAGTGGGGGCGTCGGCAGTACCGCCGTg  >  1:909480/1‑239 (MQ=255)
                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                              
CAGGACGGCAAAACTCTCGCATCAGTACAGACTCTGGACGAAAGTCGCCTGCCGGAGGGCGATGTCACGGTCGATGTTCACTGGTCGAGCCTGAACTATAAA‑GATGCGCTGGCGATTACCGGTAAGGGAAAAATCATCCGTAATTTTCCGATGATTCCTGGGATCGATTTTGCCGGAACTG‑TACGCACCAGCGAAGATCCGCGTTTTCATGCCGGTCAGGAGGTGTTACTCACTGGCTGGGGCGTTGGTGAAAACCACTGGGGTGGGCTGGCGGAGCAGGCGCGAGTGAAAGGTGACTGGCTGGTTGCCATGCCGCAAGGGCTGGACGCGCGTAAAGCAATGATTATCGGTACTGCCGGTTTTACCGCCATGCTGTGTGTGATGGCGCTGGAAGATGCCGGTGTTCGCCCGCAGGACGGGGAGATTGTCGTGACGGGTGCCAGTGGTGGCGTCGGCAGTACCGCCGTG  >  NC_000913/3403506‑3403973

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