Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,760,722 T→A 6.2% E176D (GAA→GAT)  selA ← selenocysteine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,760,7220TA 6.2% 224.4 / 3.2 75E176D (GAA→GAT) selAselenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base T (32/38);  new base A (3/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATG  >  NC_000913/3760576‑3760865
                                                                                                                                                  |                                                                                                                                               
aCTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTc                                                                                                                                               >  2:213660/1‑149 (MQ=255)
 cTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGTGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTc                                                                                                                                               <  1:112164/148‑1 (MQ=255)
 cTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCATACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTc                                                                                                                                               <  5:21043/148‑1 (MQ=255)
      atcaACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTcgccgc                                                                                                                                          >  8:178170/1‑148 (MQ=255)
       tcaACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCga                                                                                                                                        <  1:213660/149‑1 (MQ=255)
        caACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTcgccgc                                                                                                                                          >  8:86985/1‑146 (MQ=255)
         aaCAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGCAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGTGCCGCCAATCTCCACCAGTTcgccgc                                                                                                                                          >  7:104086/1‑145 (MQ=255)
                cGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCga                                                                                                                                        >  8:272531/1‑140 (MQ=255)
                cGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCga                                                                                                                                        <  7:272531/140‑1 (MQ=255)
                 ggTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATaccacca                                                                                                                              >  1:18241/1‑149 (MQ=255)
                      tttCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACctct                                                                                                                         >  3:115792/1‑149 (MQ=255)
                       ttcattcaCCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACctct                                                                                                                         <  5:80336/148‑1 (MQ=255)
                              aCCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCt                                                                                                                 >  1:106876/1‑149 (MQ=255)
                              aCCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCt                                                                                                                 <  7:263501/149‑1 (MQ=255)
                               ccGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTg                                                                                                                >  1:389935/1‑149 (MQ=255)
                                cGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTgg                                                                                                               <  8:87764/149‑1 (MQ=255)
                                      aCGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGGAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTg                                                                                                          <  2:389935/148‑1 (MQ=255)
                                        gATAATCATTCGCGTGCGTGCTGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAAAATCGGGAATAAGAAACGCGCCGCCAATCTCAACCAGATAGCAGCGATATACCACCACCTCTTATCCGCTGGCAgatct                                                                                                        >  2:330351/1‑144 (MQ=255)
                                        gATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCTCCGCCAATCTCCACCAGTTCGACGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGgcc                                                                                                       >  5:138860/1‑149 (MQ=255)
                                        gATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGgcc                                                                                                       <  7:86985/149‑1 (MQ=255)
                                           aaTCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGATCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGgccgc                                                                                                     >  2:30567/1‑148 (MQ=255)
                                           aaTCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGATCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGgccgc                                                                                                     >  2:37577/1‑148 (MQ=255)
                                             tCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCaa                                                                                                  <  8:214325/149‑1 (MQ=255)
                                             tCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGACTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCaa                                                                                                  <  1:69521/149‑1 (MQ=255)
                                               aTTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACa                                                                                                >  7:262730/1‑149 (MQ=255)
                                               aTTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACa                                                                                                >  3:59481/1‑149 (MQ=255)
                                                   gcgtgcgtGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCaa                                                                                            <  2:106876/149‑1 (MQ=255)
                                                    cgtgcgtgCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAAt                                                                                           <  8:262730/149‑1 (MQ=255)
                                                          tgcgGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACAt                                                                                               >  1:186573/1‑139 (MQ=255)
                                                          tgcgGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACAt                                                                                               <  2:186573/139‑1 (MQ=255)
                                                                     gtCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAtt                                                                          >  5:438172/1‑149 (MQ=255)
                                                                     gtCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAt                                                                           <  4:115792/148‑1 (MQ=255)
                                                                     gtCCCTACTTCGTGGAGGGTGCAGCCAGCCTGCCGCATAACAACGGGAATACGACACCCGCCGCCAATCTCCACCACTTCGCCGCGAGATACAACCCCCTCTTTTCCGCTCGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAtt                                                                          <  3:247819/149‑1 (MQ=255)
                                                                       cccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACCTCGGGAATCCGAAACGCGCCGCCAATCTCCCCCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATACCACCGCCGCCGCagggt                                                                        >  2:275067/1‑145 (MQ=255)
                                                                        ccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGt                                                                        >  6:332415/1‑148 (MQ=255)
                                                                        ccTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGt                                                                        >  6:161493/1‑148 (MQ=255)
                                                                         cTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTa                                                                      <  4:59481/149‑1 (MQ=255)
                                                                         cTACTTCGGGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTGCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTa                                                                      <  1:401406/149‑1 (MQ=255)
                                                                            cTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATa                                                                                          >  3:168677/1‑126 (MQ=255)
                                                                            cTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATa                                                                                          <  4:168677/126‑1 (MQ=255)
                                                                            cTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTAt                                                                     <  5:161493/147‑1 (MQ=255)
                                                                             ttCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGgccgc                                                                                                     >  5:169103/1‑114 (MQ=255)
                                                                             ttCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGgccgc                                                                                                     <  6:169103/114‑1 (MQ=255)
                                                                               cGTGTAGGGTTCAGCCTCCCTGACTCATAACATCGGGAATACGAAACCCGCCGCCAATCTCCACCAGTTCTCCGCGATATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGa                                                                  <  4:134335/147‑1 (MQ=255)
                                                                               cGTGTAGGGTGCAGCCTGCCTGACGCATATCATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGa                                                                  <  6:438172/147‑1 (MQ=255)
                                                                               cGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGa                                                                  <  3:178580/147‑1 (MQ=255)
                                                                               cGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGAc                                                                 <  7:178170/148‑1 (MQ=255)
                                                                               cGTGCAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGTGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGa                                                                  <  8:104086/147‑1 (MQ=255)
                                                                                gtgtAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGa                                                               <  7:420722/149‑1 (MQ=255)
                                                                                gtgtAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGa                                                               <  5:281833/149‑1 (MQ=255)
                                                                                      ggTGCAGCCTGCCTGACGCATAACCTCGGGAAAACGAAACGCGCCGCCAAGCTCCACCCGTTCGCCGCGAGATACCTCCACCTCTTATCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAg                                                          <  6:138860/148‑1 (MQ=255)
                                                                                      ggTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAgg                                                         <  5:332415/149‑1 (MQ=255)
                                                                                            gcctgcctGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAtt                                                                          >  1:451918/1‑126 (MQ=255)
                                                                                            gcctgcctGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAtt                                                                          <  2:451918/126‑1 (MQ=255)
                                                                                             cctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCaa                                                                                                  >  2:59792/1‑101 (MQ=255)
                                                                                             cctgcctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCaa                                                                                                  <  1:59792/101‑1 (MQ=255)
                                                                                                  ctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAg                                                          <  5:48489/136‑1 (MQ=255)
                                                                                                  ctgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAg                                                          >  6:48489/1‑136 (MQ=255)
                                                                                                   tgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATaa                                                                                         <  5:401351/104‑1 (MQ=255)
                                                                                                   tgACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATaa                                                                                         >  6:401351/1‑104 (MQ=255)
                                                                                                    gACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGccccc                                           >  3:421927/1‑149 (MQ=255)
                                                                                                       gCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTg                                                                   <  6:411603/122‑1 (MQ=255)
                                                                                                       gCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTg                                                                   >  5:411603/1‑122 (MQ=255)
                                                                                                       gCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGt                                         >  7:6234/1‑148 (MQ=255)
                                                                                                       gCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCCCCGCCGCATTGTAATTGACCATACAGGCATCTTCCGCCCCCGt                                         >  4:215258/1‑148 (MQ=255)
                                                                                                         aTAACATCGGGAATATGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTaa                                       >  8:246474/1‑148 (MQ=255)
                                                                                                           aaCATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGtt                                                                       >  6:170984/1‑114 (MQ=255)
                                                                                                           aaCATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGtt                                                                       <  5:170984/114‑1 (MQ=255)
                                                                                                               tCGGGAATACGAAACGCGCCGCCAATCTCCACCAGATCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGca                                <  1:37577/149‑1 (MQ=255)
                                                                                                               tCGGGAATACGAAACGCGCCGCCAATCTCCACCAGATCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGca                                <  1:30567/149‑1 (MQ=255)
                                                                                                                 gggAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGcaca                              <  2:18241/149‑1 (MQ=255)
                                                                                                                         gAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTgc                      >  8:365791/1‑149 (MQ=255)
                                                                                                                                   gccAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCAGTGTTAGTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGccc            >  2:459473/1‑149 (MQ=255)
                                                                                                                                            caccaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGAt   <  8:6234/149‑1 (MQ=255)
                                                                                                                                             accaGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATg  >  3:220621/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
ACTTTCATCAACAGTGCGGTATTTTCATTCACCGCCTGACGATAATCATTCGCGTGCGTGCGGTTGGTGGTCCCTACTTCGTGTAGGGTGCAGCCTGCCTGACGCATAACATCGGGAATACGAAACGCGCCGCCAATCTCCACCAGTTCGCCGCGAGATACCACCACCTCTTTTCCGCTGGCAGTGGCCGCCAACATCAATAACACCGCCGCCGCATTGTTATTGACGATACAGGCATCTTCCGCCCCCGTAATACGGCACAGCAGCTGCGCCAGCGCCCGATCGCGATG  >  NC_000913/3760576‑3760865

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: ‑