Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,335 G→T 15.5% intergenic (+33/+27) 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,3350GT15.5% 330.6 / 41.8 129intergenic (+33/+27)ptsG/fhuEfused glucose‑specific PTS enzymes: IIB component/IIC component/ferric‑rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (60/49);  new base T (6/14);  total (66/63)
Fisher's exact test for biased strand distribution p-value = 5.17e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.29e-01

Only 100 of 130 total aligned reads displayed.

GGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGG  >  NC_000913/1159236‑1159432
                                                                                                   |                                                                                                 
ggCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGCGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGGCTAAGGCAGCCAg                                                                                                   >  2:2505712/1‑100 (MQ=255)
 gCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGa                                                                                                  <  2:208890/100‑1 (MQ=255)
   gATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATg                                                                                                >  1:1868510/1‑100 (MQ=255)
    aTTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATgg                                                                                               <  1:2006879/100‑1 (MQ=255)
      tttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCt                                                                                             >  2:1202229/1‑100 (MQ=255)
      tgtgggTACTAAATCCGGTAGCCTGAAAGTCGGGAGGGATGATTACAGCCGTGACCACTAATCCGTAAGGCGGTGGGGAGACGAAGGCAGCCATCTGGc                                                                                              <  1:2466785/95‑1 (MQ=255)
         cGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGcc                                                                                          <  2:2549967/100‑1 (MQ=255)
         cGGTACTAAATCCCATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGTGAGACTAAGGCAGCCAGATGGCTGcc                                                                                          <  1:1219747/100‑1 (MQ=255)
          ggTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCt                                                                                         <  2:2442332/100‑1 (MQ=255)
          ggTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCt                                                                                         >  2:923344/1‑100 (MQ=255)
           gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtt                                                                                        <  1:140190/100‑1 (MQ=255)
           gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtt                                                                                        <  2:1635634/100‑1 (MQ=255)
           gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtt                                                                                        <  2:2761663/100‑1 (MQ=255)
            tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGCCTAAGGCAGCCAGATGGCTGCCttt                                                                                       <  1:526075/100‑1 (MQ=255)
            tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGTCTGCCtta                                                                                       <  1:104292/100‑2 (MQ=255)
            tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                       <  2:911843/100‑2 (MQ=255)
            tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtt                                                                                        <  1:732554/99‑1 (MQ=255)
            tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttt                                                                                       <  1:2670324/100‑1 (MQ=255)
             attAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttt                                                                                      <  1:1764790/98‑1 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                      <  2:16577/100‑3 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                      <  2:2175591/100‑3 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                      <  2:651994/100‑3 (MQ=255)
              cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttagt                                                                                     <  2:384758/100‑4 (MQ=255)
              cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                     >  1:1088338/1‑100 (MQ=255)
              cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGAAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                     >  1:2905462/1‑100 (MQ=255)
               tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttttt                                                                                    <  2:690594/100‑1 (MQ=255)
               tAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCtttttt                                                                                    >  2:1729286/1‑100 (MQ=255)
                aaaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTa                                                                                   <  2:1092359/100‑1 (MQ=255)
                aaaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTa                                                                                   >  2:2810909/1‑100 (MQ=255)
                aaaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTa                                                                                   >  1:1673502/1‑100 (MQ=255)
                 aaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTCACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAc                                                                                  >  1:1591743/1‑100 (MQ=255)
                  aTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACa                                                                                 >  1:846431/1‑100 (MQ=255)
                    ccGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAgg                                                                               <  2:2845188/100‑1 (MQ=255)
                     cGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgt                                                                              >  1:2134278/1‑100 (MQ=255)
                      gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtg                                                                             <  2:1054226/100‑1 (MQ=255)
                      gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtg                                                                             <  2:2632263/100‑1 (MQ=255)
                      gATAACCGGAAAACCGAGATGGATGAGTACAGCCGTAACCACTAAGCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtg                                                                             <  1:181003/100‑1 (MQ=255)
                         aaCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTa                                                                          <  1:119144/100‑1 (MQ=255)
                          aCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAt                                                                         >  2:590493/1‑100 (MQ=255)
                          aCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAt                                                                         >  2:1292987/1‑100 (MQ=255)
                           ccTGAAAACCGAGATGGATTAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTa                                                                          <  1:2845706/98‑1 (MQ=255)
                           ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                        >  2:2109879/1‑100 (MQ=255)
                            cTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTc                                                                       >  2:896768/1‑100 (MQ=255)
                             tGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGGTGGGGAGCATAAAGCAGCCAGATGGCTGCGTTTTTTACAGgtgtga                                                                          >  1:1974226/1‑94 (MQ=255)
                                aaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGaa                                                                   <  2:427162/100‑1 (MQ=255)
                                 aaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAt                                                                  <  2:2216974/100‑1 (MQ=255)
                                  aCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAtt                                                                 <  2:1542786/100‑1 (MQ=255)
                                  aCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAtt                                                                 <  2:86086/100‑1 (MQ=255)
                                   ccGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTg                                                                <  1:2810891/100‑1 (MQ=255)
                                     gagaTGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGAt                                                              <  2:1673509/100‑1 (MQ=255)
                                     gagaTGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGAt                                                              <  1:1729278/100‑1 (MQ=255)
                                         tggatgAGTACATCCGTAACAACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACTGGTGTTCTTCAGAATTGATCCGt                                                          >  1:2659848/1‑100 (MQ=255)
                                          ggatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTg                                                         >  1:218173/1‑100 (MQ=255)
                                           gatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGc                                                        >  2:270962/1‑100 (MQ=255)
                                            atgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGcc                                                       <  2:702028/100‑1 (MQ=255)
                                              gAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCgg                                                     <  1:460870/100‑1 (MQ=255)
                                               aGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGt                                                    >  1:2300286/1‑100 (MQ=255)
                                               aGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGt                                                    <  2:567941/100‑1 (MQ=255)
                                                 tACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTaa                                                  <  2:1263158/100‑1 (MQ=255)
                                                  aCATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAAt                                                 <  2:1465121/100‑1 (MQ=255)
                                                   cATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATg                                                >  1:1708927/1‑100 (MQ=255)
                                                    aTCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGc                                               <  2:2266262/100‑1 (MQ=255)
                                                        gTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGaa                                           <  2:2034252/100‑1 (MQ=255)
                                                         tAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGaaa                                          >  2:2334427/1‑100 (MQ=255)
                                                         tAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGaaa                                          >  1:2709226/1‑100 (MQ=255)
                                                           aCCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAAtt                                        >  2:2246149/1‑100 (MQ=255)
                                                           aCCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAAtt                                        >  1:1187288/1‑100 (MQ=255)
                                                              aCTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcg                                     <  2:1088343/100‑1 (MQ=255)
                                                               cTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGATTTGATACGTGCCGGTAATGCTGAAATTa                                       >  1:2382170/1‑97 (MQ=255)
                                                               cTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcgc                                    >  2:946452/1‑100 (MQ=255)
                                                                 aaTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCgg                                  <  2:686137/100‑1 (MQ=255)
                                                                 aaTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTACAGGTAATACTGAAATTa                                       >  1:2325589/1‑95 (MQ=255)
                                                                    ccGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGgtgt                               >  2:117822/1‑100 (MQ=255)
                                                                     cGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGtgtg                              <  2:647262/100‑1 (MQ=255)
                                                                      gTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGc                             >  2:2107319/1‑100 (MQ=255)
                                                                       tAAGACGTTTGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                            <  1:2576244/100‑1 (MQ=255)
                                                                        aaGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCg                           >  1:2170910/1‑100 (MQ=255)
                                                                            cGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGa                       >  1:352377/1‑100 (MQ=255)
                                                                             gTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGan                      <  1:347568/100‑2 (MQ=255)
                                                                               tGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGa                    >  2:422534/1‑100 (MQ=255)
                                                                                ggggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGAt                   >  2:1096551/1‑100 (MQ=255)
                                                                                  tgtaaaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAg                 >  2:1046170/8‑100 (MQ=38)
                                                                                  ggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAg                 <  1:2623629/100‑1 (MQ=255)
                                                                                  ggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAg                 <  1:1202224/100‑1 (MQ=255)
                                                                                     aaaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGaa               >  2:2764483/5‑99 (MQ=39)
                                                                                       cTAAGGCAGCCAGATGGCTGCCTTTTTTTTTTTTTTTTTTAACAAAAGACACCTGCCGGTAATTCTCCAATTCCCCGGTGTACACTAGACGacatccc              >  1:268338/1‑92 (MQ=255)
                                                                                       cTAAGGCAGCCAGATGGCGGCCTACTTTACAGGTGTTATGCAGAATTGATACGCGCCGGTAATGCGGAAATTACGTGGTGTGCCGGAGACGATAGAAc              >  1:2117328/1‑98 (MQ=255)
                                                                                        tAAGGCAGCCAGATGGCTGACTTTTTTAAAGGTGGAATTCAGACTTGATACGGTCCTGTAATGTTGAAATTACGGGGTGTGAGGTAGACGATAGAAc              >  1:2011154/1‑97 (MQ=255)
                                                                                           ggCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa        >  1:328526/1‑100 (MQ=255)
                                                                                           ggCAGCCAGATGGCTGCCTTTTTTACAAGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa        >  1:328488/1‑100 (MQ=255)
                                                                                            gCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCAc       <  2:2686812/100‑1 (MQ=255)
                                                                                             ctgCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACg      >  2:1979008/3‑100 (MQ=255)
                                                                                              aGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGt     >  1:186744/1‑100 (MQ=255)
                                                                                               gCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACg      <  1:987150/98‑1 (MQ=255)
                                                                                               gCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGtt    >  1:2575042/1‑100 (MQ=255)
                                                                                               gCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGtt    <  1:740563/100‑1 (MQ=255)
                                                                                                 cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg  >  1:1172788/1‑100 (MQ=255)
                                                                                                 cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg  >  1:1130315/1‑100 (MQ=255)
                                                                                                 cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg  <  2:1550076/100‑1 (MQ=255)
                                                                                                 cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTgg  <  2:965859/100‑1 (MQ=255)
                                                                                                   |                                                                                                 
GGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGG  >  NC_000913/1159236‑1159432

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: