Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,967,686 G→T 9.1% S250* (TCG→TAG)  ptsP ← PEP‑protein phosphotransferase enzyme I; GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,967,6860GT 9.1% 302.1 / 15.9 88S250* (TCG→TAG) ptsPPEP‑protein phosphotransferase enzyme I; GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base G (41/39);  new base T (5/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

AATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGA  >  NC_000913/2967543‑2967834
                                                                                                                                               |                                                                                                                                                    
aaTCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGc                                                                                                                                                 >  1:574124/1‑149 (MQ=255)
  tCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAg                                                                                                                                               <  3:915360/149‑1 (MQ=255)
  tCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAg                                                                                                                                               >  2:673395/1‑149 (MQ=255)
  tCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCAGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAg                                                                                                                                               <  4:763387/149‑1 (MQ=255)
        gCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGa                                                                                                                                         <  3:1863/149‑1 (MQ=255)
           cGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTa                                                                                                                                      <  1:86966/149‑1 (MQ=255)
            gCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTaa                                                                                                                                     <  2:186572/149‑1 (MQ=255)
            gCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTaa                                                                                                                                     >  1:777592/1‑149 (MQ=255)
                tttGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGAt                                                                                                                                 >  2:691603/1‑149 (MQ=255)
                   gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTaaa                                                                                                                                    <  4:86766/143‑1 (MQ=255)
                   gagaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTaaa                                                                                                                                    >  3:86766/1‑143 (MQ=255)
                    agaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGaa                                                                                                                             >  3:611734/1‑149 (MQ=255)
                     gaTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGaaa                                                                                                                            >  2:799662/1‑149 (MQ=255)
                      aTAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGaaaa                                                                                                                           >  3:187048/1‑149 (MQ=255)
                       tAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAAt                                                                                                                          >  2:524678/1‑149 (MQ=255)
                          ttATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAgc                                                                                                                       <  4:244651/149‑1 (MQ=255)
                          ttATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAgc                                                                                                                       >  4:108255/1‑149 (MQ=255)
                            aTCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAgccg                                                                                                                     >  4:52870/1‑149 (MQ=255)
                                 ttAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGtt                                                                                                                >  4:641339/1‑149 (MQ=255)
                                  tAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGttt                                                                                                               <  1:772538/149‑1 (MQ=255)
                                  tAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGttt                                                                                                               <  4:498903/149‑1 (MQ=255)
                                    gcgcGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCt                                                                                                             <  2:5512/149‑1 (MQ=255)
                                        ggCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCttttt                                                                                                         <  3:269890/149‑1 (MQ=255)
                                        ggCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCttttt                                                                                                         <  3:52870/149‑1 (MQ=255)
                                        ggCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCttttt                                                                                                         <  3:904656/149‑1 (MQ=255)
                                          cAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTtgt                                                                                                       <  1:799662/149‑1 (MQ=255)
                                              cTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCTAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCAc                                                                                                   <  3:182392/149‑1 (MQ=255)
                                                 ttCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTTCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCttttt                                                                                                         >  4:536637/1‑140 (MQ=255)
                                                 ttCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCttttt                                                                                                         <  3:536637/140‑1 (MQ=255)
                                                     gCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCggcggc                                                                                            >  1:418189/1‑149 (MQ=255)
                                                     gCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCggcggc                                                                                            >  2:220866/1‑149 (MQ=255)
                                                      cAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCa                                                                                           >  1:584844/1‑149 (MQ=255)
                                                        aaTTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCaaa                                                                                         <  2:817553/149‑1 (MQ=255)
                                                        aaTTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCaaa                                                                                         <  4:729411/149‑1 (MQ=255)
                                                          tttttCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGc                                                                                       >  3:709640/1‑149 (MQ=255)
                                                            tttCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGt                                                                                     >  3:487346/1‑149 (MQ=255)
                                                             ttCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACGTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGtt                                                                                    <  3:762225/149‑1 (MQ=255)
                                                             ttCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGtt                                                                                    <  4:675141/149‑1 (MQ=255)
                                                              tCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGttt                                                                                   <  2:418189/149‑1 (MQ=255)
                                                               cAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTg                                                                                  >  3:175881/1‑149 (MQ=255)
                                                                aaTGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGc                                                                                 <  2:524967/149‑1 (MQ=255)
                                                                 aTGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCt                                                                                >  1:148942/1‑149 (MQ=255)
                                                                 aTGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCt                                                                                >  2:252239/1‑149 (MQ=255)
                                                                 aTGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGGTTGCt                                                                                >  2:107074/1‑149 (MQ=255)
                                                                  tGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTg                                                                               <  1:524678/149‑1 (MQ=255)
                                                                    aCCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTa                                                                             >  2:392667/1‑149 (MQ=255)
                                                                    aCCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTa                                                                             >  2:721669/1‑149 (MQ=255)
                                                                     ccGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAg                                                                            >  3:673885/1‑149 (MQ=255)
                                                                     ccGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGACGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTGTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAg                                                                            <  4:590984/149‑1 (MQ=255)
                                                                        tttttACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAgcgg                                                                         >  2:57568/1‑149 (MQ=255)
                                                                        tttttACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAgcgg                                                                         >  4:406830/1‑149 (MQ=255)
                                                                        tttttACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAgcgg                                                                         >  1:32514/1‑149 (MQ=255)
                                                                              ccGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAAc                                                                   <  4:142433/149‑1 (MQ=255)
                                                                               cGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGaaa                                                                    >  3:462382/1‑147 (MQ=255)
                                                                               cGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGaaa                                                                    <  4:462382/147‑1 (MQ=255)
                                                                                gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTc                                                                 <  3:513267/149‑1 (MQ=255)
                                                                                gCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTc                                                                 >  4:331713/1‑149 (MQ=255)
                                                                                 cccACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCg                                                                >  3:492006/1‑149 (MQ=255)
                                                                                 cccACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCg                                                                <  2:651247/149‑1 (MQ=255)
                                                                                    aCTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGttt                                                             >  3:332985/1‑149 (MQ=255)
                                                                                    aCTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGttt                                                             <  3:688905/149‑1 (MQ=255)
                                                                                          ccACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGct                                                       <  2:315167/149‑1 (MQ=255)
                                                                                            aCCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGcttc                                                     >  4:163580/1‑149 (MQ=255)
                                                                                             ccGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGaaa                                                                    <  2:929024/133‑1 (MQ=255)
                                                                                             ccGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGaaa                                                                    >  1:929024/1‑133 (MQ=255)
                                                                                              cGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCTAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGcttctt                                                   >  2:709544/1‑149 (MQ=255)
                                                                                                    tttATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTtgt                                                                                                       >  3:135899/1‑91 (MQ=255)
                                                                                                    tttATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTtgt                                                                                                       <  4:135899/91‑1 (MQ=255)
                                                                                                       aTCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTcc                                                 <  1:496870/142‑1 (MQ=255)
                                                                                                       aTCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTcc                                                 >  2:496870/1‑142 (MQ=255)
                                                                                                       aTCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGccc                                          <  4:332985/149‑1 (MQ=255)
                                                                                                          aaCCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCgg                                       >  4:926184/1‑149 (MQ=255)
                                                                                                          aaCCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCgg                                       >  3:120896/1‑149 (MQ=255)
                                                                                                            ccTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTc                                     >  2:679632/1‑149 (MQ=255)
                                                                                                                    aaCAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTc                             <  1:158368/149‑1 (MQ=255)
                                                                                                                    aaCAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTc                             >  1:453151/1‑149 (MQ=255)
                                                                                                                    aaCAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTc                             <  1:673395/149‑1 (MQ=255)
                                                                                                                      cAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGATGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTcgc                           <  4:611734/149‑1 (MQ=255)
                                                                                                                        aTTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGt                         >  2:762518/1‑149 (MQ=255)
                                                                                                                         ttCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGtt                        >  4:506453/1‑149 (MQ=255)
                                                                                                                            gcgACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCa                     >  4:31530/1‑149 (MQ=255)
                                                                                                                              gACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCaga                   <  3:926184/149‑1 (MQ=255)
                                                                                                                                cGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGc                 <  4:175881/149‑1 (MQ=255)
                                                                                                                                   aGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCgg              <  2:584844/149‑1 (MQ=255)
                                                                                                                                            tCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGt     >  3:185742/1‑149 (MQ=255)
                                                                                                                                              cGAAAGCAGGTGCGATTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTg   >  3:313728/1‑149 (MQ=255)
                                                                                                                                              cGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTg   <  2:32514/149‑1 (MQ=255)
                                                                                                                                              cGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTg   <  4:187048/149‑1 (MQ=255)
                                                                                                                                               gAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGa  <  4:487346/149‑1 (MQ=255)
                                                                                                                                               gAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGa  <  2:248903/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
AATCGCCAGCCCGCTCTTTGAGATAGTTATCACTTAGCGCGGCAAACTGTTCGGCAAATTTTTCAATGACCGTTTTTACCGCCCACTCTGCCACCGAGCCTTTATCAACCTCGGCAAACAATTCGCGACGCAGCCGGGTGTCCGAAAGCAGGTGCGAGTAAAGATCGAAAATAGCCGCCGTTTCTTTTTGTGCACCGGCGGCAAAGCGTTTGCTGTAGCGGCGAAACTCGTTTGCCGCTTCTTCCAGCGCCCCGGTCAGTCGTTCGCGTTCCAGAGCCGGATCCAGCGTTGA  >  NC_000913/2967543‑2967834

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