Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,421,967 A→G 5.1% K87R (AAG→AGG)  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,9670AG 5.1% 248.5 / 3.1 80K87R (AAG→AGG) ulaCL‑ascorbate‑specific enzyme IIA component of PTS
Reads supporting (aligned to +/- strand):  ref base A (40/34);  new base G (2/2);  total (42/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACG  >  NC_000913/4421820‑4422111
                                                                                                                                                   |                                                                                                                                                
tgGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAg                                                                                                                                                 <  8:271467/149‑1 (MQ=255)
  gTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGcc                                                                                                                                               >  5:123915/1‑149 (MQ=255)
     gagCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGc                                                                                                                                             <  5:103023/148‑1 (MQ=255)
        ccGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGa                                                                                                                                         <  4:259050/149‑1 (MQ=255)
             ttACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAGGCCGCTGGAGTTc                                                                                                                                     >  8:358640/1‑148 (MQ=255)
               actacCAGGCGATTCTGGCTGGCGTTGAACAGTTCGGTCCTTATTTCGTTCTCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCaa                                                                                                                                   <  8:166529/148‑1 (MQ=255)
                  acCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCAc                                                                                                                               >  6:242423/1‑149 (MQ=255)
                    cAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                                                                                                             >  6:211962/1‑149 (MQ=255)
                    cAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCCTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGa                                                                                                                             <  1:173642/149‑1 (MQ=255)
                     aGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCAGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTATGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGAt                                                                                                                            >  2:286154/1‑149 (MQ=255)
                         gATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAGGCCGCTGGAGTTCAACCACGATGAc                                                                                                                         >  8:266607/1‑148 (MQ=255)
                            tCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGCCAAc                                                                                                                      >  6:206634/1‑148 (MQ=255)
                               ggATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATAGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGCTGACAACGATc                                                                                                                  >  2:391730/1‑149 (MQ=255)
                                    gCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAGGCCGCTGGAGTTCAACCACGATGACAACGATCCGGt                                                                                                              <  7:358640/148‑1 (MQ=255)
                                     cGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTgg                                                                                                            <  6:123915/149‑1 (MQ=255)
                                          caCAGTTCCGTCCGTATTTCGTTATCGAGCCGGGCCTGGCAATGCCGCACGGGCGGCCGGCAGCGGGCGTTAAGAACACCGGTTTCGCTCTGGTAACGCTGAAAACGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATc                                                                                                       <  2:244872/148‑1 (MQ=255)
                                          aaCAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTg                                                                                                             >  1:54775/1‑143 (MQ=255)
                                          aaCAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTg                                                                                                             <  2:54775/143‑1 (MQ=255)
                                          aaCAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATc                                                                                                       <  6:290270/149‑1 (MQ=255)
                                           aCAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATcc                                                                                                      >  6:390737/1‑149 (MQ=255)
                                            cAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCt                                                                                                     >  1:317083/1‑149 (MQ=255)
                                             aGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCtc                                                                                                    <  5:295017/149‑1 (MQ=255)
                                             aGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCtc                                                                                                    <  1:321874/149‑1 (MQ=255)
                                             aGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCt                                                                                                     >  6:295017/1‑148 (MQ=255)
                                               ttCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGa                                                                                                                    <  2:147571/131‑1 (MQ=255)
                                               ttCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGa                                                                                                                    >  1:147571/1‑131 (MQ=255)
                                                  ggTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGat                                                                                                          >  7:157839/1‑138 (MQ=255)
                                                  ggTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGat                                                                                                          <  8:157839/138‑1 (MQ=255)
                                                    tCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAGGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCa                                                                                             <  7:266607/149‑1 (MQ=255)
                                                      cTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATg                                                                                           >  3:70764/1‑149 (MQ=255)
                                                      cTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATg                                                                                           >  1:272881/1‑149 (MQ=255)
                                                       ttATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGt                                                                                                              >  7:110963/1‑129 (MQ=255)
                                                       ttATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGt                                                                                                              <  8:110963/129‑1 (MQ=255)
                                                          tttCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGCGTTCACCCACGATGACAGCGATCCGGTGGATATCCTCATCACCATggcgg                                                                                       >  8:385735/1‑149 (MQ=255)
                                                          tttCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATggcgg                                                                                       >  6:401025/1‑149 (MQ=255)
                                                            tCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATggcggc                                                                                      <  5:211962/148‑1 (MQ=255)
                                                              gTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGt                                                                                   >  5:37132/1‑149 (MQ=255)
                                                               ttATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTc                                                                                  <  2:272881/149‑1 (MQ=255)
                                                                        cGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCaa                                                                          <  5:390737/148‑1 (MQ=255)
                                                                        cGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGCGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGAGGATATCCTCATGACCATGGCGGCGGTCGATGCCaa                                                                          >  7:165101/1‑148 (MQ=255)
                                                                         gggCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATa                                                                        <  2:317083/149‑1 (MQ=255)
                                                                              tGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGCTCGATGCCAATACTCAc                                                                   >  4:41149/1‑149 (MQ=255)
                                                                               ggCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCAc                                                                   >  7:21633/1‑148 (MQ=255)
                                                                                  aaTGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGTGGTCGATGCCAATACTCACCAg                                                                >  6:158026/1‑148 (MQ=255)
                                                                                      ccGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGt                                                           >  5:217189/1‑149 (MQ=255)
                                                                                        gCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGTGGTCGATGCCAATACTCACCAGGAAGTgg                                                         <  5:158026/149‑1 (MQ=255)
                                                                                          aCGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTggg                                                        >  4:71022/1‑148 (MQ=255)
                                                                                           cGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGca                                                      <  3:369311/149‑1 (MQ=255)
                                                                                                gTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATg                                                 >  3:90916/1‑149 (MQ=255)
                                                                                                 tCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATg                                                 >  5:188813/1‑148 (MQ=255)
                                                                                                 tCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGc                                                <  5:242423/149‑1 (MQ=255)
                                                                                                   cGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCa                                               <  6:217189/148‑1 (MQ=255)
                                                                                                    ggAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGAt                                                                               >  8:134693/1‑115 (MQ=255)
                                                                                                    ggAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGAt                                                                               <  7:134693/115‑1 (MQ=255)
                                                                                                    ggAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATg                                                 >  4:614/1‑145 (MQ=255)
                                                                                                          gggCGTTAAGAAAACCGGTTTCTGTCTGGTAACGCTGAATATGCCGCGGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCTTCACCATCGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCa                                       <  8:165101/149‑1 (MQ=255)
                                                                                                          gggCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCa                                       <  5:401025/149‑1 (MQ=255)
                                                                                                            gCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATGTCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAAc                                     <  3:244902/149‑1 (MQ=255)
                                                                                                                 aaGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGt                                 <  3:41149/148‑1 (MQ=255)
                                                                                                                   gAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGcat                                                     <  8:376098/126‑1 (MQ=255)
                                                                                                                   gAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGcat                                                     >  7:376098/1‑126 (MQ=255)
                                                                                                                   gAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGttt                               >  7:375707/1‑148 (MQ=255)
                                                                                                                        ccGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaagat                         <  6:188813/149‑1 (MQ=255)
                                                                                                                         cGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaagatg                        <  6:37132/149‑1 (MQ=255)
                                                                                                                           gTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaagatgaa                      <  3:71022/149‑1 (MQ=255)
                                                                                                                             ttCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCGGTTtgaagatgaa                      >  7:270948/1‑147 (MQ=255)
                                                                                                                                  tctGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATTGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGt                                 <  2:406628/131‑1 (MQ=255)
                                                                                                                                  tctGGTAACGCTGAAAAAGCCGCTGGAGTTAAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGt                                 >  1:406628/1‑131 (MQ=255)
                                                                                                                                    tGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATg                                                 >  1:48468/1‑113 (MQ=255)
                                                                                                                                    tGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATg                                                 <  2:48468/113‑1 (MQ=255)
                                                                                                                                     ggTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAAtttt             >  6:236966/1‑148 (MQ=255)
                                                                                                                                     ggTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAAtttt             <  5:236966/148‑1 (MQ=255)
                                                                                                                                      gTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGt                                 >  7:246386/1‑127 (MQ=255)
                                                                                                                                      gTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGt                                 <  8:246386/127‑1 (MQ=255)
                                                                                                                                        aaCGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGAcc         >  3:89053/1‑149 (MQ=255)
                                                                                                                                           gCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGtc      >  4:41218/1‑148 (MQ=255)
                                                                                                                                             tGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaagatgaag                     >  8:15626/1‑132 (MQ=255)
                                                                                                                                             tGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTtgaagatgaag                     <  7:15626/132‑1 (MQ=255)
                                                                                                                                             tGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTa    >  8:246837/1‑149 (MQ=255)
                                                                                                                                             tGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTAACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTa    >  7:119232/1‑149 (MQ=255)
                                                                                                                                               aaaaaGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGAc          >  2:365925/1‑141 (MQ=255)
                                                                                                                                               aaaaaGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGAc          <  1:365925/141‑1 (MQ=255)
                                                                                                                                                aaaaGCAGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTAcg  <  4:89053/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TGGTAGAGCCGCGTTACTACCAGGCGATTCTGGATGGCGTTGAACAGTTCGGTCCTTATTTCGTTATCGCTCCGGGCCTGGCAATGCCGCACGGGCGTCCGGAAGAGGGCGTTAAGAAAACCGGTTTCTCTCTGGTAACGCTGAAAAAGCCGCTGGAGTTCAACCACGATGACAACGATCCGGTGGATATCCTCATCACCATGGCGGCGGTCGATGCCAATACTCACCAGGAAGTGGGCATCATGCAGATCGTCAACCTGTTTGAAGATGAAGAGAATTTTGACCGTTTACG  >  NC_000913/4421820‑4422111

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