Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,336 A→C 16.0% intergenic (+34/+26) ptsG → / ← fhuE fused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,3360AC16.0% 366.9 / 65.7 155intergenic (+34/+26)ptsG/fhuEfused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (56/74);  new base C (12/13);  total (68/87)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

Only 100 of 157 total aligned reads displayed.

CGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTAT  >  NC_000913/1159238‑1159435
                                                                                                  |                                                                                                   
cGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGa                                                                                                     >  1:415596/1‑99 (MQ=255)
     ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                               <  1:3518983/100‑1 (MQ=255)
     ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                               <  2:823183/100‑1 (MQ=255)
     ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                               >  2:3035053/1‑100 (MQ=255)
       cGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGcc                                                                                             <  2:2218328/100‑1 (MQ=255)
        ggTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCt                                                                                            <  1:3246542/100‑1 (MQ=255)
        ggTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCt                                                                                            <  1:2183077/100‑1 (MQ=255)
        ggTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCt                                                                                            >  2:2896910/1‑100 (MQ=255)
         ggaCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtt                                                                                           <  1:3427122/98‑1 (MQ=255)
         gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTATGACGTAGGGGAGACTAAGGCAGTCAGATGGCTGCCt                                                                                            >  1:2767481/1‑99 (MQ=255)
          tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                          <  2:2549663/100‑2 (MQ=255)
          tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                          <  2:1062058/100‑2 (MQ=255)
          tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                          <  2:3070882/100‑2 (MQ=255)
          tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                          <  1:1343020/100‑2 (MQ=255)
           aCTAAATCCTATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                         <  2:1399290/100‑3 (MQ=255)
           aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                         <  1:2047877/100‑3 (MQ=255)
           aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttt                                                                                         >  2:3270940/1‑100 (MQ=255)
           aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttt                                                                                         >  1:1344722/1‑100 (MQ=255)
            cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttagt                                                                                        <  2:718929/100‑4 (MQ=255)
            cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                        <  2:2241462/100‑1 (MQ=255)
            cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                        >  2:565347/1‑100 (MQ=255)
             tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttttt                                                                                       <  1:84541/100‑1 (MQ=255)
             tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttttt                                                                                       <  2:1891024/100‑1 (MQ=255)
             tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCAAGATGACTGCCtttttt                                                                                       >  1:1244106/1‑100 (MQ=255)
              aaaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTa                                                                                      <  1:2046709/100‑1 (MQ=255)
                  ccGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAgg                                                                                  <  2:1857773/100‑1 (MQ=255)
                    gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGCAACTAAGGCAGCCAGATGGCTGCCTTTTTTACAg                                                                                   <  1:976472/97‑1 (MQ=255)
                      tAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGtt                                                                              <  1:2795377/100‑1 (MQ=255)
                       aaCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTa                                                                             <  1:1774321/100‑1 (MQ=255)
                         ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                           >  2:1119509/1‑100 (MQ=255)
                         ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                           >  2:3453093/1‑100 (MQ=255)
                         ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                           <  2:932483/100‑1 (MQ=255)
                          cTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTc                                                                          <  1:2977348/100‑1 (MQ=255)
                            gAAAACCGAGATGGATTATGACATCCGTAACCACTAATCCGTAAGACGTTGGGGATACTAAGGCAGCCACATTGCTGCCTTTTTTACAGGTGCTATTCAg                                                                        <  1:3120899/100‑1 (MQ=255)
                             aaaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGa                                                                       <  2:2001083/100‑1 (MQ=255)
                              naaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGATGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGaa                                                                      >  1:1965088/2‑100 (MQ=255)
                              aaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGaa                                                                      <  2:959626/100‑1 (MQ=255)
                               ttccgatctggatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAt                                                                     <  2:2292446/92‑1 (MQ=255)
                               aaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAt                                                                     >  2:2787308/1‑100 (MQ=255)
                                    ttccgatctaGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATa                                                                <  2:1224691/91‑1 (MQ=255)
                                    agaTGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATa                                                                >  2:3269114/1‑100 (MQ=255)
                                      tccgatctgTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACg                                                              <  1:1576346/92‑1 (MQ=255)
                                      atggatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACg                                                              <  2:815840/100‑1 (MQ=255)
                                        ggatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATAGAtcggaa                                                            >  1:2292432/1‑95 (MQ=255)
                                          atgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGcc                                                          <  2:2290566/100‑1 (MQ=255)
                                            gAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCgg                                                        <  2:752519/100‑1 (MQ=255)
                                             aGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATaagatcggaa                                                       >  1:1224684/1‑91 (MQ=255)
                                                aCATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAAt                                                    >  2:1467960/1‑100 (MQ=255)
                                                  aTCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCCGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGc                                                  >  2:1392520/1‑100 (MQ=255)
                                                   tCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCt                                                 <  1:3292905/100‑1 (MQ=255)
                                                     cGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGa                                               <  2:2943706/100‑1 (MQ=255)
                                                       tAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGaaa                                             <  2:3475125/100‑1 (MQ=255)
                                                        aaCCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAAt                                            <  2:65115/100‑1 (MQ=255)
                                                         aCCACTAATCCTTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTTTTCAGAATTGATACTTGCCGGTAATGCTGAAAt                                            >  1:1121136/1‑99 (MQ=255)
                                                          ccACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTa                                          <  1:223840/100‑1 (MQ=255)
                                                          ccACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTa                                          >  1:310314/1‑100 (MQ=255)
                                                          ccACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTa                                          >  2:2880933/1‑100 (MQ=255)
                                                           cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAc                                         >  2:1168808/1‑100 (MQ=255)
                                                           cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAc                                         >  1:2572716/1‑100 (MQ=255)
                                                             cTAATCCGTAAGACGTTGGGGAGACTAAGGGAGACAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCCGAAATTAcgc                                       >  1:2339117/1‑100 (MQ=255)
                                                             cTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcgc                                       <  2:27847/100‑1 (MQ=255)
                                                             cTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcgc                                       >  2:1613112/1‑100 (MQ=255)
                                                             cTAATCCGTAAGAAGATGGCGAGACTAAGGCAGCCAGCTGGCAGCCTTGTTAACAGGTGTTCTGCAAAATTGATACATGCCGGCAACGCCGAAATTAcgc                                       >  1:3522966/1‑100 (MQ=255)
                                                               aaTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCgg                                     >  1:448327/1‑100 (MQ=255)
                                                                aTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGgt                                    <  1:3352056/100‑1 (MQ=255)
                                                                  ccGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGgtgt                                  <  1:1613679/100‑1 (MQ=255)
                                                                   cGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGtgtg                                 >  1:2523650/1‑100 (MQ=255)
                                                                      aaGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCg                              <  2:2097409/100‑1 (MQ=255)
                                                                        gACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                            <  2:2079473/100‑1 (MQ=255)
                                                                        gACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                            >  1:1831979/1‑100 (MQ=255)
                                                                         aCGTTGGGGAGGCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGGAATACTGAAATTACGCGGAGTTCCg                              >  1:3103894/1‑97 (MQ=255)
                                                                           gTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGAc                         <  2:2321070/100‑1 (MQ=255)
                                                                           gTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGAc                         <  1:2449503/100‑1 (MQ=255)
                                                                           gTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGAc                         <  2:64741/100‑1 (MQ=255)
                                                                           gTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGAc                         <  2:178517/100‑1 (MQ=255)
                                                                            ttGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGAc                         >  1:219197/1‑99 (MQ=255)
                                                                            ttGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACg                        <  2:3314121/100‑1 (MQ=255)
                                                                               gggAGACTAAGGCAGCCAGATGGCTGCCTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATa                     <  1:1701399/99‑1 (MQ=255)
                                                                                tgtaacaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAg                    >  1:2944956/8‑100 (MQ=39)
                                                                                   gaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc                 <  2:2237836/100‑1 (MQ=255)
                                                                                   aaaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc                 >  2:3153738/5‑100 (MQ=39)
                                                                                    aCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCTGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGa                   <  1:672196/97‑1 (MQ=255)
                                                                                     ntAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt               >  1:2201402/2‑100 (MQ=255)
                                                                                     cTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt               >  1:2215692/1‑100 (MQ=255)
                                                                                     aaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt               >  2:1615179/3‑100 (MQ=255)
                                                                                      tAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCtt              <  2:2013214/100‑1 (MQ=255)
                                                                                      aaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCtt              >  2:1678949/2‑100 (MQ=255)
                                                                                       aaGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc             >  2:438068/1‑100 (MQ=255)
                                                                                        aGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTcc            <  2:3191556/100‑1 (MQ=255)
                                                                                        aGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTcc            <  1:1159047/100‑1 (MQ=255)
                                                                                         ggCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa           >  1:3245449/1‑100 (MQ=255)
                                                                                         ggCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa           >  2:896740/1‑100 (MQ=255)
                                                                                         ggCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa           <  1:2422502/100‑1 (MQ=255)
                                                                                           cAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACg         >  2:1534150/1‑100 (MQ=255)
                                                                                             gCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGtt       >  1:2045721/1‑100 (MQ=255)
                                                                                             gCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGtt       >  2:722790/1‑100 (MQ=255)
                                                                                              ccAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTg      >  1:1961741/1‑100 (MQ=255)
                                                                                               cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg     <  2:1831988/100‑1 (MQ=255)
                                                                                               cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg     <  2:2572733/100‑1 (MQ=255)
                                                                                                  aTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTAt  <  1:896734/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTAT  >  NC_000913/1159238‑1159435

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