Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,717,751 C→T L148L (CTG→TTG)  ghrB → glyoxylate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,717,7510CT98.4% 230.1 / ‑6.8 61L148L (CTG→TTG) ghrBglyoxylate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/28);  minor base G (0/1);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

ACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATT  >  NC_000913/3717517‑3717977
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
acaTCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGc                                                                                                                                                                                                                               <  1:688209/240‑1 (MQ=255)
acaTCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGc                                                                                                                                                                                                                               <  1:196429/240‑1 (MQ=255)
          tctcCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGa                                                                                                                                                                                                                     <  1:109515/240‑1 (MQ=255)
          tctcCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGa                                                                                                                                                                                                                     <  1:8069/240‑1 (MQ=255)
                          cAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGc                                                                                                                                                                                                                               <  2:81618/214‑1 (MQ=255)
                           aaTTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACTGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCAt                                                                                                                                                                                                     <  2:303512/239‑1 (MQ=255)
                              tttGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGgcg                                                                                                                                                                                                 <  1:56344/240‑1 (MQ=255)
                                       gATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACaa                                                                                                                                                                                        >  1:135109/1‑240 (MQ=255)
                                               tACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTgcgc                                                                                                                                                                                 >  1:731325/1‑239 (MQ=255)
                                                 ccGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACt                                                                                                                                                                              <  1:597625/240‑1 (MQ=255)
                                                                tgctgATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGc                                                                                                                                                               <  1:701134/240‑1 (MQ=255)
                                                                  ctgATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGccc                                                                                                                                                             <  1:589218/240‑1 (MQ=255)
                                                                   tgATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGccc                                                                                                                                                             <  2:502435/239‑1 (MQ=255)
                                                                    gATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCAt                                                                                                                                                           <  2:55579/240‑1 (MQ=255)
                                                                     aTGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATc                                                                                                                                                          >  1:554239/1‑240 (MQ=255)
                                                                                      aTTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCAt                                                                                                                                                                                                     <  2:64450/180‑1 (MQ=255)
                                                                                       ttAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGcgccgc                                                                                                                                        <  2:598238/240‑1 (MQ=255)
                                                                                       ttAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGcgccgc                                                                                                                                        >  1:661392/1‑240 (MQ=255)
                                                                                          aCAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGccac                                                                                                                                     <  1:116744/240‑1 (MQ=255)
                                                                                             gAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACgcgcgc                                                                                                                                           >  2:320715/1‑231 (MQ=255)
                                                                                             gAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACgcgcgc                                                                                                                                           <  1:320720/231‑1 (MQ=255)
                                                                                               aaCCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATaa                                                                                                                                >  1:555867/1‑240 (MQ=255)
                                                                                                 ccGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATaaa                                                                                                                               >  1:520678/1‑239 (MQ=255)
                                                                                                    tCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAaagaa                                                                                                                            >  1:475466/1‑239 (MQ=255)
                                                                                                      gCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATaa                                                                                                                                                  >  2:539685/1‑215 (MQ=255)
                                                                                                      gCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGc                                                                                                                          >  2:562468/1‑239 (MQ=255)
                                                                                                                        cTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGcc                                                                                                       >  1:98895/1‑240 (MQ=255)
                                                                                                                                   tACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGa                                                                                            <  1:715272/240‑1 (MQ=255)
                                                                                                                                     ccGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTAc                                                                                                 >  2:463749/1‑233 (MQ=255)
                                                                                                                                      cGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGAttt                                                                                         >  2:117431/1‑240 (MQ=255)
                                                                                                                                        cTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTgg                                                                                       >  1:399632/1‑240 (MQ=255)
                                                                                                                                           gtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGccacc                                                                                                                                    >  2:245610/1‑192 (MQ=255)
                                                                                                                                           gtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGAt                                                                                     >  2:694687/1‑239 (MQ=255)
                                                                                                                                           gtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGAt                                                                                     <  2:484909/239‑1 (MQ=255)
                                                                                                                                             cgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATAc                                                                                   <  1:48292/239‑1 (MQ=255)
                                                                                                                                             cgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATAc                                                                                   <  2:487266/239‑1 (MQ=255)
                                                                                                                                             cgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATAc                                                                                   <  2:731312/239‑1 (MQ=255)
                                                                                                                                              gggTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATAct                                                                                  <  2:551470/239‑1 (MQ=255)
                                                                                                                                                         gTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGcgcc                                                                                                                                          >  2:453636/1‑172 (MQ=255)
                                                                                                                                                         gTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATAc                                                                                   >  1:624211/1‑227 (MQ=255)
                                                                                                                                                         gTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAagag                                                                      >  1:213525/1‑240 (MQ=255)
                                                                                                                                                              aGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGa                                                                 >  1:514931/1‑240 (MQ=255)
                                                                                                                                                                     gTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTc                                                                    >  2:452949/1‑230 (MQ=255)
                                                                                                                                                                      tAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGa                                                                                            <  2:549749/205‑1 (MQ=255)
                                                                                                                                                                                  aaTGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACAGTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATc                                                 <  1:445980/236‑1 (MQ=255)
                                                                                                                                                                                                     cccGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGc                           <  2:369013/239‑1 (MQ=255)
                                                                                                                                                                                                      ccGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGcaca                                                                                                                                                                                         <  1:532556/80‑1 (MQ=255)
                                                                                                                                                                                                      ccGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGcaca                                                                                                                                                                                         >  2:532547/1‑80 (MQ=255)
                                                                                                                                                                                                        ggACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGcatc                        >  1:707972/1‑239 (MQ=255)
                                                                                                                                                                                                         gACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTAc                                                                           <  2:236557/187‑1 (MQ=255)
                                                                                                                                                                                                         gACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTAc                                                                           >  1:236560/1‑187 (MQ=255)
                                                                                                                                                                                                         gACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACtgatga                                >  2:647303/1‑230 (MQ=255)
                                                                                                                                                                                                              gTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATActc                                                                                 >  2:732274/1‑176 (MQ=255)
                                                                                                                                                                                                              gTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGAc                             >  2:414123/1‑228 (MQ=255)
                                                                                                                                                                                                              gTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGtt                 <  1:722942/240‑1 (MQ=255)
                                                                                                                                                                                                                    cACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAaga                                                                       >  1:682302/1‑180 (MQ=255)
                                                                                                                                                                                                                    cACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGgcg            >  2:589450/1‑239 (MQ=255)
                                                                                                                                                                                                                         aCGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGtt                 <  2:211659/229‑1 (MQ=255)
                                                                                                                                                                                                                         aCGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGaa       >  2:249408/1‑239 (MQ=255)
                                                                                                                                                                                                                              cACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGAt                                                                >  2:319977/1‑177 (MQ=255)
                                                                                                                                                                                                                              cACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAAtt  >  1:682448/1‑239 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
ACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATT  >  NC_000913/3717517‑3717977

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