Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 602,046 C→G 5.0% L30V (CTG→GTG)  pheP → phenylalanine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913602,0460CG 5.0% 277.0 / 5.0 80L30V (CTG→GTG) phePphenylalanine transporter
Reads supporting (aligned to +/- strand):  ref base C (37/39);  new base G (2/2);  total (39/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.44e-01

GACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGT  >  NC_000913/601898‑602194
                                                                                                                                                    |                                                                                                                                                    
gacgaTGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAAc                                                                                                                                                      >  3:36299/1‑149 (MQ=255)
  cgaTGGTCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGTAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTg                                                                                                                                                    <  2:22287/149‑1 (MQ=255)
     tGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGAtt                                                                                                                                                 <  4:228879/149‑1 (MQ=255)
     tGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGAtt                                                                                                                                                 <  3:82669/149‑1 (MQ=255)
     tGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGAtt                                                                                                                                                 <  3:249258/149‑1 (MQ=255)
       ggaaCTAAATGTTAAAAGGTGCCCCTCAACAAATATCACACACAGGGGAAAGGCGTGAAAACCGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTAATATTCAACTGAtt                                                                                                                                                 <  8:75026/144‑1 (MQ=255)
        gCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCg                                                                                                                                              >  1:20071/1‑149 (MQ=255)
        gCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCg                                                                                                                                              >  6:91668/1‑149 (MQ=255)
          aCTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGTAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGtt                                                                                                                                            <  4:344742/149‑1 (MQ=255)
           cTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTAACGGAAGATACTGAGTCGCATCAAGAGCCGACGCTTCATCCCGGATTACATAACCGTCATATTCAACTGATTTCGTTg                                                                                                                                           >  5:152951/1‑149 (MQ=255)
            tAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTgg                                                                                                                                          >  5:144027/1‑149 (MQ=255)
             aaaTGTTAAAAGGGGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTggg                                                                                                                                         <  1:120957/149‑1 (MQ=255)
                 gTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTggg                                                                                                                                         >  4:178412/1‑145 (MQ=255)
                 gTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTggg                                                                                                                                         <  3:178412/145‑1 (MQ=255)
                   tAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGgcg                                                                                                                                    >  4:335043/1‑148 (MQ=255)
                    aaaaGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGgcgc                                                                                                                                   >  5:333131/1‑148 (MQ=255)
                     aaaGGGCCCCCTCAACAAAACGGACACACAGGGGAAAGGCGTCAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCaa                                                                                                                                 <  8:18728/149‑1 (MQ=255)
                      aaGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCTTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAAt                                                                                                                                <  5:156486/149‑1 (MQ=255)
                      aaGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAAt                                                                                                                                >  5:349959/1‑149 (MQ=255)
                        gttGCCCCTCAACAAAAAACACACACAGGGGAAATGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTg                                                                                                                              <  3:334300/147‑1 (MQ=255)
                        ggTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGCTACTGCGTCGAATCAAGAGCCGACGCTTCGTCGCGGATTACATAACCGTCATATTCAACTGATTGCGGTGCGTGGCGCAAGTg                                                                                                                              >  3:74170/1‑149 (MQ=255)
                        ggTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTg                                                                                                                              <  6:333131/149‑1 (MQ=255)
                        ggTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTg                                                                                                                              <  6:349959/149‑1 (MQ=255)
                           gCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTa                                                                                                                           <  6:144027/149‑1 (MQ=255)
                           gCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTa                                                                                                                           >  7:238618/1‑149 (MQ=255)
                            ccccTCAAAAAAAAAGTCCCAATGGGTCATTGCGTTAAAAACTCTTCAACCGTCGCGTAAGATACTGGGTCTAATCACGATCCGACGCTTCATCGCGGATTTCATAACCTTCATATTATACTGATTGCTTTGGGTGGCTCAATTGGTa                                                                                                                           <  4:278533/148‑1 (MQ=255)
                              ccTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTg                                                                                                                                       >  7:321359/1‑134 (MQ=255)
                               cTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAAGTGATTGCGTTGGGTGGCGCAATTGGTACTgg                                                                                                                       >  7:215389/1‑149 (MQ=255)
                               cTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAAGTGATTGCGTTGGGTGGCGCAATTGGTACTgg                                                                                                                       >  7:220909/1‑149 (MQ=255)
                                tCAACAAAAAAGCCACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATAGCGTTGGGTGGCGCAATTGGTACTGGt                                                                                                                      >  3:333939/1‑149 (MQ=255)
                                tCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGATTCATCGCGGATTACATAAGCGTAATATTCAACTGAGTGCGTTGGGTGTCGCAATTGGTACTGGt                                                                                                                      >  5:300693/1‑149 (MQ=255)
                                   acaaAAAAGACACACAGGGGAAAGGCTTGAAAAACGCGTCACCCCTAACGGAAGATACTGCGTCGAATCAACCCCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTg                                                                                                                   <  6:230985/149‑1 (MQ=255)
                                       aaaaGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTc                                                                                                               <  8:238618/149‑1 (MQ=255)
                                             cacacaGGTTAAAGGCGTGAAAAACGCGTCAACCTTATCGGAAGATACTGCGTCTAATAATTATCCGACGCTTCATCGCTGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCa                                                                                                         <  6:152951/149‑1 (MQ=255)
                                              acacaGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTa                                                                                                                           >  4:356170/1‑130 (MQ=255)
                                              acacaGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTa                                                                                                                           <  3:356170/130‑1 (MQ=255)
                                                acaGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGTATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTg                                                                                                      >  5:197944/1‑149 (MQ=255)
                                                 caGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTgg                                                                                                     >  2:366415/1‑149 (MQ=255)
                                                    gtgAAAGGCGGGAAAAACGCGGAAACCTTAGCGGAAGATACTGCGTCGAAGCACGAGCCGACGCTTCATCGCGGATTACATAACCGGCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGccc                                                                                                  <  4:248542/147‑1 (MQ=255)
                                                    gggAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGCGCCGACGCTTCATCGCGGATTACATAACCGGCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGccc                                                                                                  >  2:9126/1‑149 (MQ=255)
                                                    gggAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGccc                                                                                                  >  4:124604/1‑149 (MQ=255)
                                                       aaaGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGc                                                                                               <  5:91668/149‑1 (MQ=255)
                                                            cGTGAAAAACGCGTCAACCGTATCGAAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTc                                                                                          <  1:289068/149‑1 (MQ=255)
                                                               gAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAAtt                                                                                                                               >  5:168156/1‑109 (MQ=255)
                                                                 aaaaCGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATg                                                                                     <  2:20071/149‑1 (MQ=255)
                                                                   aaCGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGc                                                                                   <  3:166844/149‑1 (MQ=255)
                                                                     cgcgTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCgg                                                                                 <  3:297846/149‑1 (MQ=255)
                                                                          cAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCg                                                                            >  1:138880/1‑149 (MQ=255)
                                                                             ccGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCt                                                                         <  3:335043/149‑1 (MQ=255)
                                                                              cgttacGGAAGATACTGCGTCGACTCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTGGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTg                                                                        <  6:197944/144‑1 (MQ=255)
                                                                               gTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAAGTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGt                                                                       <  8:220909/149‑1 (MQ=255)
                                                                               gTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAAGTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGt                                                                       <  8:215389/149‑1 (MQ=255)
                                                                                tATCGGAATATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTa                                                                      <  1:366415/149‑1 (MQ=255)
                                                                                  tCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGc                                                                                                          >  1:248729/1‑111 (MQ=255)
                                                                                  tCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGc                                                                                                          <  2:248729/111‑1 (MQ=255)
                                                                                          tACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACAGAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGc                                                             >  6:241686/1‑148 (MQ=255)
                                                                                               cGTCGACTCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGc                                                       <  1:339984/149‑1 (MQ=255)
                                                                                                gTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCg                                                      >  8:337961/1‑149 (MQ=255)
                                                                                                 tCGAATCAAGAGCCGACTCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTCGTCTGTTTCTTGGCATTGGCCCGCCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGt                                                     >  4:108845/1‑149 (MQ=255)
                                                                                                 tCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGt                                                     >  5:14291/1‑149 (MQ=255)
                                                                                                  cGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTc                                                    >  5:228313/1‑149 (MQ=255)
                                                                                                      tCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCg                                                <  3:56857/149‑1 (MQ=255)
                                                                                                         agagCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGa                                             >  7:387243/1‑149 (MQ=255)
                                                                                                             ccGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCgg                                                                           >  4:157888/1‑115 (MQ=255)
                                                                                                             ccGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCgg                                                                           <  3:157888/115‑1 (MQ=255)
                                                                                                              cGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGatcatc                                        <  6:14291/149‑1 (MQ=255)
                                                                                                               gACGCTTCTTCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGGCTGTTGCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCg                                       <  4:74170/149‑1 (MQ=255)
                                                                                                                aCGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGc                                      >  1:60227/1‑149 (MQ=255)
                                                                                                                      cATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTc                                  <  2:60227/147‑1 (MQ=255)
                                                                                                                        tCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCACAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGa                              >  1:40117/1‑149 (MQ=255)
                                                                                                                              aTTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATgc                        <  3:124604/149‑1 (MQ=255)
                                                                                                                                  cATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATgcgc                      >  7:139221/1‑147 (MQ=255)
                                                                                                                                   aTAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAg                   >  5:189311/1‑149 (MQ=255)
                                                                                                                                    tAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGc                  >  1:25872/1‑149 (MQ=255)
                                                                                                                                     aaCCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCt                 >  4:211812/1‑149 (MQ=255)
                                                                                                                                        cGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTgg              >  6:349141/1‑149 (MQ=255)
                                                                                                                                         gTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGc             <  7:337961/149‑1 (MQ=255)
                                                                                                                                           cATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCg            >  2:357536/1‑148 (MQ=255)
                                                                                                                                                 cAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAAtggt     >  5:31065/1‑149 (MQ=255)
                                                                                                                                                  aaCTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAAtggtg    <  6:31065/149‑1 (MQ=255)
                                                                                                                                                   aCTGATTGCGTTGGGGGGCGCAATTGGTACTGGTCTGTTTCTTGGCAGTGGCCCGGCGATTCAGATGGCGGGGCCGGCGGTATTGCTGGGCTACGGCGGCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAAtggtgg   <  8:285081/149‑1 (MQ=255)
                                                                                                                                                    cTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATggtggt  <  5:349141/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GACGATGGGCACTAAATGTTAAAAGGTGCCCCTCAACAAAAAAGACACACAGGGGAAAGGCGTGAAAAACGCGTCAACCGTATCGGAAGATACTGCGTCGAATCAAGAGCCGACGCTTCATCGCGGATTACATAACCGTCATATTCAACTGATTGCGTTGGGTGGCGCAATTGGTACTGGTCTGTTTCTTGGCATTGGCCCGGCGATTCAGATGGCGGGTCCGGCTGTATTGCTGGGCTACGGCGTCGCCGGGATCATCGCTTTCCTGATTATGCGCCAGCTTGGCGAAATGGTGGT  >  NC_000913/601898‑602194

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: