Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,827,198 G→C 5.5% A269A (GCG→GCC srlE → glucitol/sorbitol‑specific enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,827,1980GC 5.5% 514.0 / 15.0 144A269A (GCG→GCCsrlEglucitol/sorbitol‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base G (67/69);  new base C (4/4);  total (71/73)
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.71e-01

Only 100 of 149 total aligned reads displayed.

TCTGTTCCTTCCCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTAT  >  NC_000913/2827050‑2827345
                                                                                                                                                    |                                                                                                                                                   
tCTGTTCCTTCCCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAggcg                                                                                                                                                     >  2:870006/1‑149 (MQ=255)
    ttccttccCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCt                                                                                                                                                 <  3:913526/149‑1 (MQ=255)
    ttccttccCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCt                                                                                                                                                 >  2:1005016/1‑149 (MQ=255)
      ccttccCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGc                                                                                                                                               <  1:1005016/149‑1 (MQ=255)
            cACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTc                                                                                                                                         >  3:1154604/1‑149 (MQ=255)
            cACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTc                                                                                                                                         >  3:339192/1‑149 (MQ=255)
             aCTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCa                                                                                                                                        >  1:748733/1‑149 (MQ=255)
                 cTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTc                                                                                                                                         >  4:664165/1‑144 (MQ=255)
                 cTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTc                                                                                                                                         <  3:664165/144‑1 (MQ=255)
                  tttCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGAc                                                                                                                                            <  4:900206/140‑1 (MQ=255)
                  tttCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACtt                                                                                                                                          >  3:898425/1‑142 (MQ=255)
                  tttCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATc                                                                                                                                      >  3:46718/1‑146 (MQ=255)
                  tttCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATc                                                                                                                                      <  4:46718/146‑1 (MQ=255)
                     cACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtc                                                                                                                                <  1:945842/149‑1 (MQ=255)
                      aCCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGAATTGATTGGCGTGAAGATTGGTCTCGGCAATATTCCGCCGCAGCTGGCTTTACCGGCACTGTTTGCCATCAGCGCGCAGGCGGCCTGCGACTTCACCCcggtcg                                                                                                                               >  3:912574/1‑149 (MQ=255)
                        cTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcggt                                                                                                                             >  3:249062/1‑149 (MQ=255)
                        cTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcggt                                                                                                                             <  4:1154604/149‑1 (MQ=255)
                         tttCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGtt                                                                                                                            >  3:194496/1‑149 (MQ=255)
                         tttCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGtt                                                                                                                            >  3:1190090/1‑149 (MQ=255)
                         tttCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGtt                                                                                                                            <  4:249062/149‑1 (MQ=255)
                             cTCGGCCCAGGCGCAGTTATCGCAGAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTc                                                                                                                        <  2:954456/149‑1 (MQ=255)
                              tCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCg                                                                                                                       >  2:1411946/1‑149 (MQ=255)
                               cGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGc                                                                                                                      <  4:194496/149‑1 (MQ=255)
                                  cccATGCGCAGTTATCGCACAGGTTATCGGCGTATTGATAGGCGTGCCGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGTACTGTTTTCCATCAACGCGCAGGCGGCCTGCGACTTGATCCCGGTCGGTTTGTCGCTgg                                                                                                                   <  4:912574/149‑1 (MQ=255)
                                   ccAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGc                                                                                                                  >  3:240032/1‑149 (MQ=255)
                                     aGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGttt                                                                                                                           >  4:325095/1‑138 (MQ=255)
                                       gcgcAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcggt                                                                                                                             >  3:315504/1‑134 (MQ=255)
                                       gcgcAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcggt                                                                                                                             <  4:315504/134‑1 (MQ=255)
                                           aGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGccc                                                                                                          <  4:677260/149‑1 (MQ=255)
                                            gTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCg                                                                                                         >  3:167295/1‑149 (MQ=255)
                                             ttATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGt                                                                                                        >  3:1221133/1‑149 (MQ=255)
                                               aTCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCa                                                                                                      <  2:547579/149‑1 (MQ=255)
                                                tCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAg                                                                                                     <  2:883516/149‑1 (MQ=255)
                                                   cacaGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACCTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGac                                                                                                  >  3:230686/1‑149 (MQ=255)
                                                     caGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGacac                                                                                                <  2:299019/149‑1 (MQ=255)
                                                       ggTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTc                                                                                                                                         >  4:1352824/1‑106 (MQ=255)
                                                       ggTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACCTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACgg                                                                                              <  4:230686/149‑1 (MQ=255)
                                                        gTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGt                                                                                             >  2:563580/1‑149 (MQ=255)
                                                        gTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGt                                                                                             <  2:311949/149‑1 (MQ=255)
                                                          tATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTc                                                                                           >  1:633262/1‑149 (MQ=255)
                                                           aTCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCCGCCTGCgcgt                                                                                                                                           <  1:953282/100‑4 (MQ=255)
                                                           aTCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCCGCCTGCgcgt                                                                                                                                           >  2:953282/1‑97 (MQ=255)
                                                            tCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTcgc                                                                                         <  2:633262/149‑1 (MQ=255)
                                                               gCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGt                                                                                                        <  3:644728/131‑1 (MQ=255)
                                                               gCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGt                                                                                                        >  4:644728/1‑131 (MQ=255)
                                                               gCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTc                                                                                      <  3:1206223/149‑1 (MQ=255)
                                                                  tATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCg                                                                                     <  3:229049/147‑1 (MQ=255)
                                                                    ttgattgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGgtgt                                                                                 >  1:760551/1‑149 (MQ=255)
                                                                      gattgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTcc                                                                               <  4:167295/149‑1 (MQ=255)
                                                                        ttgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGaa                                                                                                              >  2:213068/1‑116 (MQ=255)
                                                                        ttgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGaa                                                                                                              <  1:213068/116‑1 (MQ=255)
                                                                        ttgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGccc                                                                                                          >  2:1219129/1‑120 (MQ=255)
                                                                        ttgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTcc                                                                               >  3:571277/1‑147 (MQ=255)
                                                                        ttgGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTcc                                                                               <  4:571277/147‑1 (MQ=255)
                                                                           gCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGccc                                                                                                          <  4:1353687/117‑1 (MQ=255)
                                                                           gCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGccc                                                                                                          >  3:1353687/1‑117 (MQ=255)
                                                                              tGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTa                                                                       <  3:269527/149‑1 (MQ=255)
                                                                                 aGATTGGTCACGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTg                                                                    >  1:1255198/1‑149 (MQ=255)
                                                                                  gATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATcc                                                                                                                                     >  2:182918/1‑83 (MQ=255)
                                                                                   aTTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcggt                                                                                                                             <  3:1452785/90‑1 (MQ=255)
                                                                                   aTTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcggt                                                                                                                             >  4:1452785/1‑90 (MQ=255)
                                                                                   aTTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcg                                                                                                                               >  3:453456/1‑88 (MQ=255)
                                                                                   aTTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCcggtcg                                                                                                                               <  4:453456/88‑1 (MQ=255)
                                                                                   aTTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGt                                                                  >  4:168986/1‑149 (MQ=255)
                                                                                   aTTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGt                                                                  <  1:563580/149‑1 (MQ=255)
                                                                                   aTTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGt                                                                  <  2:760551/149‑1 (MQ=255)
                                                                                     tGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGa                                                                <  1:1411946/149‑1 (MQ=255)
                                                                                        tctcGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGa                                                                >  1:1030271/1‑146 (MQ=255)
                                                                                        tctcGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGa                                                                <  2:1030271/146‑1 (MQ=255)
                                                                                        tctcGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGcc                                                             >  4:792504/1‑149 (MQ=255)
                                                                                             gCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCttt                                                        >  3:1036418/1‑149 (MQ=255)
                                                                                                    tCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCg                                                 <  4:353227/149‑1 (MQ=255)
                                                                                                       gccgcATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCt                                                                          >  4:1190863/1‑121 (MQ=255)
                                                                                                        ccgcATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGgcgc                                             <  4:1400421/149‑1 (MQ=255)
                                                                                                             tCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGa                                        <  3:168986/149‑1 (MQ=255)
                                                                                                                ggCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTg                                     >  4:1195109/1‑149 (MQ=255)
                                                                                                                   tttACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGttt                                                                                                                           >  4:1064869/1‑60 (MQ=255)
                                                                                                                   tttACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGttt                                                                                                                           <  3:1064869/60‑1 (MQ=255)
                                                                                                                   tttACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTAc                                  >  4:223326/1‑149 (MQ=255)
                                                                                                                        cGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATc                             <  2:748733/149‑1 (MQ=255)
                                                                                                                         ggCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCg                            <  3:1195109/149‑1 (MQ=255)
                                                                                                                          gCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGc                           >  1:46749/1‑149 (MQ=255)
                                                                                                                           cACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGcc                          >  2:1324157/1‑149 (MQ=255)
                                                                                                                           cACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGcc                          <  3:792504/149‑1 (MQ=255)
                                                                                                                              tGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCtttt                                                       >  1:1301493/1‑117 (MQ=255)
                                                                                                                              tGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCtttt                                                       <  2:1301493/117‑1 (MQ=255)
                                                                                                                                tttGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTAc                                                                      >  1:629940/1‑100 (MQ=255)
                                                                                                                                tttGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTAc                                                                      <  2:629940/100‑1 (MQ=255)
                                                                                                                                 ttGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGttt                    <  4:1036418/149‑1 (MQ=255)
                                                                                                                                     cATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCgg                                                                                    <  3:1244434/81‑1 (MQ=255)
                                                                                                                                     cATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGt                                                                  >  1:1113106/1‑99 (MQ=255)
                                                                                                                                      aTCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGtctc               <  4:120646/149‑1 (MQ=255)
                                                                                                                                       tCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTGGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTcc              <  1:1324157/149‑1 (MQ=255)
                                                                                                                                         aaCGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCg                                                 >  4:254759/1‑112 (MQ=255)
                                                                                                                                          agtcgcAGGCCGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGc                                                              <  3:1357467/95‑1 (MQ=255)
                                                                                                                                           gtcgcAGGCCGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGa                               <  4:1104724/126‑1 (MQ=255)
                                                                                                                                            tcgcAGGCCGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGttt         >  3:1431012/2‑149 (MQ=255)
                                                                                                                                             cgcAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGtttt        <  1:258585/149‑1 (MQ=255)
                                                                                                                                              gcAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTta       <  4:1221133/149‑1 (MQ=255)
                                                                                                                                                   cggcCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTtatctat  <  4:239213/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TCTGTTCCTTCCCACTGCTTTCACCTTTCCTCGGCCCAGGCGCAGTTATCGCACAGGTTATCGGCGTATTGATTGGCGTGCAGATTGGTCTCGGCAATATTCCGCCGCATCTGGCTTTACCGGCACTGTTTGCCATCAACGCGCAGGCGGCCTGCGACTTCATCCCGGTCGGTTTGTCGCTGGCGGAAGCCCGTCAGGACACGGTTCGCGTCGGTGTCCCTTCTGTACTGGTGAGCCGCTTTTTAACCGGCGCGCCGACTGTACTGATCGCCTGGTTTGTCTCCGGTTTTATCTAT  >  NC_000913/2827050‑2827345

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: