Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,159,335 G→T 22.2% 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,3350GT22.2% 372.0 / 94.0 162intergenic (+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 (64/62);  new base T (13/23);  total (77/85)
Fisher's exact test for biased strand distribution p-value = 1.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

Only 100 of 166 total aligned reads displayed.

GGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTA  >  NC_000913/1159236‑1159434
                                                                                                   |                                                                                                   
ggCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAg                                                                                                     >  1:1614727/1‑100 (MQ=255)
  cGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGAt                                                                                                   >  2:909147/1‑100 (MQ=255)
   gATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATg                                                                                                  <  1:4120281/100‑1 (MQ=255)
    aTTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATgg                                                                                                 >  2:238557/1‑100 (MQ=255)
    aTTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATgg                                                                                                 <  1:1475584/100‑1 (MQ=255)
     ttttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGc                                                                                                >  1:593009/1‑100 (MQ=255)
       ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                              <  1:294849/100‑1 (MQ=255)
       ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                              <  2:4027855/100‑1 (MQ=255)
       ttCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                              <  1:3539134/100‑1 (MQ=255)
         cGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTg                                                                                              <  1:3181539/98‑1 (MQ=255)
         cGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGcc                                                                                            <  2:2926897/100‑1 (MQ=255)
         cGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGcc                                                                                            >  1:3462069/1‑100 (MQ=255)
           gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGccct                                                                                          <  1:2820362/100‑3 (MQ=255)
           gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtt                                                                                          <  1:1683751/100‑1 (MQ=255)
           gTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtt                                                                                          <  1:2329835/100‑1 (MQ=255)
            tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                         <  2:2099235/100‑2 (MQ=255)
            tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                         <  2:2570323/100‑2 (MQ=255)
            tACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCtta                                                                                         <  2:2750579/100‑2 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                        <  1:4153216/100‑3 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                        <  1:358393/100‑3 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                        <  1:3437278/100‑3 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                        <  2:76612/100‑3 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                        <  1:1040486/100‑3 (MQ=255)
             aCTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCATCTGGCTGCCttag                                                                                        <  1:1061390/100‑3 (MQ=255)
              cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                       >  2:3933293/1‑100 (MQ=255)
              cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                       >  2:1148341/1‑100 (MQ=255)
              cTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCttttt                                                                                       >  2:2678925/1‑100 (MQ=255)
                 aaTCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGCTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTAc                                                                                    <  1:1227819/100‑1 (MQ=255)
                   tCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAg                                                                                  <  2:524890/100‑1 (MQ=255)
                     cGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgt                                                                                >  1:665111/1‑100 (MQ=255)
                      gATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtg                                                                               <  1:1858641/100‑1 (MQ=255)
                       aTAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGgtg                                                                               <  2:1848285/99‑1 (MQ=255)
                        tAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGtt                                                                             <  1:1883545/100‑1 (MQ=255)
                          aCCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAt                                                                           >  2:2261422/1‑100 (MQ=255)
                           ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                          >  2:3976745/1‑100 (MQ=255)
                           ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                          >  1:4104031/1‑100 (MQ=255)
                           ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                          >  1:3092197/1‑100 (MQ=255)
                           ccTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTAtt                                                                          >  2:2071277/1‑100 (MQ=255)
                            cTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTc                                                                         <  2:560118/100‑1 (MQ=255)
                             tGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCa                                                                        >  1:1507663/1‑100 (MQ=255)
                                aaaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGaa                                                                     >  1:1260702/1‑100 (MQ=255)
                                 aaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAt                                                                    <  1:3894756/100‑1 (MQ=255)
                                 aaCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAt                                                                    <  1:332617/100‑1 (MQ=255)
                                  aCCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAAtt                                                                   <  1:1080595/100‑1 (MQ=255)
                                   ccGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTg                                                                  >  1:3344320/1‑100 (MQ=255)
                                     gagaTGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGAt                                                                <  1:3203030/100‑1 (MQ=255)
                                       gatgtatgAGTACACCCGTAACCACTAATCCGTCAGACTTTGTGGAGACTAAGGCTGCTAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATa                                                               <  2:556375/99‑1 (MQ=255)
                                        atggatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACg                                                             <  1:790718/100‑1 (MQ=255)
                                        atggatgAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACg                                                             <  1:638190/100‑1 (MQ=255)
                                              gAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATAGATACg                                                             >  2:335619/1‑94 (MQ=255)
                                               aGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGt                                                      >  1:3297201/1‑100 (MQ=255)
                                                 tACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGt                                                      <  2:3439461/98‑1 (MQ=255)
                                                   cATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATg                                                  >  1:3827669/1‑100 (MQ=255)
                                                    aTCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGc                                                 <  1:3183425/100‑1 (MQ=255)
                                                    aTCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGc                                                 <  1:3858251/100‑1 (MQ=255)
                                                      ccGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTg                                               <  2:267854/100‑1 (MQ=255)
                                                            ccACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTa                                         >  2:3125799/1‑100 (MQ=255)
                                                             cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAc                                        >  2:3774066/1‑100 (MQ=255)
                                                             cACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAc                                        >  1:3241535/1‑100 (MQ=255)
                                                               cTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTAcg                                       >  2:2354110/1‑99 (MQ=255)
                                                                 aaTCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCgg                                    >  1:1768924/1‑100 (MQ=255)
                                                                     cGTAAGACGTTGGGGAGACTAAGGTAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGCATTGATACGTGCCGGTAATGCTGAAATTACGCGGtgtg                                <  2:2583150/100‑1 (MQ=255)
                                                                     cGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGtgtg                                >  1:1364543/1‑100 (MQ=255)
                                                                       tAAGACGTTGGGGAGACTAAGGCAGCTAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                              >  1:2544441/1‑100 (MQ=255)
                                                                       tAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGTTAATGCTGAAATTACGCGGTGTGc                               <  2:2999760/99‑1 (MQ=255)
                                                                       tAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGcc                              <  1:2123738/100‑1 (MQ=255)
                                                                        aaGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACCGGAGTGATTCCGAATTGATACGCGCCGGTAATTCTGAAATTACGCGTTTTGc                               >  2:3858636/1‑98 (MQ=255)
                                                                         aGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGt                            <  1:2489324/100‑1 (MQ=255)
                                                                          gACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTa                           <  1:3105346/100‑1 (MQ=255)
                                                                           aCGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAg                          <  1:314240/100‑1 (MQ=255)
                                                                           aCGTTGGGGAGACTAAGGCAGCCAGATGGCTGACTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAg                          <  1:2066674/100‑1 (MQ=255)
                                                                            cGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGa                         >  2:167307/1‑100 (MQ=255)
                                                                                ggggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGAt                     <  1:204711/100‑1 (MQ=255)
                                                                                ggggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGAt                     >  1:1062431/1‑100 (MQ=255)
                                                                                ggggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGAt                     <  1:1844752/100‑1 (MQ=255)
                                                                                  ggAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAg                   <  1:3013995/100‑1 (MQ=255)
                                                                                   gagaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGACGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGa                  <  2:3732108/100‑1 (MQ=255)
                                                                                    taaaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCAGTGTGCCGTAGAGGATAGaa                 >  2:1553917/6‑100 (MQ=38)
                                                                                    agaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGaa                 <  1:484377/100‑1 (MQ=255)
                                                                                     gaCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGaa                 >  2:2646972/1‑99 (MQ=255)
                                                                                     aaaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc                >  2:83823/5‑100 (MQ=39)
                                                                                     aaaaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAAc                >  1:937996/5‑100 (MQ=255)
                                                                                      aCTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGCATTGATACGTGCCGGTAATGCTGCAATTACGCGGTGTGCCGTACCCGATa                    >  2:2954997/1‑95 (MQ=255)
                                                                                       cTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt              <  1:1010953/100‑1 (MQ=255)
                                                                                       cTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCt              >  2:19930/1‑100 (MQ=255)
                                                                                       aaaaGGCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGAGGTGTGCCGTAGACGATAGAACCt              >  2:1028287/3‑100 (MQ=255)
                                                                                        tAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCtt             >  2:1844054/1‑100 (MQ=255)
                                                                                        tAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCtt             >  1:937205/1‑100 (MQ=255)
                                                                                         aaGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc            <  2:4013457/100‑1 (MQ=255)
                                                                                         aaGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTc            >  1:395746/1‑100 (MQ=255)
                                                                                          aGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTcc           >  1:2960416/1‑100 (MQ=255)
                                                                                           ggCAGCCATCTGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa          >  2:3162934/1‑100 (MQ=255)
                                                                                           ggCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCTGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa          >  2:3592638/1‑100 (MQ=255)
                                                                                           ggCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCa          >  2:3105525/1‑100 (MQ=255)
                                                                                             cAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACg        >  2:1314250/1‑100 (MQ=255)
                                                                                               gCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGtt      >  2:2501442/1‑100 (MQ=255)
                                                                                               gCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGtt      >  2:1725222/1‑100 (MQ=255)
                                                                                                 cAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTg     <  2:697928/99‑1 (MQ=255)
                                                                                                   gATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTa  >  2:3926054/1‑100 (MQ=255)
                                                                                                   gATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTa  >  1:3233138/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GGCGATTTTCGGTACTAAATCCGATAACCTGAAAACCGAGATGGATGAGTACATCCGTAACCACTAATCCGTAAGACGTTGGGGAGACTAAGGCAGCCAGATGGCTGCCTTTTTTACAGGTGTTATTCAGAATTGATACGTGCCGGTAATGCTGAAATTACGCGGTGTGCCGTAGACGATAGAACCTTCCACGTTGGTA  >  NC_000913/1159236‑1159434

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: