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.5% 249.5 / NA 66L148L (CTG→TTG) ghrBglyoxylate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (33/32);  minor base A (1/0);  total (34/32)
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.98e-01

AACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAA  >  NC_000913/3717522‑3717983
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
aaCGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGt                                                                                                                                                                                                                                <  1:200002/240‑1 (MQ=255)
     tctcCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGa                                                                                                                                                                                                                           <  1:687779/240‑1 (MQ=255)
               cTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGAt                                                                                                                                                                                                                 <  1:708696/240‑1 (MQ=255)
                tATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGAt                                                                                                                                                                                                                 <  1:452911/239‑1 (MQ=255)
                 aTGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGa                                                                                                                                                                                                                      <  2:375191/233‑1 (MQ=255)
                      aaTTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCAt                                                                                                                                                                                                           <  1:374202/239‑1 (MQ=255)
                            gATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGc                                                                                                                                                                                                             <  1:345174/231‑1 (MQ=255)
                                      cgcTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTg                                                                                                                                                                                          <  1:577326/240‑1 (MQ=255)
                                        cTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTgcg                                                                                                                                                                                        >  2:670740/1‑240 (MQ=255)
                                                 cGATAAATTCTGCTGATGCACACGCCTACCGTATTAACAGAAAACGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGTCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGc                                                                                                                                                                               <  1:307051/240‑1 (MQ=255)
                                                 cGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGGTGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTgg                                                                                                                                                                                >  1:109396/1‑239 (MQ=255)
                                                 cGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTgg                                                                                                                                                                                >  2:397269/1‑239 (MQ=255)
                                                  gAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGc                                                                                                                                                                               <  2:123614/239‑1 (MQ=255)
                                                           tgctgATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGc                                                                                                                                                                     <  1:79273/240‑1 (MQ=255)
                                                                          cAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGcac                                                                                                                                                                                                >  2:487307/1‑198 (MQ=255)
                                                                                  ttAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGa                                                                                                                                                                                                                      <  2:721164/168‑1 (MQ=255)
                                                                                  ttAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGcgccgc                                                                                                                                              >  1:56407/1‑240 (MQ=255)
                                                                                  ttAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGcgccgc                                                                                                                                              <  1:615952/240‑1 (MQ=255)
                                                                                  ttAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGcgccgc                                                                                                                                              >  1:538104/1‑240 (MQ=255)
                                                                                        gAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGccacca                                                                                                                                         >  1:125963/1‑239 (MQ=255)
                                                                                          aaCCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATgg                                                                                                                                                                                                         >  1:279372/1‑173 (MQ=255)
                                                                                           aCCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCctc                                                                                                                                                             >  2:278124/1‑216 (MQ=255)
                                                                                               tCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAaagaa                                                                                                                                  >  1:446095/1‑239 (MQ=255)
                                                                                                          cTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGccc                                                                                                                                                                   <  1:477878/195‑1 (MQ=255)
                                                                                                                   cTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGcc                                                                                                             <  1:153191/240‑1 (MQ=255)
                                                                                                                   cTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGcc                                                                                                             <  1:54908/240‑1 (MQ=255)
                                                                                                                            tCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGc                                                                                                    <  1:419342/240‑1 (MQ=255)
                                                                                                                                 cGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACAGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGAttt                                                                                               >  1:192510/1‑240 (MQ=255)
                                                                                                                                   cTCGTCGGGTTGTGGAGGTGGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGATTAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGccac                                                                                                                                           <  2:88128/194‑1 (MQ=255)
                                                                                                                                   cTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTgg                                                                                             >  1:587482/1‑240 (MQ=255)
                                                                                                                                   cTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTg                                                                                              >  2:648764/1‑239 (MQ=255)
                                                                                                                                     cgtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGa                                                                                            <  1:730321/239‑1 (MQ=255)
                                                                                                                                      gtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATa                                                                                          >  2:98902/1‑240 (MQ=255)
                                                                                                                                      gtcgGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATa                                                                                          >  1:700703/1‑239 (MQ=255)
                                                                                                                                         gggTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATAct                                                                                        <  1:231870/239‑1 (MQ=255)
                                                                                                                                          ggTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAaaga                                                                                                                                   <  2:660897/195‑1 (MQ=255)
                                                                                                                                                     tAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAaga                                                                             >  2:651793/1‑238 (MQ=255)
                                                                                                                                                     tAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGt                                                                           <  1:578169/240‑1 (MQ=255)
                                                                                                                                                      aGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGAAGTTCACCATAAAAGAAAGGGCATTGTCGGGAGGGGACGGATCGGCATGGGGCTGGCACAACGTgc                                                                                                                                                                                         >  2:260695/1‑129 (MQ=255)
                                                                                                                                                              gggTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCg                                                                  <  1:651796/240‑1 (MQ=255)
                                                                                                                                                              gggTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCg                                                                  <  2:366878/240‑1 (MQ=255)
                                                                                                                                                                  aaaaGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATa                                                                                          <  1:580834/212‑1 (MQ=255)
                                                                                                                                                                                 gACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGtt                                               <  2:538103/240‑1 (MQ=255)
                                                                                                                                                                                   ccGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTaa                                             <  2:647578/240‑1 (MQ=255)
                                                                                                                                                                                      cgAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACtg                                          >  2:701709/1‑240 (MQ=255)
                                                                                                                                                                                             aGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATgaga                                    <  2:56407/239‑1 (MQ=255)
                                                                                                                                                                                             aGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATGGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCATTAACTGATgaga                                    <  2:561059/239‑1 (MQ=255)
                                                                                                                                                                                               gCCCGGACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGc                                 <  1:241517/240‑1 (MQ=255)
                                                                                                                                                                                                    gACTGGTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAaga                                                                             >  1:463427/1‑191 (MQ=255)
                                                                                                                                                                                                         gTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGGGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGtt                       >  1:305603/1‑240 (MQ=255)
                                                                                                                                                                                                         gTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCatcatc                           >  1:377688/1‑236 (MQ=255)
                                                                                                                                                                                                         gTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGtt                       >  1:151700/1‑240 (MQ=255)
                                                                                                                                                                                                         gTACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGtt                       >  1:432783/1‑240 (MQ=255)
                                                                                                                                                                                                          tACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATa                                                                                          >  1:169033/1‑172 (MQ=255)
                                                                                                                                                                                                          tACGGCACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATa                                                                                          <  2:169033/172‑1 (MQ=255)
                                                                                                                                                                                                               cACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGtt                       >  2:132270/1‑234 (MQ=255)
                                                                                                                                                                                                               cACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGtt                       <  1:132270/234‑1 (MQ=255)
                                                                                                                                                                                                               cACTGACGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGgcg                  >  1:271230/1‑239 (MQ=255)
                                                                                                                                                                                                                    aCGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGca                                >  1:536585/1‑220 (MQ=255)
                                                                                                                                                                                                                    aCGTTCACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGaa             >  1:3408/1‑239 (MQ=255)
                                                                                                                                                                                                                         cACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGgc                   >  2:571489/1‑228 (MQ=255)
                                                                                                                                                                                                                         cACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTc       >  1:496223/1‑240 (MQ=255)
                                                                                                                                                                                                                         cACCATAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTc       >  1:661512/1‑240 (MQ=255)
                                                                                                                                                                                                                              tAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGcc                                                  >  2:182510/1‑192 (MQ=255)
                                                                                                                                                                                                                              tAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCTaa  >  1:734523/1‑240 (MQ=255)
                                                                                                                                                                                                                              tAAAACATTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCaa  >  1:226472/1‑240 (MQ=255)
                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                        
AACGATCTCCGTCGGCTATGACAATTTTGATGTCGATGCGCTTACCGCCCGAAAAATTCTGCTGATGCACACGCCAACCGTATTAACAGAAACCGTCGCCGATACGCTGATGGCGCTGGTGTTGTCTACCGCTCGTCGGGTTGTGGAGGTAGCAGAACGGGTAAAAGCAGGCGAATGGACCGCGAGCATAGGCCCGGACTGGTACGGCACTGACGTTCACCATAAAACACTGGGCATTGTCGGGATGGGACGGATCGGCATGGCGCTGGCACAACGTGCGCACTTTGGCTTCAACATGCCCATCCTCTATAACGCGCGCCGCCACCATAAAGAAGCAGAAGAACGCTTCAACGCCCGCTACTGCGATTTGGATACTCTGTTACAAGAGTCAGATTTCGTTTGCCTGATCCTGCCGTTAACTGATGAGACGCATCATCTGTTTGGCGCAGAACAATTCGCCAA  >  NC_000913/3717522‑3717983

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