Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,157,636 G→C 15.1% intergenic (+62/‑233) ycfH → / → ptsG putative DNase/fused glucose‑specific PTS enzymes: IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,157,6360GC15.1% 332.1 / 50.0 125intergenic (+62/‑233)ycfH/ptsGputative DNase/fused glucose‑specific PTS enzymes: IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (53/53);  new base C (9/10);  total (62/63)
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.20e-02

Only 100 of 136 total aligned reads displayed.

TTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACA  >  NC_000913/1157537‑1157735
                                                                                                   |                                                                                                   
ttCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTg                                                                                                     >  1:1624039/1‑100 (MQ=255)
ttCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTg                                                                                                     >  1:171636/1‑100 (MQ=255)
ttCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTg                                                                                                     <  1:2276585/100‑1 (MQ=255)
  cacaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTggg                                                                                                   <  1:744221/100‑1 (MQ=255)
  cacaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTggg                                                                                                   <  2:272868/100‑1 (MQ=255)
   acaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGc                                                                                                  <  1:3614296/100‑1 (MQ=255)
   acaTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGc                                                                                                  <  1:2048720/100‑1 (MQ=255)
    caTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCg                                                                                                 >  2:1069706/1‑100 (MQ=255)
    caTCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCg                                                                                                 >  1:2532307/1‑100 (MQ=255)
      tCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGt                                                                                               <  2:155257/100‑1 (MQ=255)
       cGACGCTTCCCGCCTTGAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTg                                                                                              <  2:2127344/100‑1 (MQ=255)
       cGACGCTTCCCGCCTTCAATCCATCCGTTGACTGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTg                                                                                              <  1:3173538/100‑1 (MQ=255)
       cGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTg                                                                                              <  2:2765986/100‑1 (MQ=255)
       cGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTg                                                                                              <  2:2850969/100‑1 (MQ=255)
        gACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAACGTTCGTACTTATTGGCTAAAACGAGTAAAGTTTACCGCCGCTAATTGGGCGGTGa                                                                                             >  2:93027/1‑100 (MQ=255)
        gACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGa                                                                                             >  1:3331699/1‑100 (MQ=255)
        gACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGa                                                                                             <  1:3244564/100‑1 (MQ=255)
        gACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGa                                                                                             >  2:3092023/1‑100 (MQ=255)
          cGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAAt                                                                                           >  2:2443418/1‑100 (MQ=255)
          cGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaac                                                                                           <  1:2131034/100‑2 (MQ=38)
          cGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaac                                                                                           <  1:3092046/100‑2 (MQ=38)
           gCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATa                                                                                          <  1:3549062/100‑1 (MQ=255)
           gCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaact                                                                                          <  2:2477276/100‑3 (MQ=38)
            cTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATaa                                                                                         >  1:184489/1‑100 (MQ=255)
            cTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGGAAAGTTCACAGCCGAAAATTGGGCGGTGAATaa                                                                                         >  2:2572681/1‑100 (MQ=255)
               cccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCa                                                                                      >  1:869277/1‑100 (MQ=255)
               cccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCa                                                                                      >  2:3383634/1‑100 (MQ=255)
               cccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaactttac                                                                                      <  2:573628/100‑7 (MQ=37)
                ccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCAc                                                                                     >  2:1805022/1‑100 (MQ=255)
                ccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaactttact                                                                                     <  1:1286558/100‑8 (MQ=25)
                ccGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCCAATTTTCGGCGGTGaactttact                                                                                     <  2:2698998/100‑8 (MQ=25)
                 cGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACg                                                                                    <  2:1881795/100‑1 (MQ=255)
                 cGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACg                                                                                    >  1:2969366/1‑100 (MQ=255)
                   ccTTCAATCCATCCGTTGAATGAGTTTTTTTAACGCTCGTAATTAATGGCTAAAACGCGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCAc                                                                                     >  2:112406/1‑97 (MQ=255)
                   ccTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACTGCCGAAAATTGGGCGGTGAATAACCACGtt                                                                                  >  2:2590029/1‑100 (MQ=255)
                      tCAATCCATCAGTTGAATGAGTACTTTTAAACCTCCTAATAAATGGCTAAAACGAGCAAAGTTCAACGCCGAAAATCGGACGAAGAATAACCACCTTCGa                                                                               >  2:2755990/1‑100 (MQ=255)
                       cAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGaa                                                                              >  1:803498/1‑100 (MQ=255)
                        aaTCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGaaa                                                                             >  2:2318778/1‑100 (MQ=255)
                        aaTCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAATTCACCGCCGAAAATTGGGCGGTGAATAACAACGTTTGaaa                                                                             >  2:3157731/1‑100 (MQ=255)
                         atccatccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAat                                                                            <  1:934244/100‑1 (MQ=255)
                            catccGTTGAATGAGTTTTTTTAAAGATCGTAATTAATGGCTAAAAAGAGTAAAGTTCACCGCCGAAAATTGGACGGTGAATAACCACGTTTGAAATAtt                                                                         >  2:2470832/1‑100 (MQ=255)
                             atccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtg                                                                        >  1:2259453/1‑100 (MQ=255)
                              tccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgt                                                                       >  1:1776529/1‑100 (MQ=255)
                              tccGTTGAATGAGGTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtg                                                                        <  1:3070065/99‑1 (MQ=255)
                               ccGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATtgtg                                                                      >  1:2534843/1‑100 (MQ=255)
                                cGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGGTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGa                                                                     <  2:2990580/100‑1 (MQ=255)
                                cGTCGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGa                                                                     >  2:2318120/1‑100 (MQ=255)
                                   tGAATGAGTTTTTTTTTTTTTTTTATTTAATGGCTAAAACGAGTAAAGTTTATCGTTGAATATTGGGCGGTGAATAACCACTTTTGAATTATTGTGACat                                                                  >  2:1155226/1‑100 (MQ=255)
                                   tGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACat                                                                  <  1:1805030/100‑1 (MQ=255)
                                         aGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGttt                                                            <  2:1716608/100‑1 (MQ=255)
                                          gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGtttt                                                           >  2:3066374/1‑100 (MQ=255)
                                          gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGtttt                                                           >  1:1138926/1‑100 (MQ=255)
                                          gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGtttt                                                           <  2:1263180/100‑1 (MQ=255)
                                          gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGtttt                                                           >  1:3011598/1‑100 (MQ=255)
                                          gTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGtttt                                                           >  2:10676/1‑100 (MQ=255)
                                           tttttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTg                                                          <  2:2009918/100‑1 (MQ=255)
                                           tttttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTg                                                          <  1:3383659/100‑1 (MQ=255)
                                            ttttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGtttt                                                           >  2:2184247/1‑98 (MQ=255)
                                              ttttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCa                                                       >  1:3255272/1‑100 (MQ=255)
                                               tttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCaa                                                      >  1:2906420/1‑100 (MQ=255)
                                               tttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCaa                                                      <  2:199707/100‑1 (MQ=255)
                                                ttAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCaaa                                                     >  1:565458/1‑100 (MQ=255)
                                                 tAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCaaaa                                                    >  1:1933297/1‑100 (MQ=255)
                                                  aaaGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAt                                                   <  2:898420/100‑1 (MQ=255)
                                                    aGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtgt                                                 >  1:1366184/1‑100 (MQ=255)
                                                     gCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtgtg                                                >  1:2296144/1‑100 (MQ=255)
                                                     gCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAAtgtg                                                <  1:2648252/100‑1 (MQ=255)
                                                       tCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCa                                              <  1:1145082/100‑1 (MQ=255)
                                                          taattaatGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACt                                           <  2:657119/100‑1 (MQ=255)
                                                             ttaatGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTtct                                        >  2:2236244/1‑100 (MQ=255)
                                                               aatGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTcc                                      <  2:1171807/100‑1 (MQ=255)
                                                                atGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCa                                     <  1:410178/100‑1 (MQ=255)
                                                                atGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCa                                     <  2:1212891/100‑1 (MQ=255)
                                                                 tGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCaa                                    >  2:481015/1‑100 (MQ=255)
                                                                    cTAAAACGAGTAAAGTTTACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGACATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGa                                 >  1:3094065/1‑100 (MQ=255)
                                                                      aaaaCAAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATc                               <  2:352206/100‑1 (MQ=255)
                                                                       aaaCGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACAGATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCt                              <  2:2797506/100‑1 (MQ=255)
                                                                         aCGGGTAAAGTTTACCGCCGAGAATTGGGCTGTGAATAACGACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCATTGATc                               <  2:1631292/97‑1 (MQ=255)
                                                                           gAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGa                            <  2:2272952/98‑1 (MQ=255)
                                                                           gAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAg                          <  1:1069714/100‑1 (MQ=255)
                                                                            gtggttattCACCGCCCAATTTTCAGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGt                         >  2:576293/8‑100 (MQ=37)
                                                                             gTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGtt                        >  1:3390282/1‑100 (MQ=255)
                                                                             gTAAAGTTAAACGACGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGAAATATGTTTTGTCAAAATGTGCAACTTCTCAAATGATCTGAAGtt                        >  2:742884/1‑100 (MQ=255)
                                                                              tAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTg                       <  2:646710/100‑1 (MQ=255)
                                                                               gttattCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGa                      >  1:1361817/5‑100 (MQ=37)
                                                                               gttattCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGCCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGa                      >  2:2832663/5‑100 (MQ=37)
                                                                                ttattCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGaa                     >  2:406286/4‑100 (MQ=38)
                                                                                aaGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGaa                     <  2:1776521/100‑1 (MQ=255)
                                                                                  attCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTATCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAAc                   >  2:2282926/2‑100 (MQ=38)
                                                                                   ttCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACg                  >  1:3611230/1‑100 (MQ=38)
                                                                                   ttCACCGCCCAATTTTCGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACg                  >  1:704996/1‑100 (MQ=38)
                                                                                      aCCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGa               >  2:692889/1‑100 (MQ=255)
                                                                                         gccgAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAg            >  2:3195766/1‑100 (MQ=255)
                                                                                             aaaaTTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTGTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGcc          <  2:1478015/98‑1 (MQ=255)
                                                                                              aaaTTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTc       <  1:2622688/100‑1 (MQ=255)
                                                                                               aaTTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTca      <  1:1546957/100‑1 (MQ=255)
                                                                                                 ttGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTcaaa    <  2:1366176/100‑1 (MQ=255)
                                                                                                 ttGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTcaaa    <  1:914285/100‑1 (MQ=255)
                                                                                                  tGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTcaaac   <  2:458270/100‑1 (MQ=255)
                                                                                                   gggCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTcaaaca  <  2:3563905/100‑1 (MQ=255)
                                                                                                   |                                                                                                   
TTCACATCGACGCTTCCCGCCTTCAATCCATCCGTTGAATGAGTTTTTTTAAAGCTCGTAATTAATGGCTAAAACGAGTAAAGTTCACCGCCGAAAATTGGGCGGTGAATAACCACGTTTGAAATATTGTGACATATGTTTTGTCAAAATGTGCAACTTCTCCAATGATCTGAAGTTGAAACGTGATAGCCGTCAAACA  >  NC_000913/1157537‑1157735

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