Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,534,388 A→T 5.9% Y524F (TAC→TTC)  pck → phosphoenolpyruvate carboxykinase [ATP]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,534,3880AT 5.9% 236.2 / 4.8 68Y524F (TAC→TTC) pckphosphoenolpyruvate carboxykinase [ATP]
Reads supporting (aligned to +/- strand):  ref base A (29/35);  new base T (2/2);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCCTTCTTTTGTC  >  NC_000913/3534240‑3534531
                                                                                                                                                    |                                                                                                                                               
atGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATa                                                                                                                                                 >  4:120636/2‑149 (MQ=255)
  gCCGATGTTTAACCTGGCGATCCCAACCGAGCTGCCGGCCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATaca                                                                                                                                               <  2:547197/149‑1 (MQ=255)
   ccGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATacac                                                                                                                                              <  2:218672/149‑1 (MQ=255)
    cGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATacac                                                                                                                                              <  2:414514/148‑1 (MQ=255)
    cGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATacac                                                                                                                                              >  1:414514/1‑148 (MQ=255)
      aTGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGa                                                                                                                                           <  3:119264/149‑1 (MQ=255)
         tttAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGacac                                                                                                                                        <  3:120636/149‑1 (MQ=255)
               cTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGc                                                                                                                                  >  4:273906/1‑149 (MQ=255)
               cTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGc                                                                                                                                  >  3:62126/1‑149 (MQ=255)
                tGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATTCACCGACACCCCTGCg                                                                                                                                 >  4:186350/1‑149 (MQ=255)
                tGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATTCACCGACACCCCTGCg                                                                                                                                 >  4:78579/1‑149 (MQ=255)
                      tCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTTcc                                                                                                                           >  2:19417/1‑149 (MQ=255)
                      tCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGcc                                                                                                                           >  1:455738/1‑149 (MQ=255)
                      tCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGcc                                                                                                                           <  1:216124/149‑1 (MQ=255)
                        ccAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCcgc                                                                                                                         >  4:88625/1‑149 (MQ=255)
                          aaCCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCgcgc                                                                                                                       >  1:318118/1‑149 (MQ=255)
                                  tGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCg                                                                                                               <  2:184992/149‑1 (MQ=255)
                                    ccGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGc                                                                                                             >  1:11724/1‑149 (MQ=255)
                                      gggCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTg                                                                                                           <  1:56630/149‑1 (MQ=255)
                                      gggCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTg                                                                                                           >  3:147933/1‑149 (MQ=255)
                                      gggCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTg                                                                                                           <  4:147933/149‑1 (MQ=255)
                                              acacGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAACTACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGaaa                                                                                                   >  4:45416/1‑149 (MQ=255)
                                                  gAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACAcc                                                                                                                                             >  3:105130/1‑103 (MQ=255)
                                                  gAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACAcc                                                                                                                                             <  4:105130/103‑1 (MQ=255)
                                                    aGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTaa                                                                                             <  2:263163/149‑1 (MQ=255)
                                                        tctcGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATacac                                                                                                                                              <  3:87500/96‑1 (MQ=255)
                                                        tctcGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATacac                                                                                                                                              >  4:87500/1‑96 (MQ=255)
                                                                 gcgTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTg                                                                                >  3:286705/1‑149 (MQ=255)
                                                                  cgTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTgg                                                                               <  2:373399/149‑1 (MQ=255)
                                                                   gTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGga                                                                              <  4:62126/149‑1 (MQ=255)
                                                                     aaCACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGgaga                                                                            >  3:9036/1‑149 (MQ=255)
                                                                        acCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAata                                                                         <  3:107678/149‑1 (MQ=255)
                                                                         ccTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAatat                                                                        >  2:532333/1‑149 (MQ=255)
                                                                         ccTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAatat                                                                        <  2:11724/149‑1 (MQ=255)
                                                                          cTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGAtt                                                                                        <  2:149952/132‑1 (MQ=255)
                                                                          cTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGAtt                                                                                        >  1:149952/1‑132 (MQ=255)
                                                                           tACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAAtatct                                                                      >  2:376383/1‑149 (MQ=255)
                                                                            aCGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAAtatcta                                                                     >  4:356039/1‑149 (MQ=255)
                                                                                ttCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCa                                                                 >  4:456639/1‑149 (MQ=255)
                                                                                 tctcCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAg                                                                >  4:114245/1‑149 (MQ=255)
                                                                                    ccGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATTCACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTAt                                                             <  3:78579/149‑1 (MQ=255)
                                                                                    ccGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATTCACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTAt                                                             <  3:186350/149‑1 (MQ=255)
                                                                                      ggtaCAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTACCt                                                           <  3:389101/146‑1 (MQ=255)
                                                                                       gAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAAtatctat                                                                    <  4:96903/139‑1 (MQ=255)
                                                                                       gAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAAtatctat                                                                    >  3:96903/1‑139 (MQ=255)
                                                                                        aaCAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTta                                                         <  1:152476/149‑1 (MQ=255)
                                                                                        aaCAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTta                                                         <  3:99297/149‑1 (MQ=255)
                                                                                        aaCAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTta                                                         <  3:550230/149‑1 (MQ=255)
                                                                                        aaCAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGAGTTGAAGCTGGAGAATATCTATCCAGTATCTta                                                         <  4:214054/149‑1 (MQ=255)
                                                                                          cAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTtata                                                       <  1:69823/149‑1 (MQ=255)
                                                                                              ggCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGaaa                                                   <  2:318118/149‑1 (MQ=255)
                                                                                               gCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAg                                                  >  3:154995/1‑149 (MQ=255)
                                                                                                 aGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAAtatctat                                                                    <  3:427347/129‑1 (MQ=255)
                                                                                                 aGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAAtatctat                                                                    >  4:427347/1‑129 (MQ=255)
                                                                                                   gAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCaaa                                              >  1:273042/1‑149 (MQ=255)
                                                                                                          ccGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAg                                       <  1:19417/149‑1 (MQ=255)
                                                                                                           cGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAgg                                      <  3:189899/149‑1 (MQ=255)
                                                                                                                 ccTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCtt                                <  3:273906/149‑1 (MQ=255)
                                                                                                                  cTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGcgaaggt                               <  2:97161/149‑7 (MQ=255)
                                                                                                                  cTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTc                               <  3:45416/149‑1 (MQ=255)
                                                                                                                        aaaCTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTcc                         >  3:405621/1‑149 (MQ=255)
                                                                                                                         aaCTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTccc                        <  4:172881/149‑1 (MQ=255)
                                                                                                                         aaCTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTccc                        <  3:88625/149‑1 (MQ=255)
                                                                                                                              tttATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGttta                   >  4:504892/1‑149 (MQ=255)
                                                                                                                                tATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGtttatt                 <  3:504892/149‑1 (MQ=255)
                                                                                                                                  tCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAg                                                                >  2:261482/1‑100 (MQ=255)
                                                                                                                                  tCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAg                                                                <  1:261482/100‑1 (MQ=255)
                                                                                                                                   cGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTAc              >  1:245467/1‑149 (MQ=255)
                                                                                                                                          ttCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGga                                                                              <  4:1217/78‑1 (MQ=255)
                                                                                                                                          ttCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGga                                                                              >  3:1217/1‑78 (MQ=255)
                                                                                                                                          ttCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCcttctt       >  1:60674/1‑149 (MQ=255)
                                                                                                                                          atCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGCAGAATATCTATCCACTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTCTATTTACCcttctt       >  4:507280/2‑149 (MQ=255)
                                                                                                                                            cGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCCTTCtttt     <  2:273042/149‑1 (MQ=255)
                                                                                                                                             gATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCCTTCTTTTg    <  3:456639/149‑1 (MQ=255)
                                                                                                                                               tAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCCTTCTTTTGTc  <  3:522238/149‑1 (MQ=255)
                                                                                                                                                    acacCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTAccc            >  1:325558/1‑134 (MQ=255)
                                                                                                                                                    acacCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTAccc            <  2:325558/134‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAACACCTACGCTTCTCCGGAACAGTGGCAGGAAAAAGCCGAAACCCTGGCGAAACTGTTTATCGACAACTTCGATAAATACACCGACACCCCTGCGGGTGCCGCGCTGGTAGCGGCTGGTCCGAAACTGTAATGATTTGAAGCTGGAGAATATCTATCCAGTATCTTATAGAAAGCAAAACGGGAGGCACCTTCGCCTCCCGTTTATTTACCCTTCTTTTGTC  >  NC_000913/3534240‑3534531

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑