Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,717,410 A→C 6.1% Q34P (CAA→CCA)  ghrB → glyoxylate/hydroxypyruvate reductase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,717,4100AC 6.1% 200.3 / 4.8 66Q34P (CAA→CCA) ghrBglyoxylate/hydroxypyruvate reductase B
Reads supporting (aligned to +/- strand):  ref base A (33/29);  new base C (2/2);  total (35/31)
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.09e-01

TCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGAT  >  NC_000913/3717269‑3717558
                                                                                                                                             |                                                                                                                                                    
tCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGt                                                                                                                                                <  5:231295/148‑1 (MQ=255)
tCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTc                                                                                                                                               >  5:16967/1‑149 (MQ=255)
            gATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATg                                                                                                                                     >  7:160365/1‑147 (MQ=255)
             aTCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATgca                                                                                                                                   <  8:163921/148‑1 (MQ=255)
                gCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCaa                                                                                                                               >  5:116479/1‑149 (MQ=255)
                         gTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGcc                                                                                                                       >  8:229858/1‑148 (MQ=255)
                         gTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGcc                                                                                                                       >  8:172304/1‑148 (MQ=255)
                         gTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGcc                                                                                                                       >  3:57427/1‑148 (MQ=255)
                            aaTGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAg                                                                                                                   <  5:210370/149‑1 (MQ=255)
                                 agagaAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCCCCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTg                                                                                                                >  5:123307/1‑147 (MQ=255)
                                  gagaAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAg                                                                                                             <  3:37363/149‑1 (MQ=255)
                                   agaAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAgg                                                                                                            <  4:197578/149‑1 (MQ=255)
                                      aGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGttt                                                                                                         <  5:34710/149‑1 (MQ=255)
                                          tGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACCAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTg                                                                                                     >  4:94328/1‑149 (MQ=255)
                                          tGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACCAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTg                                                                                                     >  4:90931/1‑149 (MQ=255)
                                            aaGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTgg                                                                                                    <  6:16967/148‑1 (MQ=255)
                                              ggcgtccgtTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCTCAAACCGTCGAACATAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGtt                                                                                                 <  6:4952/147‑1 (MQ=255)
                                              gCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGtt                                                                                                 <  7:229858/149‑1 (MQ=255)
                                                    cgtTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAAc                                                                                            <  7:17849/148‑1 (MQ=255)
                                                    cgtTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAAc                                                                                            <  8:106622/148‑1 (MQ=255)
                                                         tCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAAt                                                                                      >  1:92212/1‑149 (MQ=255)
                                                           ctctACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATg                                                                                     >  5:237657/1‑148 (MQ=255)
                                                            tctACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGATTGt                                                                                    >  7:198022/1‑148 (MQ=255)
                                                               aCAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAAt                                                                                >  7:248401/1‑149 (MQ=255)
                                                                cAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATg                                                                               >  1:92693/1‑149 (MQ=255)
                                                                 aaaGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGc                                                                              >  1:5962/1‑149 (MQ=255)
                                                                       ttACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCAGCAtt                                                                        >  1:202654/1‑149 (MQ=255)
                                                                           cTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTg                                                                    <  1:136826/149‑1 (MQ=255)
                                                                             gatgatTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTgg                                                                   >  4:55340/1‑148 (MQ=255)
                                                                                  tttACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTg                                                                    <  1:136458/142‑1 (MQ=255)
                                                                                  tttACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTg                                                                    >  2:136458/1‑142 (MQ=255)
                                                                                  tttACTGCAACGCCTGCAAGAGCATGTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAAGGGGAATGTAACTGCCGCATTGCTGGaaaaaa                                                             >  6:186678/1‑149 (MQ=255)
                                                                                         cAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGa                                                       >  2:7181/1‑148 (MQ=255)
                                                                                              ccTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACCAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGc                                                 <  3:94328/149‑1 (MQ=255)
                                                                                              ccTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACCAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGc                                                 <  3:90931/149‑1 (MQ=255)
                                                                                               cTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATg                                                           >  5:151015/1‑138 (MQ=255)
                                                                                               cTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATg                                                           <  6:151015/138‑1 (MQ=255)
                                                                                                tGCAAGAGCATTTCCCCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTg                                              >  2:160359/1‑149 (MQ=255)
                                                                                                 gCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTg                                              >  6:36983/1‑149 (MQ=255)
                                                                                                  cAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGATTGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGt                                               <  8:198022/147‑1 (MQ=255)
                                                                                                   aaGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCCTTGCTGGAAAACATGCCGAAACTGCGTGc                                             >  5:62701/1‑148 (MQ=255)
                                                                                                   aaGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGc                                             >  3:60764/1‑148 (MQ=255)
                                                                                                   aaGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGAAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGca                                            >  6:63583/1‑148 (MQ=255)
                                                                                                     gagCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCca                                           >  4:82193/1‑148 (MQ=255)
                                                                                                     gagCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCca                                           >  5:209773/1‑148 (MQ=255)
                                                                                                     gagCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCca                                           >  4:85548/1‑148 (MQ=255)
                                                                                                        cATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATc                                       <  6:209773/149‑1 (MQ=255)
                                                                                                         aTTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATg                                                           <  7:85938/128‑1 (MQ=255)
                                                                                                         aTTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATg                                                           >  8:85938/1‑128 (MQ=255)
                                                                                                             cACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGa                                  >  8:154184/1‑149 (MQ=255)
                                                                                                               ccGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGAtc                                <  8:248401/149‑1 (MQ=255)
                                                                                                                  ttCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTcc                             >  7:107268/1‑149 (MQ=255)
                                                                                                                     aCCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTc                          >  1:67128/1‑149 (MQ=255)
                                                                                                                         ggTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCt                      <  5:36983/149‑1 (MQ=255)
                                                                                                                          gTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTa                     <  5:63583/149‑1 (MQ=255)
                                                                                                                           tGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTAt                    <  7:197002/149‑1 (MQ=255)
                                                                                                                             gCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTAACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGa                  <  6:153288/149‑1 (MQ=255)
                                                                                                                                aaCCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACaa               <  3:109990/149‑1 (MQ=255)
                                                                                                                                 aCCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAAt              <  7:172304/149‑1 (MQ=255)
                                                                                                                                  ccTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCTACATCAACGATCTCCGTCGGCTATGACAAtt             >  5:88755/1‑149 (MQ=255)
                                                                                                                                   cTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAAttt            <  8:107268/149‑1 (MQ=255)
                                                                                                                                   atCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAAttt            >  7:175678/2‑149 (MQ=255)
                                                                                                                                       gCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGAt        <  8:160365/149‑1 (MQ=255)
                                                                                                                                       gCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGAt        <  2:92693/149‑1 (MQ=255)
                                                                                                                                         ccACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTTATGt      <  3:85548/149‑1 (MQ=255)
                                                                                                                                         ccACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGt      <  3:82193/149‑1 (MQ=255)
                                                                                                                                           acaAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGTTTTAGTGGTTTCAAACGAGAATGTAAATGCCGCATTCCTGAAAAAAATGCCGAAACGGCGTCCCACAACAACGATGTCCTGCGGCTATGACAATTTGGATGTc     >  8:52478/1‑148 (MQ=255)
                                                                                                                                           acaAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTCCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGACGGCTATGACAAATTTGATGTc     >  3:73244/1‑148 (MQ=255)
                                                                                                                                             aaaCCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCCTCGGCTATGACAATTTTGATGTCGAt  <  6:123307/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                    
TCAGGCTAAGGTGATCGCCTTATCAGTGAATGGAGAGAAGCATGAAGCCGTCCGTTATCCTCTACAAAGCCTTACCTGATGATTTACTGCAACGCCTGCAAGAGCATTTCACCGTTCACCAGGTGGCAAACCTCAGCCCACAAACCGTCGAACAAAATGCAGCAATTTTTGCCGAAGCTGAAGGTTTACTGGGTTCAAACGAGAATGTAAATGCCGCATTGCTGGAAAAAATGCCGAAACTGCGTGCCACATCAACGATCTCCGTCGGCTATGACAATTTTGATGTCGAT  >  NC_000913/3717269‑3717558

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