Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 851,094 G→A 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,0940GA100.0% 310.4 / NA 92P502S (CCG→TCG) opgEOPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (48/44);  total (48/44)

CGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGC  >  NC_000913/851001‑851190
                                                                                             |                                                                                                
cGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGc                                                                                            <  2:524877/100‑1 (MQ=255)
  gATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCtt                                                                                          >  2:1083774/1‑100 (MQ=255)
  gATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCtt                                                                                          <  2:821461/100‑1 (MQ=255)
  gATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCtt                                                                                          <  2:764567/100‑1 (MQ=255)
     ttttGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTttct                                                                                       <  1:207491/100‑1 (MQ=255)
        tGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTtcttct                                                                                    >  1:1167306/1‑100 (MQ=255)
        tGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTtcttct                                                                                    <  2:3379256/100‑1 (MQ=255)
         gTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTtctc                                                                                   >  1:1572516/1‑100 (MQ=255)
         gTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTtctc                                                                                   <  1:2236904/100‑1 (MQ=255)
          tttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCa                                                                                  >  2:273143/1‑100 (MQ=255)
          tttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCa                                                                                  <  1:1910375/100‑1 (MQ=255)
           ttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAg                                                                                 <  1:1630666/100‑1 (MQ=255)
           ttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAg                                                                                 <  2:2704065/100‑1 (MQ=255)
            tttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGa                                                                                >  2:1738979/1‑100 (MQ=255)
            tttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGa                                                                                >  2:2945963/1‑100 (MQ=255)
            tttAAGGTTTCGGGTCGAAAATAACGGTTCCGAGATGGTTGTAATCCAGCTGCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCa                                                                                  >  1:837389/1‑98 (MQ=255)
             ttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGat                                                                               >  1:1119959/1‑100 (MQ=255)
               aaGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGatat                                                                             <  2:66051/100‑1 (MQ=255)
                aGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGAtata                                                                            <  1:1452335/100‑1 (MQ=255)
                   tttCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAAc                                                                         <  1:253002/100‑1 (MQ=255)
                    ttCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGGTATAAACg                                                                        >  1:2358237/1‑100 (MQ=255)
                      cGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACggg                                                                      <  1:1942771/100‑1 (MQ=255)
                        ggTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTa                                                                    >  2:219906/1‑100 (MQ=255)
                         gTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTAt                                                                   >  2:656542/1‑100 (MQ=255)
                          tCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTAtt                                                                  <  1:1346312/100‑1 (MQ=255)
                           cGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACg                                                                        >  1:2983426/1‑93 (MQ=255)
                           cGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTAttt                                                                 >  2:372402/1‑100 (MQ=255)
                           cGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTAttt                                                                 >  2:2838561/1‑100 (MQ=255)
                            gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                                <  2:417303/100‑1 (MQ=255)
                            gAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTg                                                                >  2:639221/1‑100 (MQ=255)
                               aaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg                                                             <  2:1342060/100‑1 (MQ=255)
                               aaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg                                                             <  2:2214846/100‑1 (MQ=255)
                               aaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg                                                             <  2:1915778/100‑1 (MQ=255)
                                atatCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGGTGTATATCGACTCGGCTTTCTTGTCAGATATAAAGGGGGATTTGCTGt                                                            >  1:2533133/1‑100 (MQ=255)
                                 tatCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGtt                                                           >  2:1512037/1‑100 (MQ=255)
                                  atCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTa                                                          <  1:999149/100‑1 (MQ=255)
                                  atCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTa                                                          >  1:242703/1‑100 (MQ=255)
                                    cGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAt                                                        <  1:1673395/100‑1 (MQ=255)
                                     ggTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAtt                                                       >  2:3099092/1‑100 (MQ=255)
                                      gTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAttt                                                      <  2:861580/100‑1 (MQ=255)
                                          cGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCtt                                                  <  2:1890108/100‑1 (MQ=255)
                                          cGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTAGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCtt                                                  >  1:1860611/1‑100 (MQ=255)
                                           gagaTGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTc                                                 <  1:2795271/100‑1 (MQ=255)
                                               tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                             >  2:1377970/1‑100 (MQ=255)
                                                 gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                           >  1:2787612/1‑100 (MQ=255)
                                                 gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                           >  1:2760124/1‑100 (MQ=255)
                                                  ttGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTaa                                          <  1:2669807/100‑1 (MQ=255)
                                                     tAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcc                                       <  2:2388526/100‑1 (MQ=255)
                                                      aaTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATccc                                      <  2:616005/100‑1 (MQ=255)
                                                      aaTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATcc                                       >  2:2781909/1‑99 (MQ=255)
                                                      aaTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATacc                                      >  1:3157915/1‑100 (MQ=255)
                                                          cACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGt                                   <  1:1735463/99‑1 (MQ=255)
                                                               tCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTg                             <  2:1454087/100‑1 (MQ=255)
                                                                cTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTgg                            >  2:1673416/1‑100 (MQ=255)
                                                                 tGTAACTGGAAGTTGGTGATGTATATCGACTCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCTCTTTAATCCCCGTCCACTgg                            <  1:2024975/99‑1 (MQ=255)
                                                                 tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                           >  2:1232383/1‑100 (MQ=255)
                                                                 tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                           >  2:1480593/1‑100 (MQ=255)
                                                                  gTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTgggn                          <  1:3465975/100‑2 (MQ=255)
                                                                  gTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGga                          >  2:2936060/1‑100 (MQ=255)
                                                                  gTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGga                          <  2:1731927/100‑1 (MQ=255)
                                                                   tAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgag                         >  2:268535/1‑100 (MQ=255)
                                                                     acgggAAGTTGGTGGTGTATATCGACGCGTCTTTCTTCTCAGATATAAACAGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCTTTCACTGGgaga                        >  1:687413/4‑99 (MQ=255)
                                                                     aCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaa                       >  1:315579/1‑100 (MQ=255)
                                                                     aCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaa                       >  1:1191286/1‑100 (MQ=255)
                                                                     aCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaa                       >  1:913255/1‑100 (MQ=255)
                                                                        ggAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaaaaa                    >  2:1142687/1‑100 (MQ=255)
                                                                        ggAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaaaaa                    >  2:1766205/1‑100 (MQ=255)
                                                                        ggAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaaaaa                    >  1:651986/1‑100 (MQ=255)
                                                                         gAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAg                   <  2:1685600/100‑1 (MQ=255)
                                                                         gAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAg                   <  2:1126137/100‑1 (MQ=255)
                                                                         gAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAg                   <  1:1387171/100‑1 (MQ=255)
                                                                           aGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGcc                 <  2:2218675/100‑1 (MQ=255)
                                                                           aGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGcc                 <  2:1715247/100‑1 (MQ=255)
                                                                           aGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGcc                 <  1:2163318/100‑1 (MQ=255)
                                                                             ttGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCtn               <  1:1377950/100‑2 (MQ=255)
                                                                               ggagATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTAAAAAGAAAGAAAAAGCCTaag             <  1:219903/97‑1 (MQ=255)
                                                                               ggTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCTTCCACTGGGCGAAAACGCCTaag             >  1:2783597/1‑100 (MQ=255)
                                                                                gTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaaga            <  2:98346/100‑1 (MQ=255)
                                                                                 tGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaagaa           <  2:2760168/100‑1 (MQ=255)
                                                                                    tGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGt         >  2:481721/1‑99 (MQ=255)
                                                                                    tGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTc        >  2:2516507/1‑100 (MQ=255)
                                                                                     gTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCa       >  2:3256122/1‑100 (MQ=255)
                                                                                       atatCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCAtt     >  1:946874/1‑100 (MQ=255)
                                                                                       atatCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCAtt     >  1:999493/1‑100 (MQ=255)
                                                                                        tatCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTg    >  1:2172776/1‑100 (MQ=255)
                                                                                        tatCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTg    >  2:882333/1‑100 (MQ=255)
                                                                                        tatCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTg    >  2:2543592/1‑100 (MQ=255)
                                                                                         atCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTgg   <  1:2786272/100‑1 (MQ=255)
                                                                                         atCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTgg   >  1:2578856/1‑100 (MQ=255)
                                                                                          tCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGc  <  2:3210673/100‑1 (MQ=255)
                                                                                          tCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGc  <  1:1813585/100‑1 (MQ=255)
                                                                                          tCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGc  >  1:1574480/1‑100 (MQ=255)
                                                                                             |                                                                                                
CGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGC  >  NC_000913/851001‑851190

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