Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,058 C→T 10.0% P241L (CCC→CTC)  thrA → Bifunctional aspartokinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,0580CT10.0% 170.8 / 8.4 60P241L (CCC→CTC) thrABifunctional aspartokinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base C (22/32);  new base T (3/3);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.34e-01

TGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTC  >  NC_000913/918‑1201
                                                                                                                                            |                                                                                                                                               
tGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGa                                                                                                                                         <  8:367055/149‑1 (MQ=255)
  aaaaaGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAgg                                                                                                                                       >  1:176615/1‑149 (MQ=255)
    aaaGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGgtt                                                                                                                                     >  5:398437/1‑149 (MQ=255)
    aaaGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGgtt                                                                                                                                     <  5:30387/149‑1 (MQ=255)
      aGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGgttgt                                                                                                                                   <  2:123512/149‑1 (MQ=255)
      aGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGgttgt                                                                                                                                   <  2:16384/149‑1 (MQ=255)
       ggCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGgttgtt                                                                                                                                  <  1:119229/149‑1 (MQ=255)
          gAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGaa                                                                                                                               <  1:372799/149‑1 (MQ=255)
                gtggtgCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCTCGATGCGAGGTTGTTGAAGTCGAt                                                                                                                         >  6:334403/1‑149 (MQ=255)
                gtggtgCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGAt                                                                                                                         <  2:205699/149‑1 (MQ=255)
                      cTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGGTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTa                                                                                                                   <  7:107532/149‑1 (MQ=255)
                         ggAAGCAACGGTTCCGACTACTCTGCTGCGGTGATGGCTGCCTGTTTACGCGACTATTGTTGCGAGTTTTGGACGGACGTTGACGGGGTCTATACATGCTAACCGCGTCAGGTGCTAGATGCGATGTTGTTTAAGTAGATTTCCTAACa                                                                                                                >  2:111247/1‑149 (MQ=255)
                          gACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAg                                                                                                               <  6:214114/149‑1 (MQ=255)
                               aaCGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGa                                                                                                                          <  4:275310/133‑1 (MQ=255)
                               aaCGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGa                                                                                                                          >  3:275310/1‑133 (MQ=255)
                                   gTTCCGACTAGTGTGCTGCGGTTGTGGTTGCCTGTTTTCGCGCAGATTGTTGGTAGATTTGGACGGTCGTTGACGGGGTCTATTCCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCTATGTCCTACCAGGAAGCGATg                                                                                                      <  6:398437/149‑1 (MQ=255)
                                      ccGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAg                                                                                                   >  1:156074/1‑149 (MQ=255)
                                         actactCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCtt                                                                                                <  1:338805/149‑1 (MQ=255)
                                          ctactCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGTGTCTATACCTGCGACCCGCGTCAGTTGTCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCttt                                                                                               <  8:361256/149‑1 (MQ=255)
                                             ctctGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTcc                                                                                             <  6:285662/148‑1 (MQ=255)
                                              tctGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCt                                                                                            >  6:179268/1‑148 (MQ=255)
                                              tctGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTa                                                                                           >  5:145816/1‑149 (MQ=255)
                                              gctgctgCGGTGCTGGCTGCCTGTTTACGCGCCGATTGGTGCGAGATCTCGCCGCACGTTGACGGGGGCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTa                                                                                           <  5:200117/148‑1 (MQ=255)
                                                  ctgcGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTc                                                                                       >  3:285543/1‑149 (MQ=255)
                                                  ctgcGGTGATGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGa                                                                                                    >  1:273306/1‑136 (MQ=255)
                                                  cggcGGTGCTGGCTGCCTGTTTACGCTCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGAGTCAGGTGCCCGATGCTAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGa                                                                                                    <  2:273306/134‑1 (MQ=255)
                                                     cGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCg                                                                                      <  6:145816/147‑1 (MQ=255)
                                                          ctggctgCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCTCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCgg                                                                                     <  4:35540/143‑1 (MQ=255)
                                                          ctggctgCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCTCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCgg                                                                                     >  3:35540/1‑143 (MQ=255)
                                                           tggctgCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTaaa                                                                              <  4:30192/149‑1 (MQ=255)
                                                           tggctgCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTa                                                                                >  4:79305/1‑147 (MQ=255)
                                                             gctgCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTaaa                                                                              <  8:23300/147‑1 (MQ=255)
                                                                 ccTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGttct                                                                         <  5:179268/148‑1 (MQ=255)
                                                                 ccTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATCGAGCTTTCCTACTTCGGCGCTAAAGttctt                                                                        >  5:353945/1‑149 (MQ=255)
                                                                  cTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGttctt                                                                        >  5:117064/1‑148 (MQ=255)
                                                                   tGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGttcttga                                                                      >  2:311716/1‑147 (MQ=255)
                                                                     tttACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCAc                                                                     <  8:43396/148‑1 (MQ=255)
                                                                     tttACGCGCCGATTGGTGCGAGATTTGGACGGACGTTGACGGTGTCTCTACCTGCGACCCGCGTCAGGTGCTCGATGCGAGGTTGTTCAATTCGATGTCCTAACAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCAc                                                                     <  1:111247/148‑1 (MQ=255)
                                                                          gcgcCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAACGTTCTTCACCCCCGc                                                               >  8:168819/1‑149 (MQ=255)
                                                                            gcCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTAc                                                                                                                  <  5:129453/96‑1 (MQ=255)
                                                                            gcCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTAc                                                                                                                  >  6:129453/1‑96 (MQ=255)
                                                                             ccGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGAGGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCAcc                                                            <  2:117611/149‑1 (MQ=255)
                                                                              cGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGGCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCa                                                           <  2:176615/149‑1 (MQ=255)
                                                                                    ttgCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTAcc                                                      <  6:117064/148‑1 (MQ=255)
                                                                                                      agtggACGGGGTCTATTCCTGCGACCCGCGTCAGGGGCCCGATGCGAGGTTGTTGTAGTCGATGTCCGACCACGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCAGCCCCCGAATCGCACAGTTCCAGa                                   <  6:356790/145‑1 (MQ=255)
                                                                                                                 tCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCCCCCCCGCACCATTACGCCCAGCGCCCAGGTCCAGATCCCGGGCCt                         >  4:300875/1‑148 (MQ=255)
                                                                                                                 tCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCGCCATTACCCCCCTCGCCCAGTTCGAGATCCCTTGCCt                         >  2:378439/1‑148 (MQ=255)
                                                                                                                     taCCTGCGACCCGTGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACGCCCGCCGCCTTCCCCGCCGCGCCCAGTGCCAGATGCCGTGCCGGAtga                    >  8:406961/1‑147 (MQ=255)
                                                                                                                     taCCTGCGACCCGCGTCAGGTGCTCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGAtt                     <  5:334403/148‑1 (MQ=255)
                                                                                                                     taCCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACGCCCGGCGCAGGGGCCCCCTCGCCCAGTGCCAGCTGGCGTGCCTgctg                     >  8:51700/1‑145 (MQ=255)
                                                                                                                     taCCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGa                       <  3:300875/146‑1 (MQ=255)
                                                                                                                        cTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCCCCGGCGCACCAGTGCCCCCCGCGGCCCCGGCCAGATCCCGTGCCTGaggaca                  >  5:157530/1‑143 (MQ=255)
                                                                                                                        cTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTaaaa                 <  1:187033/149‑1 (MQ=255)
                                                                                                                         tGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCAc                                                                     >  7:421117/1‑96 (MQ=255)
                                                                                                                         tGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCAc                                                                     <  8:421117/96‑1 (MQ=255)
                                                                                                                            gACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATAc             <  4:96234/149‑1 (MQ=255)
                                                                                                                             aCCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATAcc            <  3:79305/149‑1 (MQ=255)
                                                                                                                                 gcgTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCCCCCGCGCACCCTTCCCCCCAGGGCCCTGGTCCAGATCCCTGGACGAGTTAAAAAAACCGGa         >  3:267751/1‑148 (MQ=255)
                                                                                                                                 gcgTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTTCCTGATTAAAAATACCGGaa        <  1:382598/149‑1 (MQ=255)
                                                                                                                                        gTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCGCCAGGGCCCCCCTCGCCCAGTGCCAGATCCCAGGCCAGAGTAAAAAAACCCGAAATcaac  >  2:159350/1‑145 (MQ=255)
                                                                                                                                            |                                                                                                                                               
TGAAAAAGGCGAACTGGTGGTGCTTGGACGCAACGGTTCCGACTACTCTGCTGCGGTGCTGGCTGCCTGTTTACGCGCCGATTGTTGCGAGATTTGGACGGACGTTGACGGGGTCTATACCTGCGACCCGCGTCAGGTGCCCGATGCGAGGTTGTTGAAGTCGATGTCCTACCAGGAAGCGATGGAGCTTTCCTACTTCGGCGCTAAAGTTCTTCACCCCCGCACCATTACCCCCATCGCCCAGTTCCAGATCCCTTGCCTGATTAAAAATACCGGAAATCCTC  >  NC_000913/918‑1201

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