Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,421,878 T→C 5.5% Y57Y (TAT→TAC ulaC → L‑ascorbate‑specific enzyme IIA component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,421,8780TC 5.5% 232.7 / 3.4 72Y57Y (TAT→TACulaCL‑ascorbate‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base T (41/27);  new base C (2/2);  total (43/29)
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.90e-01

CCTGCAGGCTGAAGCAGAGACATGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCG  >  NC_000913/4421749‑4422025
                                                                                                                                 |                                                                                                                                                   
ccTGCAGGCTGAAGCAGAGACATGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGc                                                                                                                                   >  5:53278/1‑148 (MQ=255)
ccTGCAGGCTGAAGCAGAGACATGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGc                                                                                                                                   >  7:264966/1‑148 (MQ=255)
   gCACGCTGAAGCAGAGACATGGCAGGAAGCGCTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGATATCGCTCCGGGCCTg                                                                                                                                >  8:191774/1‑148 (MQ=255)
     aGGCTGAAGCAGAGACATGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGTCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTgg                                                                                                                               >  3:84919/1‑147 (MQ=255)
          gAAGCAGAGACATGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGc                                                                                                                         >  5:208072/1‑148 (MQ=255)
                gagaCATGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACg                                                                                                                   >  4:232412/1‑148 (MQ=255)
                    cATGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCg                                                                                                               >  5:123249/1‑148 (MQ=255)
                      tGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCt                                                                                                                                         >  7:312105/1‑120 (MQ=255)
                      tGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCt                                                                                                                                         <  8:312105/120‑1 (MQ=255)
                      tGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTc                                                                                                             <  2:298100/148‑1 (MQ=255)
                            gAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAag                                                                                                       <  8:264966/148‑1 (MQ=255)
                              aGCGGTGACACTCGGCGTTGACCGGCTGGTGGCGGCAGAGGTGGTAGAGCCGCGTTCCTACCAGGCGATTCTGGATGGCGTTGCACAGTGCGGGCCTTATTTCGTTATCGCGCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAagag                                                                                                     <  5:56778/148‑1 (MQ=255)
                              aGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCg                                                                                                                       <  2:247139/130‑1 (MQ=255)
                              aGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCg                                                                                                                       >  1:247139/1‑130 (MQ=255)
                                cGGGGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAggg                                                                                                   <  3:20017/148‑1 (MQ=255)
                                  gTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCg                                                                                                 <  8:42746/148‑1 (MQ=255)
                                  gTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATGCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCg                                                                                                 <  8:290889/148‑1 (MQ=255)
                                     aaaaTCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGtt                                                                                               >  5:312503/1‑147 (MQ=255)
                                               ttGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGTACAGTTCGGTCCTTATTTCGATATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCgg                                                                                    <  7:191774/148‑1 (MQ=255)
                                                tGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTAt                                                                                                                                                     >  5:58899/1‑82 (MQ=255)
                                                tGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTAt                                                                                                                                                     <  6:58899/82‑1 (MQ=255)
                                                 gACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGtt                                                                                  <  6:53278/148‑1 (MQ=255)
                                                   ccTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTtc                                                                                >  5:272352/1‑148 (MQ=255)
                                                      gctgGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAag                                                                                                       >  3:268719/1‑122 (MQ=255)
                                                      gctgGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAag                                                                                                       <  4:268719/122‑1 (MQ=255)
                                                      gctgGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTtctct                                                                             >  3:30087/1‑148 (MQ=255)
                                                       ctgGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTtctctc                                                                            <  7:90990/148‑1 (MQ=255)
                                                         ggtggCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCCCGGGCGTCCGGAAGAGGGCGTTAAGCAACCCGGGTTCTCTCTg                                                                          >  8:125720/1‑148 (MQ=255)
                                                         ggtggCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTg                                                                          <  6:208072/148‑1 (MQ=255)
                                                           tggCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGcc                                                                                                                                  <  5:77422/90‑1 (MQ=255)
                                                           tggCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGcc                                                                                                                                  >  6:77422/1‑90 (MQ=255)
                                                           tggCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTgg                                                                         >  1:210486/1‑147 (MQ=255)
                                                              cggcAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAAc                                                                     >  8:187560/1‑148 (MQ=255)
                                                                gcAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAag                                                                                                       >  7:334528/1‑112 (MQ=255)
                                                                gcAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAag                                                                                                       <  8:334528/112‑1 (MQ=255)
                                                                gcAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGc                                                                   >  8:252915/1‑148 (MQ=255)
                                                                gcAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGc                                                                   >  7:300372/1‑148 (MQ=255)
                                                                          tAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGcc                                                         >  4:293465/1‑148 (MQ=255)
                                                                             agCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCt                                                      <  5:116157/148‑1 (MQ=255)
                                                                              gCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGa                                                                <  5:199400/137‑1 (MQ=255)
                                                                              gCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGATGa                                                                >  6:199400/1‑137 (MQ=255)
                                                                                         acCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAACCCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCt                                                      <  5:251111/136‑1 (MQ=255)
                                                                                         acCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCt                                                      >  6:251111/1‑136 (MQ=255)
                                                                                           cAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCAc                                         >  7:84714/1‑147 (MQ=255)
                                                                                            aGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCTCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                       >  2:44337/1‑148 (MQ=255)
                                                                                            aGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                       <  6:123249/148‑1 (MQ=255)
                                                                                            aGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                       >  2:141785/1‑148 (MQ=255)
                                                                                              gCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCa                                              >  7:65953/1‑139 (MQ=255)
                                                                                              gCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCa                                              <  8:65953/139‑1 (MQ=255)
                                                                                                      ggATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGAt                             >  7:327356/1‑148 (MQ=255)
                                                                                                       gATGGCGTTGATCAGTTCGGTCCTTATTTCGTTATGGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCTGAAGAGTGCGTTAAGTAAACCTGTTTCTCTCTGGTAACGCTGAAATAGACGCTGGAGTTCAACCACGATGACAACGAt                             <  4:84919/147‑1 (MQ=255)
                                                                                                       gATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGTCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATc                            >  1:174786/1‑148 (MQ=255)
                                                                                                            cGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTg                       <  2:210486/148‑1 (MQ=255)
                                                                                                               tGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGat                    >  2:211932/1‑148 (MQ=255)
                                                                                                               tGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGa                     >  6:62875/1‑147 (MQ=255)
                                                                                                               tGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGa                     >  2:207822/1‑147 (MQ=255)
                                                                                                               tGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGa                     >  4:117951/1‑147 (MQ=255)
                                                                                                                  aCAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATc                 <  8:300372/148‑1 (MQ=255)
                                                                                                                  aCAGTTCGGTCCTTACTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                       >  4:121921/1‑126 (MQ=255)
                                                                                                                  aCAGTTCGGTCCTTACTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                       <  3:121921/126‑1 (MQ=255)
                                                                                                                  aCAGTTCGGTCCTTACTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                       <  3:27249/126‑1 (MQ=255)
                                                                                                                  aCAGTTCGGTCCTTACTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                       >  4:27249/1‑126 (MQ=255)
                                                                                                                    aGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCt               >  2:20456/1‑148 (MQ=255)
                                                                                                                       tCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTTAACCACGATGACAACGATCCGGTGGATATCCtcat            >  4:241533/1‑148 (MQ=255)
                                                                                                                       tCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCtcat            >  2:324640/1‑148 (MQ=255)
                                                                                                                         ggTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCtcatca          <  3:117951/148‑1 (MQ=255)
                                                                                                                         ggTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCtcatca          <  6:312503/148‑1 (MQ=255)
                                                                                                                          gTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCAc         >  3:173504/1‑148 (MQ=255)
                                                                                                                          gTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCAc         >  1:161611/1‑148 (MQ=255)
                                                                                                                             cTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAag                                                                                                       <  7:97553/51‑1 (MQ=255)
                                                                                                                             cTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAag                                                                                                       >  8:97553/1‑51 (MQ=255)
                                                                                                                               tATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCg                          <  5:201193/126‑1 (MQ=255)
                                                                                                                               tATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCg                          >  6:201193/1‑126 (MQ=255)
                                                                                                                               tATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATgg    >  5:107025/1‑148 (MQ=255)
                                                                                                                                 tttCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATggcg  >  8:43567/1‑148 (MQ=255)
                                                                                                                                 tttCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATggc   <  6:272352/147‑1 (MQ=255)
                                                                                                                                 |                                                                                                                                                   
CCTGCAGGCTGAAGCAGAGACATGGCAGGAAGCGGTGAAAATCGGCGTTGACCTGCTGGTGGCGGCAGATGTGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCG  >  NC_000913/4421749‑4422025

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