Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,147,304 G→C 5.0% R56P (CGG→CCG)  frwD → putative enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,147,3040GC 5.0% 558.8 / 14.1 156R56P (CGG→CCG) frwDputative enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (71/77);  new base C (3/5);  total (74/82)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.61e-01

Only 100 of 162 total aligned reads displayed.

AACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCC  >  NC_000913/4147158‑4147451
                                                                                                                                                  |                                                                                                                                                   
aaCCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCggg                                                                                                                                                   >  1:806113/1‑149 (MQ=255)
    gCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGa                                                                                                                                               <  4:1266107/149‑1 (MQ=255)
    gCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGa                                                                                                                                               <  3:109722/149‑1 (MQ=255)
    gCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGa                                                                                                                                               <  1:1240345/149‑1 (MQ=255)
    gCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCACTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGa                                                                                                                                               <  1:234565/149‑1 (MQ=255)
     ccTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGc                                                                                                                                                  >  1:787853/1‑145 (MQ=255)
      cTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAta                                                                                                                                             <  1:6946/149‑2 (MQ=255)
           tAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatcctc                                                                                                                                        <  1:104406/149‑7 (MQ=255)
            aaGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGctc                                                                                                                                       >  2:1058913/1‑149 (MQ=255)
                ggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGtt                                                                                                                                   <  1:1058913/149‑1 (MQ=255)
                ggcgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGtt                                                                                                                                   <  4:254775/149‑1 (MQ=255)
                  cgTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGa                                                                                                                                 <  1:224810/149‑1 (MQ=255)
                   gTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGAt                                                                                                                                <  4:1076970/149‑1 (MQ=255)
                   gTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGAt                                                                                                                                <  4:1423329/149‑1 (MQ=255)
                    tGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGAtt                                                                                                                               <  1:334513/149‑1 (MQ=255)
                     ggCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTa                                                                                                                              <  4:461052/149‑1 (MQ=255)
                      gcgcATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTAc                                                                                                                             <  4:1253373/149‑1 (MQ=255)
                          aTACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGat                                                                                                                         >  2:877888/1‑149 (MQ=255)
                          aTACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGat                                                                                                                         <  4:1012906/149‑1 (MQ=255)
                             cTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGa                                                                                                                          <  1:428238/145‑1 (MQ=255)
                             cTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGa                                                                                                                          >  2:428238/1‑145 (MQ=255)
                              ttATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCg                                                                                                                     >  1:594319/1‑149 (MQ=255)
                               tataTGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGa                                                                                                                    <  4:702905/149‑1 (MQ=255)
                                         cggAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGt                                                                                                          <  4:123882/149‑1 (MQ=255)
                                           gAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGc                                                                                                        <  2:1232886/149‑1 (MQ=255)
                                             aCGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCg                                                                                                      >  2:772462/1‑149 (MQ=255)
                                                gCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGc                                                                                                   >  1:969556/1‑149 (MQ=255)
                                                 cTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCg                                                                                                  >  2:698428/1‑149 (MQ=255)
                                                   ggAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAt                                                                                                >  3:898776/1‑149 (MQ=255)
                                                   ggAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAt                                                                                                <  4:215094/149‑1 (MQ=255)
                                                    gAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAg                                                                                                                   >  1:1073637/1‑129 (MQ=255)
                                                    gAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGAtt                                                                                               <  2:612603/149‑1 (MQ=255)
                                                       aaGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGa                                                                                            <  2:969556/149‑1 (MQ=255)
                                                        aGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGaa                                                                                           <  4:499844/149‑1 (MQ=255)
                                                         gTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGc                                                                                                              <  3:1605146/129‑1 (MQ=255)
                                                         gTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGc                                                                                                              >  4:1605146/1‑129 (MQ=255)
                                                              gCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATc                                                                                                                      <  2:100150/116‑1 (MQ=255)
                                                              gCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATc                                                                                                                      >  1:100150/1‑116 (MQ=255)
                                                                   ttAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTa                                                                                >  1:819543/1‑149 (MQ=255)
                                                                     agagAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGtt                                                                                                                                   <  1:569474/96‑1 (MQ=255)
                                                                     agagAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGtt                                                                                                                                   >  2:569474/1‑96 (MQ=255)
                                                                     agagAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatc                                                                                                                                           <  1:1092753/88‑4 (MQ=255)
                                                                     agagAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCCGGCGGAtatc                                                                                                                                           >  2:1092753/1‑85 (MQ=255)
                                                                      gagaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAtgt                                                                             >  1:259925/1‑149 (MQ=255)
                                                                      gagaAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCGGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTACCGGTGCCGAGCGATTTGAACATTGCCGCTAtgt                                                                             >  1:1405334/1‑149 (MQ=255)
                                                                          aGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCaa                                                                         <  3:798614/149‑1 (MQ=255)
                                                                           gTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAAt                                                                        <  1:877888/149‑1 (MQ=255)
                                                                            tGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAt                                                                               >  1:550781/1‑141 (MQ=255)
                                                                            tGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTAt                                                                               <  2:550781/141‑1 (MQ=255)
                                                                              gggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                     <  3:644737/149‑1 (MQ=255)
                                                                              gggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                     <  4:1621848/149‑1 (MQ=255)
                                                                              gggAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATgca                                                                     <  1:376596/149‑1 (MQ=255)
                                                                                  gTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCAt                                                                 >  1:415913/1‑149 (MQ=255)
                                                                                   tCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATc                                                                >  2:18654/1‑149 (MQ=255)
                                                                                     aGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTa                                                              >  2:405284/1‑149 (MQ=255)
                                                                                     aGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTa                                                              <  3:968270/149‑1 (MQ=255)
                                                                                      gCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTAc                                                             <  3:1013284/149‑1 (MQ=255)
                                                                                         ttGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCa                                                          >  2:1682451/1‑149 (MQ=255)
                                                                                         ttGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCa                                                          >  2:1501079/1‑149 (MQ=255)
                                                                                          tGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCAt                                                         <  1:1501079/149‑1 (MQ=255)
                                                                                              aCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCt                                                     <  3:653693/149‑1 (MQ=255)
                                                                                              aCTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCt                                                     <  3:1072820/149‑1 (MQ=255)
                                                                                                tCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGc                                                   >  4:103048/1‑149 (MQ=255)
                                                                                                tCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGc                                                   <  1:772462/149‑1 (MQ=255)
                                                                                                tCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGc                                                   <  3:1392587/149‑1 (MQ=255)
                                                                                                   gggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACAt                                                                                        >  4:134754/1‑109 (MQ=255)
                                                                                                   gggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTg                                                >  2:1288512/1‑149 (MQ=255)
                                                                                                   gggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTg                                                >  2:700938/1‑149 (MQ=255)
                                                                                                   gggCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTg                                                >  1:194768/1‑149 (MQ=255)
                                                                                                         gcTGGGAACGGAGAATCTTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCGCTGTTGATTACGGGTATCGAGCTTGCCGGTGCCGAGCGATTTTAACCTTGCCGCTACGTGCCATGCAGCATCTACGCATTCCTGCGTGAGCCGc                                          >  4:937541/1‑149 (MQ=255)
                                                                                                            gggAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGc                                       >  4:294056/1‑149 (MQ=255)
                                                                                                             ggAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCg                                      >  2:311217/1‑149 (MQ=255)
                                                                                                              gAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCaa                                                                         >  3:461053/1‑113 (MQ=255)
                                                                                                              gAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCgg                                     >  3:1571272/1‑149 (MQ=255)
                                                                                                                   gagaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAAt                                <  1:1670998/149‑1 (MQ=255)
                                                                                                                   gagaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAAt                                <  3:1442748/149‑1 (MQ=255)
                                                                                                                     gaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGa                              <  2:377069/149‑1 (MQ=255)
                                                                                                                     gaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGa                              <  3:543222/149‑1 (MQ=255)
                                                                                                                     gaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGa                              <  2:596725/149‑1 (MQ=255)
                                                                                                                     gaATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGa                              <  2:930822/149‑1 (MQ=255)
                                                                                                                      taTCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGTGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAg                             >  4:1092183/2‑149 (MQ=255)
                                                                                                                          gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGcgcg                         >  1:1303472/1‑149 (MQ=255)
                                                                                                                          gTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGcgcg                         >  2:1451186/1‑149 (MQ=255)
                                                                                                                            ttAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGt                       <  4:622635/149‑1 (MQ=255)
                                                                                                                            ttAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGt                       <  3:258632/149‑1 (MQ=255)
                                                                                                                            ttAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGt                       <  3:430561/149‑1 (MQ=255)
                                                                                                                            ttAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGt                       <  1:362542/149‑1 (MQ=255)
                                                                                                                               gCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgc                     <  1:1424706/148‑1 (MQ=255)
                                                                                                                               gCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTgc                     >  2:1424706/1‑148 (MQ=255)
                                                                                                                                     agcaacatCCGCCCGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCaaa                >  3:1310378/7‑147 (MQ=255)
                                                                                                                                     agcaacatCCGCCCGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCaaa                <  4:1310378/141‑1 (MQ=255)
                                                                                                                                       ggATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGc            >  3:685327/1‑149 (MQ=255)
                                                                                                                                       ggATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGc            >  4:1000627/1‑149 (MQ=255)
                                                                                                                                         atatCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCtt          <  2:483831/149‑1 (MQ=255)
                                                                                                                                          tatCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCttt         <  2:1303472/149‑1 (MQ=255)
                                                                                                                                           atCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTc        >  1:1656226/1‑149 (MQ=255)
                                                                                                                                            tCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCt       <  4:447861/149‑1 (MQ=255)
                                                                                                                                               gccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTgcg    >  1:294326/1‑149 (MQ=255)
                                                                                                                                                ccgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTgcgc   <  3:721257/149‑1 (MQ=255)
                                                                                                                                                 cgGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGcgcc  <  3:108291/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AACCGCCTGCGTAAGCGGCGTGGCGCATACTTATATGGCGGCGGAACGGCTGGAAAAGTTGTGCCTGTTAGAGAAGTGGGGAGTCAGCATTGAAACTCAGGGCGCGCTGGGAACGGAGAATCGTTTAGCGGATGAGGATATCCGCCGGGCGGATGTTGCTCTGTTGATTACGGATATCGAGCTTGCCGGTGCCGAGCGATTTGAACATTGCCGCTATGTGCAATGCAGCATCTACGCATTCCTGCGTGAGCCGCAGCGGGTAATGAGCGCGGTGCGCAAAGTGCTTTCTGCGCC  >  NC_000913/4147158‑4147451

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: