Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 851,094 G→A 67.1% P502S (CCG→TCG)  opgE ← OPG biosynthetic transmembrane phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913851,0940GA67.1% 118.8 / 86.5 109P502S (CCG→TCG) opgEOPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base G (19/17);  new base A (29/44);  total (48/61)
Fisher's exact test for biased strand distribution p-value = 2.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 111 total aligned reads displayed.

GGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTG  >  NC_000913/850997‑851192
                                                                                                 |                                                                                                  
ggggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGAcc                                                                                                  <  2:1363098/100‑1 (MQ=255)
ggggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTCGAAGTTGGTGATGTATATCGAcc                                                                                                  >  1:2864949/1‑100 (MQ=255)
 gggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGccc                                                                                                 <  1:1823246/100‑1 (MQ=255)
  ggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCg                                                                                                >  2:3225677/1‑100 (MQ=255)
  ggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCg                                                                                                >  1:35639/1‑100 (MQ=255)
  ggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCg                                                                                                >  2:2379030/1‑100 (MQ=255)
  ggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCg                                                                                                >  1:1604148/1‑100 (MQ=255)
  ggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCg                                                                                                >  1:1669258/1‑100 (MQ=255)
    cGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGc                                                                                              <  2:1796092/100‑1 (MQ=255)
      gATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCtt                                                                                            <  2:1712691/100‑1 (MQ=255)
        tttttGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTttc                                                                                          <  2:802460/100‑1 (MQ=255)
          tttGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTttctt                                                                                        >  1:2939723/1‑100 (MQ=255)
             gTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTtctc                                                                                     <  1:620220/100‑1 (MQ=255)
             gTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCATCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTtctc                                                                                     <  1:2039687/100‑1 (MQ=255)
             gTTTTTAAGGTTTCGGGTCCCAAATATCGGTTCCGATATGGATGTAATCCACCTTCTGTAACTAGCACTTGGTAGTGTATATCGGCCCGGCTTacttct                                                                                      >  1:3055483/1‑99 (MQ=255)
               ttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAg                                                                                   <  1:1913739/100‑1 (MQ=255)
                tttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGa                                                                                  <  2:1924718/100‑1 (MQ=255)
                tttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGa                                                                                  <  2:484718/100‑1 (MQ=255)
                tttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGa                                                                                  >  2:3070210/1‑100 (MQ=255)
                tttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGa                                                                                  >  2:3180668/1‑100 (MQ=255)
                tttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGa                                                                                  <  2:3208424/100‑1 (MQ=255)
                  tAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGGAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGata                                                                                >  2:481368/1‑100 (MQ=255)
                  tAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGata                                                                                <  1:2496586/100‑1 (MQ=255)
                   aaGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGatat                                                                               <  2:3392496/100‑1 (MQ=255)
                     ggTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGCAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATaa                                                                             >  1:2163450/1‑100 (MQ=255)
                        ttCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACg                                                                          <  2:117328/100‑1 (MQ=255)
                          cGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACggg                                                                        <  2:2630/100‑1 (MQ=255)
                          cGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACggg                                                                        <  2:2243348/100‑1 (MQ=255)
                            ggTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTa                                                                      >  1:163804/1‑100 (MQ=255)
                             gTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTAt                                                                     >  2:269414/1‑100 (MQ=255)
                                gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                                  >  2:2003060/1‑100 (MQ=255)
                                gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                                  >  1:3056805/1‑100 (MQ=255)
                                gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATACCGACCGGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                                  >  1:2286958/1‑100 (MQ=255)
                                gAAAATATCGGTTACGACAAGGATGTAATCCACCTTCTGTAACTAGAAGTTGGTGATGTATTTCGACTCAGGTTTCTTCTCAGAGATAAATGAGTaaat                                                                   >  1:3155161/1‑96 (MQ=255)
                                 aaaaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGa                                                                 <  1:2293176/100‑1 (MQ=255)
                                  aaaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTATCTAGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGAt                                                                <  1:2379013/100‑1 (MQ=255)
                                   aaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg                                                               <  1:3193217/100‑1 (MQ=255)
                                   aaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg                                                               <  2:2970722/100‑1 (MQ=255)
                                      atCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTa                                                            <  2:2893781/100‑1 (MQ=255)
                                        cGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAt                                                          <  1:3261434/100‑1 (MQ=255)
                                          gTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAttt                                                        <  2:2534145/100‑1 (MQ=255)
                                           ttCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTc                                                       <  2:2507359/100‑1 (MQ=255)
                                           ttCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTc                                                       <  1:1518614/100‑1 (MQ=255)
                                           ttCCGAGATGGTTGTAATCCACCTTATGTAACTGGAAGTTGGTGATGTATATCGACACGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTc                                                       >  1:1927574/1‑100 (MQ=255)
                                             cagagaTGGTTGGACTCAAACTTCTGTAACTGGAAGTGGGTGATGGATATCGGCCAGGCTTTCTTCTCATATATAAACGGGTATTTTATGTTAATTTc                                                       <  1:216984/96‑1 (MQ=255)
                                              ngagaTGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCtt                                                    >  2:908497/2‑100 (MQ=255)
                                              cGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCtt                                                    >  2:334865/1‑100 (MQ=255)
                                              cGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCtt                                                    <  2:1434033/100‑1 (MQ=255)
                                               gagaTGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTAGATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTc                                                   <  1:2669847/100‑1 (MQ=255)
                                                   tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                               >  2:2513644/1‑100 (MQ=255)
                                                   tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                               >  1:980163/1‑100 (MQ=255)
                                                   tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                               <  2:2268701/100‑1 (MQ=255)
                                                   tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                               <  2:829595/100‑1 (MQ=255)
                                                     gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                               >  1:2184854/1‑98 (MQ=255)
                                                     gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                             <  2:1011016/100‑1 (MQ=255)
                                                     gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                             <  2:138299/100‑1 (MQ=255)
                                                     gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                             <  1:108739/100‑1 (MQ=255)
                                                      ttGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTaa                                            <  2:276280/100‑1 (MQ=255)
                                                      ttGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTaa                                            >  2:3309757/1‑100 (MQ=255)
                                                       tGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAAt                                           >  1:1877639/1‑100 (MQ=255)
                                                        gTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATc                                          <  2:3358152/100‑1 (MQ=255)
                                                         tAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcc                                         >  2:469413/1‑100 (MQ=255)
                                                         tAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcc                                         <  2:1787036/100‑1 (MQ=255)
                                                         tAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATc                                          <  1:1249441/99‑1 (MQ=255)
                                                          aaTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATccc                                        >  1:1923273/1‑100 (MQ=255)
                                                           aTCCACCTTCTGTAACTGGAAGTTGGTTATGTATATCGACCCCGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcc                                         >  1:1561039/1‑98 (MQ=255)
                                                           aTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcccc                                       >  1:2607527/1‑100 (MQ=255)
                                                           aTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcccc                                       >  2:488045/1‑100 (MQ=255)
                                                           aTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcccc                                       >  1:3067313/1‑100 (MQ=255)
                                                            tCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCg                                      >  2:2238718/1‑100 (MQ=255)
                                                            tCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCg                                      >  2:463733/1‑100 (MQ=255)
                                                            tCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCg                                      <  1:458147/100‑1 (MQ=255)
                                                               aCCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTcc                                   <  2:2029873/100‑1 (MQ=255)
                                                               aCCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTcc                                   <  1:44238/100‑1 (MQ=255)
                                                                  ttCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACt                                >  2:3277075/1‑100 (MQ=255)
                                                                    cTGTATCTGGAATTTGGTGATGTGTAGCGGCCGGGATTTCTTCGCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTgt                              <  1:896240/100‑2 (MQ=255)
                                                                     tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTgg                              >  1:2984105/1‑99 (MQ=255)
                                                                      gTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGga                            >  1:369878/1‑100 (MQ=255)
                                                                        aaCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgaga                          <  1:587688/100‑1 (MQ=255)
                                                                          cTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaaa                        >  2:3066833/1‑100 (MQ=255)
                                                                          cTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaaa                        >  2:3232680/1‑100 (MQ=255)
                                                                           tGGAAGTTGGTGTTGTATATCTGCCCGGCTTTCTTCTCAGATATAAGCGTGTAGTTGATGTTAATTTCCTTCGCGTTAATCCCCTTCCACTGGga                            <  1:2015815/95‑1 (MQ=255)
                                                                              aaGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGc                    <  1:1314297/100‑1 (MQ=255)
                                                                                gTTGGGGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGcc                   <  1:953700/99‑1 (MQ=255)
                                                                                  tGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaa                >  2:1597956/1‑100 (MQ=255)
                                                                                      gATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaagaag            >  2:2181402/1‑100 (MQ=255)
                                                                                      gATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaagaag            <  2:867938/100‑1 (MQ=255)
                                                                                      gATGTATATCGACCCGGCTTTCTTCTCAGATATAACCGGGTCTTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaacaag            >  2:1559822/1‑100 (MQ=255)
                                                                                       aTGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGt           <  1:1453427/100‑1 (MQ=255)
                                                                                       aTGTAGATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTGTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGt           <  1:1559815/100‑1 (MQ=255)
                                                                                          tataTCGGCCCGGCTTTCTTCCCAGATATAAACGGGTATTTGATGTTAATTTCTTTCGCTTTAATCCCCGTCCACTGGGAGAAACAGCCTACGAAGTCAt        >  1:2927587/1‑100 (MQ=255)
                                                                                          tataTCGATCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCa         <  1:2820034/99‑1 (MQ=255)
                                                                                          tataTCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCAt        >  2:756538/1‑100 (MQ=255)
                                                                                            tatCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTg      <  2:150641/100‑1 (MQ=255)
                                                                                            tatCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTg      <  2:102861/100‑1 (MQ=255)
                                                                                             atCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTgg     <  2:2247749/100‑1 (MQ=255)
                                                                                             atCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTgg     <  1:2717596/100‑1 (MQ=255)
                                                                                              tCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGc    >  1:276241/1‑100 (MQ=255)
                                                                                              tCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGc    <  1:2567965/100‑1 (MQ=255)
                                                                                                ggCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTg  >  2:3247560/1‑100 (MQ=255)
                                                                                                 |                                                                                                  
GGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGCTG  >  NC_000913/850997‑851192

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: