Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,533,559 G→A 6.3% G248S (GGC→AGC)  pck → phosphoenolpyruvate carboxykinase [ATP]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,533,5590GA 6.3% 200.5 / 3.0 63G248S (GGC→AGC) pckphosphoenolpyruvate carboxykinase [ATP]
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

TGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCG  >  NC_000913/3533413‑3533705
                                                                                                                                                  |                                                                                                                                                  
tGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGc                                                                                                                                                  <  1:185530/149‑1 (MQ=255)
  cAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCt                                                                                                                                                >  7:44090/1‑149 (MQ=255)
            gcggcACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCAc                                                                                                                                       <  5:170765/148‑1 (MQ=255)
                     ggTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAac                                                                                                                              >  7:363586/1‑148 (MQ=255)
                         cggcggCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGATGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAaccacc                                                                                                                          <  6:370437/148‑1 (MQ=255)
                             ggcgAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTc                                                                                                                     <  3:354090/149‑1 (MQ=255)
                               cgAAATGAAGAAAGGGATGTTCTCGATTATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGAGGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCa                                                                                                                   <  5:355817/149‑1 (MQ=255)
                               cgAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCa                                                                                                                   <  5:232132/149‑1 (MQ=255)
                                        gaaAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCg                                                                                                           >  4:45769/1‑148 (MQ=255)
                                             ggATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGt                                                                                                     >  3:60231/1‑149 (MQ=255)
                                               aTGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCg                                                                                                   >  6:14895/1‑149 (MQ=255)
                                                 gTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAcccacc                                                                                                                          >  6:209909/1‑124 (MQ=255)
                                                 gTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGcc                                                                                                 >  6:347326/1‑149 (MQ=255)
                                                     tcGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCg                                                                                                                                                    >  8:156944/1‑94 (MQ=255)
                                                     tcGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCg                                                                                                                                                    <  7:156944/94‑1 (MQ=255)
                                                       gatgatGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGAtt                                                                                            <  5:14895/148‑1 (MQ=255)
                                                        atgatgAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTgg                                                                                          <  5:347326/149‑1 (MQ=255)
                                                        atgatgAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGAGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTgg                                                                                          <  8:44090/149‑1 (MQ=255)
                                                           atgaACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGa                                                                                       <  3:45769/149‑1 (MQ=255)
                                                            tgaACTACCTGCTGCCGCTGAAAGGTATCGCCTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGAt                                                                                      >  7:75038/1‑149 (MQ=255)
                                                             gaACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAATGAGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCATATCCACCGACCCGAAACGTCGCCTGATTGGCGAt                                                                                      >  3:52261/1‑148 (MQ=255)
                                                             gaACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGa                                                                                       <  2:223781/147‑1 (MQ=255)
                                                             gaACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGa                                                                                       >  1:223781/1‑147 (MQ=255)
                                                              caCTACCTGCTGCCGCTGAAAGGTATCGCTTCTAAGCACTGCTCCGCCAACGTTTGTGAGAAAGGCGATGTTGCGGTGGTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATg                                                                                     <  1:244812/147‑1 (MQ=255)
                                                              aaCTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGa                                                                                    <  6:125480/149‑1 (MQ=255)
                                                               aCTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGTTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGAc                                                                                   <  8:123487/149‑1 (MQ=255)
                                                               aCTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGAc                                                                                   <  4:373440/149‑1 (MQ=255)
                                                                 tACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGAtt                                                                                            >  5:31499/1‑138 (MQ=255)
                                                                 tACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGAtt                                                                                            <  6:31499/138‑1 (MQ=255)
                                                                   ccTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGaa                                                                                <  7:63725/148‑1 (MQ=255)
                                                                       ctgcCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACgg                                                                           >  6:129936/1‑149 (MQ=255)
                                                                        tgcCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACgg                                                                           >  3:230431/1‑148 (MQ=255)
                                                                          ccGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACgg                                                                           <  4:230431/146‑1 (MQ=255)
                                                                             cTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGa                                                                     <  8:122404/149‑1 (MQ=255)
                                                                               gAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGAc                                                                    >  1:111353/1‑148 (MQ=255)
                                                                                 aaGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGAt                                                                 >  2:272645/1‑149 (MQ=255)
                                                                                  aGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATg                                                                <  5:129936/149‑1 (MQ=255)
                                                                                     tATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACg                                                             <  6:10174/149‑1 (MQ=255)
                                                                                          cTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCgtg                                                        >  8:24918/1‑149 (MQ=255)
                                                                                             cTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGtt                                                      >  6:110104/1‑148 (MQ=255)
                                                                                              tATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCAGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGttt                                                     >  7:135816/1‑148 (MQ=255)
                                                                                              tATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCAGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGttt                                                     >  7:141867/1‑148 (MQ=255)
                                                                                               aTGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTGTTCGGCCTTTCCGGCCCCGGTATAACCACCATTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGCCGATGACGGCGGGttga                                                    >  7:371539/1‑146 (MQ=255)
                                                                                                tGCACTGCTCGGCCAACGTTGGTGAGAAAGGCGATGTTGCTGTGTTCTTCGGCCTTTCCGGCACCTGTAAAACCACCCTTTCCAACGACCCGAAACGTCTCCTGATTGGCGATGACGAACACGGCTGTGACGATGACGGCGTGTTTaa                                                   <  4:21905/148‑1 (MQ=255)
                                                                                                tGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTaa                                                   <  8:363586/148‑1 (MQ=255)
                                                                                                tGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAAc                                                  <  2:111353/149‑1 (MQ=255)
                                                                                                tGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAAc                                                  >  1:360268/1‑149 (MQ=255)
                                                                                                          cGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAggcg                                        >  4:168483/1‑149 (MQ=255)
                                                                                                           gCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGAGGACGGCGTGGTTAACGTCGAAggcgg                                       >  7:205516/1‑149 (MQ=255)
                                                                                                             cAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCAGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGgct                                     <  8:135816/149‑1 (MQ=255)
                                                                                                             cAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCAGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGgct                                     <  8:141867/149‑1 (MQ=255)
                                                                                                                     gTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGATCACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGc                              <  1:272645/148‑1 (MQ=255)
                                                                                                                         gaAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAggcgg                                       <  1:125170/135‑1 (MQ=255)
                                                                                                                         gaAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAggcgg                                       >  2:125170/1‑135 (MQ=255)
                                                                                                                          aaaGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCaaaaa                         <  7:361699/148‑1 (MQ=255)
                                                                                                                          aaaGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCaaaaa                         >  8:361699/1‑148 (MQ=255)
                                                                                                                          aaaGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAAc                        >  8:343742/1‑149 (MQ=255)
                                                                                                                                gATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCaa                  <  7:24918/149‑1 (MQ=255)
                                                                                                                                       cGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAACCTATCAAGCTGTc            >  3:239382/1‑148 (MQ=255)
                                                                                                                                        ggTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCaa                  <  2:318404/141‑1 (MQ=255)
                                                                                                                                        ggTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCaa                  >  1:318404/1‑141 (MQ=255)
                                                                                                                                        ggTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCg           >  4:353649/1‑148 (MQ=255)
                                                                                                                                          tgtTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGaaa        <  4:60231/149‑1 (MQ=255)
                                                                                                                                              cttcGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAaagaa     <  8:75038/148‑1 (MQ=255)
                                                                                                                                                 cGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCg  >  2:372868/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TGCAGCTGATTGGCGGCACCTGGTACGGCGGCGAAATGAAGAAAGGGATGTTCTCGATGATGAACTACCTGCTGCCGCTGAAAGGTATCGCTTCTATGCACTGCTCCGCCAACGTTGGTGAGAAAGGCGATGTTGCGGTGTTCTTCGGCCTTTCCGGCACCGGTAAAACCACCCTTTCCACCGACCCGAAACGTCGCCTGATTGGCGATGACGAACACGGCTGGGACGATGACGGCGTGTTTAACTTCGAAGGCGGCTGCTACGCAAAAACTATCAAGCTGTCGAAAGAAGCG  >  NC_000913/3533413‑3533705

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