Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,010,797 A→T 5.2% Q257L (CAG→CTG)  hyuA → D‑stereospecific phenylhydantoinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,010,7970AT 5.2% 244.7 / 3.1 77Q257L (CAG→CTG) hyuAD‑stereospecific phenylhydantoinase
Reads supporting (aligned to +/- strand):  ref base A (34/39);  new base T (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.57e-01

TGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACTGTGGTGTGGCATCAGCGA  >  NC_000913/3010653‑3010938
                                                                                                                                                |                                                                                                                                             
tGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCAccagc                                                                                                                                             >  2:9871/1‑147 (MQ=255)
           ctGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGtt                                                                                                                                  >  6:226170/1‑147 (MQ=255)
           ctGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGtt                                                                                                                                  >  7:158896/1‑147 (MQ=255)
           ctGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCAGTGCGAATCACCAGCCAGTCTGGGtt                                                                                                                                  >  5:109443/1‑147 (MQ=255)
              taaTGCGAAGCGGAAGCGATCGCCAGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAAGGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGaa                                                                                                                               <  3:196367/146‑1 (MQ=255)
              gAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGTCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTg                                                                                                                                 <  5:11028/145‑1 (MQ=255)
                 tGCGAAGCGGAAGCCATCGCCCGCATGAGTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCt                                                                                                                           <  5:120924/148‑1 (MQ=255)
                      aGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCca                                                                                                                      <  6:173064/148‑1 (MQ=255)
                        cGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCcaca                                                                                                                    <  7:336912/148‑1 (MQ=255)
                                 tCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAAcc                                                                                                                            <  6:158095/131‑1 (MQ=255)
                                 tCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAAcc                                                                                                                            >  5:158095/1‑131 (MQ=255)
                                 tCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCt                                                                                                           <  6:31888/148‑1 (MQ=255)
                                  cGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTg                                                                                                          >  1:356866/1‑148 (MQ=255)
                                  cGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTg                                                                                                          >  4:40192/1‑148 (MQ=255)
                                  cGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTg                                                                                                          >  4:123749/1‑148 (MQ=255)
                                           gATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTATGTCTGTATTATCTGCTTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCAACAAAATCTCCTGTTGGACGaa                                                                                                 >  6:59625/1‑148 (MQ=255)
                                            aTTAACCTGGCACAAATTGCCGTTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTTAGTATCACCAGCAAGTCTGGGTTGAAACCTGCCGACAATATCTCCTGTTGGACGAAc                                                                                                >  8:231936/1‑148 (MQ=255)
                                              tAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACg                                                                                               >  5:347260/1‑147 (MQ=255)
                                                         aaaTTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATa                                                                                    >  2:212005/1‑147 (MQ=255)
                                                          aaTTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCTGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGtt                                                                                          <  8:389603/140‑1 (MQ=255)
                                                          aaTTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCTGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGtt                                                                                          >  7:383997/1‑140 (MQ=255)
                                                          aaTTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCTGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGtt                                                                                          >  7:389603/1‑140 (MQ=255)
                                                          aaTTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCTGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGtt                                                                                          <  8:383997/140‑1 (MQ=255)
                                                          aaTTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACa                                                                                  <  3:40192/148‑1 (MQ=255)
                                                          aaTTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACa                                                                                  <  8:111482/148‑1 (MQ=255)
                                                              gCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACagaa                                                                               <  1:9871/147‑1 (MQ=255)
                                                               ccGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACagaag                                                                              <  3:123749/147‑1 (MQ=255)
                                                               ccGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAAGCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACagaag                                                                              <  6:40843/147‑1 (MQ=255)
                                                                 ggTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCt                                                                                                                                       >  2:360302/1‑88 (MQ=255)
                                                                 ggTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCt                                                                                                                                       <  1:360302/88‑1 (MQ=255)
                                                                      cGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCAt                                                                       >  5:105820/1‑147 (MQ=255)
                                                                       gCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGa                                                                     <  7:69457/148‑1 (MQ=255)
                                                                       gCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCGGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGa                                                                     <  6:347260/148‑1 (MQ=255)
                                                                             ctctATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAACCCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCa                                                               >  6:103484/1‑148 (MQ=255)
                                                                              gctATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCAt                                                              <  1:238480/147‑1 (MQ=255)
                                                                               ctATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCAtt                                                             >  1:29982/1‑148 (MQ=255)
                                                                                 atatCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTc                                                                <  2:68371/143‑1 (MQ=255)
                                                                                 atatCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTc                                                                >  1:68371/1‑143 (MQ=255)
                                                                                  tatCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCAtt                                                             <  7:369078/145‑1 (MQ=255)
                                                                                  tatCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCtt                                                          <  8:188024/148‑1 (MQ=255)
                                                                                   atCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGaaa                                                                                                                              >  2:81813/1‑79 (MQ=255)
                                                                                   atCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGaaa                                                                                                                              <  1:81813/79‑1 (MQ=255)
                                                                                    tCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAg                                                        >  6:390246/1‑148 (MQ=255)
                                                                                     cGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGc                                                       >  4:99415/1‑148 (MQ=255)
                                                                                           ccTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCAccagccag                                                                                                                                          <  1:128687/59‑1 (MQ=255)
                                                                                           ccTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCAccagccag                                                                                                                                          >  2:128687/1‑59 (MQ=255)
                                                                                           ccTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCg                                                 >  6:152281/1‑148 (MQ=255)
                                                                                           acTGTCAAACGGCTTAGGTCTGGATTATCTGCGTATTGCCCGTGCGAATCACAAGCCAGTCTGGGTTGAAACCTGCCCAAAATATCTCCAGTTGGACGAACGCAGTTACGATACAGAACATTGAATGAACTTCAATCTTAGCCCAccc                                                 >  6:36334/2‑147 (MQ=255)
                                                                                              gTCTAACGTCTTAGGTCTGGATTATCTGCGTCTTGCACGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACa                                                                                  <  5:99246/112‑1 (MQ=255)
                                                                                               tCTAACGTCTTAGGTCTGGATTATATGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTAGAAACCTGCACACAATATCTCATGTTGGACGAACGCAGTTACGATACag                                                                                 >  6:99246/1‑112 (MQ=255)
                                                                                               tCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGc                                             >  5:46569/1‑148 (MQ=255)
                                                                                               tCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTAATTCTTAGCCCACCGCTGc                                             >  5:142045/1‑148 (MQ=255)
                                                                                               tCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGGTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTg                                              <  1:404020/147‑1 (MQ=255)
                                                                                                cTAACGGCTTAGCTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTACCCCACCGCTGCg                                            >  5:306811/1‑148 (MQ=255)
                                                                                                   aCGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTaa                                         <  5:390246/148‑1 (MQ=255)
                                                                                                      gCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACg                                       >  6:317251/1‑147 (MQ=255)
                                                                                                      gCTTAGGTATGGATTATCTGCGTCTCGCCCGTGCGAATCACCAGCCAGTCTTGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGt                                      <  5:59625/148‑1 (MQ=255)
                                                                                                                ggATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAg                            >  8:358202/1‑148 (MQ=255)
                                                                                                                 taTTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGATACCTGCCCACAATATCTCCTGTTGGTCGAACGCTGTTACGATACAGAATATGGCATGATGTTCATTCTTAGCGCACCGCTGCGTAACGTACGCGAGCa                             <  5:36334/145‑1 (MQ=255)
                                                                                                                 gATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCAt                                                              >  2:179509/1‑113 (MQ=255)
                                                                                                                 gATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCAt                                                              <  1:179509/113‑1 (MQ=255)
                                                                                                                 gATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCa                             <  5:152281/146‑1 (MQ=255)
                                                                                                                 gATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAgg                           <  5:317251/148‑1 (MQ=255)
                                                                                                                 gATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGGCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGc                              <  2:245908/145‑1 (MQ=255)
                                                                                                                     aTCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACa                        >  3:318761/1‑147 (MQ=255)
                                                                                                                      tCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACaa                       <  6:142045/147‑1 (MQ=255)
                                                                                                                       cTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGACACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAAc                     >  5:11674/1‑148 (MQ=255)
                                                                                                                        tGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAACCt                    >  4:199313/1‑148 (MQ=255)
                                                                                                                        tGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACt                    >  3:179719/1‑148 (MQ=255)
                                                                                                                        tGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACt                    >  2:30236/1‑148 (MQ=255)
                                                                                                                          cGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTAAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACtgt                  <  5:226170/148‑1 (MQ=255)
                                                                                                                            tCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACtgtgg                <  3:99415/148‑1 (MQ=255)
                                                                                                                                gCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACtgtggtg              <  5:103484/146‑1 (MQ=255)
                                                                                                                                gCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACtgtggt               <  8:158896/145‑1 (MQ=255)
                                                                                                                                  ccGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACTGTGGTGTGGc          <  2:211333/148‑1 (MQ=255)
                                                                                                                                     tGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACTGTGGTGTGGCAt        >  7:77520/1‑147 (MQ=255)
                                                                                                                                       cGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACTGTGGTGTGGCATCAg     <  2:356866/148‑1 (MQ=255)
                                                                                                                                       cGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACTGTGGTGTGGCATCAg     <  1:30236/148‑1 (MQ=255)
                                                                                                                                         aaTCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTTTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACTGTGGTGTGGCATCAGCg   <  6:109443/148‑1 (MQ=255)
                                                                                                                                          aTCACCAGCCAGTCTGGGTTGAAACCTGCCCACTATATCTCCTGTTGTTCTAACTCATTTACGATACAGATGATGGACTGAAGTTGATTCTTATGACACCGGTGCGTAACCTACCCGAGGAGTACAAACTGTGTTGTGTCTTCAGGGa  <  8:278978/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                             
TGAGTCGCCCTCTGGAATGCGAAGCGGAAGCCATCGCCCGCATGATTAACCTGGCACAAATTGCCGGTAACGCCCCGCTCTATATCGTGCACCTGTCTAACGGCTTAGGTCTGGATTATCTGCGTCTTGCCCGTGCGAATCACCAGCCAGTCTGGGTTGAAACCTGCCCACAATATCTCCTGTTGGACGAACGCAGTTACGATACAGAAGATGGCATGAAGTTCATTCTTAGCCCACCGCTGCGTAACGTACGCGAGCAGGACAAACTGTGGTGTGGCATCAGCGA  >  NC_000913/3010653‑3010938

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