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_000913	851,094	0	G	A	93.8%	285.5 / 6.5	96	P502S (CCG→TCG)	opgE	OPG biosynthetic transmembrane phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base G (1/5);  new base A (40/50);  total (41/55)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGC  >  NC_000913/850997‑851190                                                                                                  |                                                                                                 ggggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGGGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGAcc                                                                                                <  1:1803645/100‑1 (MQ=255) ggggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGAcc                                                                                                >  1:203449/1‑100 (MQ=255) ggggCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGCTGGTGATGTATATCGAcc                                                                                                <  2:414010/100‑1 (MQ=255)       gATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCtt                                                                                          <  1:918261/100‑1 (MQ=255)       gATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCtt                                                                                          <  2:1770815/100‑1 (MQ=255)       gATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCtt                                                                                          <  2:185467/100‑1 (MQ=255)        aTTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCttt                                                                                         >  1:839217/1‑100 (MQ=255)         tttttGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTttc                                                                                        <  1:100774/100‑1 (MQ=255)          ttttGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTttct                                                                                       <  1:107262/100‑1 (MQ=255)             tGTTTTTAAGGTTTTGGGTCGAAAATATCGGTTCCGAGATGGTGGTAATCCACCTTCTGTAACTGGAAGGGCTTGATGTAGATCGACCCTGCTTtcttct                                                                                    >  1:1134040/1‑100 (MQ=255)             tGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTtcttct                                                                                    >  2:1530562/1‑100 (MQ=255)             tGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTtcttct                                                                                    <  2:598623/100‑1 (MQ=255)              gTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTtctc                                                                                   >  2:240633/1‑100 (MQ=255)               tttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCa                                                                                  <  1:805972/100‑1 (MQ=255)               tttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCa                                                                                  <  1:755189/100‑1 (MQ=255)               tttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCa                                                                                  <  2:151774/100‑1 (MQ=255)                ttttAAGGTTTGGGCTCGAAAATATCGGTTCCGTGATGGTTGGAATCCACCTTCTGTAACTTGAAGTTGGTGATGTATATCGACCCGGCTTTCTtctctg                                                                                 <  1:156156/100‑3 (MQ=255)                ttttAAGGTTTCGGGTCGAAAATATCGTTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAg                                                                                 <  2:1928644/100‑1 (MQ=255)                ttttAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAg                                                                                 >  2:291593/1‑100 (MQ=255)                ttttAAGGGTTGGGGTCGATAAGATCGGTGCCGAGGTGGTTGTAAGCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCGCGGCTTTCTTCTCAg                                                                                 <  1:1673921/100‑1 (MQ=255)                 tttAAGGTTTCGGGTCGAAAATATCGGTTCCTAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATTTATATCGACCCGGCTTTCTTCTCAGa                                                                                <  1:1894046/100‑1 (MQ=255)                   tAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGata                                                                              <  1:348980/100‑1 (MQ=255)                   tAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGata                                                                              <  1:1861405/100‑1 (MQ=255)                   tAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGata                                                                              >  1:1871852/1‑100 (MQ=255)                    aaGGTTTCGGGTCGAAAATATCGGGTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGatat                                                                             <  1:905047/100‑1 (MQ=255)                         ttCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACg                                                                        <  2:979472/100‑1 (MQ=255)                         ttCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACg                                                                        >  1:1398000/1‑100 (MQ=255)                          tCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACgg                                                                       <  2:689191/100‑1 (MQ=255)                           cGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACggg                                                                      <  2:1056540/100‑1 (MQ=255)                             ggTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTa                                                                    <  1:732183/100‑1 (MQ=255)                             ggTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTa                                                                    <  2:385597/100‑1 (MQ=255)                              gTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTAt                                                                   >  1:987709/1‑100 (MQ=255)                                cGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTAttt                                                                 >  2:807722/1‑100 (MQ=255)                                 gAAAATATCGGTTCCGAGATGGATGTAATCAACGTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAGCGGGTATTTg                                                                >  1:1275621/1‑100 (MQ=255)                                    aaTATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATg                                                             <  1:1541550/100‑1 (MQ=255)                                     atatCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGt                                                            <  1:1106216/100‑1 (MQ=255)                                      tatCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGtt                                                           >  2:192901/1‑100 (MQ=255)                                       atCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTa                                                          <  1:806048/100‑1 (MQ=255)                                          ggTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAtt                                                       >  2:414621/1‑100 (MQ=255)                                          ggTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAtt                                                       >  2:1084173/1‑100 (MQ=255)                                           gTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCGGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAAttt                                                      >  2:1462214/1‑100 (MQ=255)                                              ccGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCt                                                   <  1:211862/100‑1 (MQ=255)                                               cGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTcc                                                    >  1:1233708/1‑98 (MQ=255)                                                gagaTGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTc                                                 >  2:191997/1‑100 (MQ=255)                                                gagaTGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCATTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTc                                                 >  2:824337/1‑100 (MQ=255)                                                 agaTGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCg                                                >  1:150933/1‑100 (MQ=255)                                                    tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                             >  2:657976/1‑100 (MQ=255)                                                    tGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCtt                                             >  2:605426/1‑100 (MQ=255)                                                     ggTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCttt                                            <  2:1407592/100‑1 (MQ=255)                                                      gTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTa                                           <  2:203450/100‑1 (MQ=255)                                                       ttGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTaa                                          >  2:1432673/1‑100 (MQ=255)                                                        tGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAAt                                         <  2:1680563/100‑1 (MQ=255)                                                        tGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAAt                                         >  1:1607765/1‑100 (MQ=255)                                                         gTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATc                                        >  1:697691/1‑100 (MQ=255)                                                           aaTCCACCTTCTGTAACTGGAAGTTGTTGATGTATATCGGCCCGGCTTTCTTCTCAGATATGAACGGGTATTTGATGTTAATTTCCTTCGCTTTAAt                                         <  1:1106501/97‑1 (MQ=255)                                                           aaTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATccc                                      >  2:329138/1‑100 (MQ=255)                                                           aaTCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATccc                                      >  2:1591634/1‑100 (MQ=255)                                                             tCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCg                                    <  1:853106/100‑1 (MQ=255)                                                             tCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCg                                    <  2:1295751/100‑1 (MQ=255)                                                                 ccTTCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCa                                <  1:1699419/100‑1 (MQ=255)                                                                   ttCTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACt                              >  2:220441/1‑100 (MQ=255)                                                                     cTGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTgg                            <  1:1941138/100‑1 (MQ=255)                                                                      tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                           <  2:697692/100‑1 (MQ=255)                                                                      tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                           <  2:580620/100‑1 (MQ=255)                                                                      tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                           <  1:1306517/100‑1 (MQ=255)                                                                      tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                           >  2:1087596/1‑100 (MQ=255)                                                                      tGTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTggg                           >  2:1070862/1‑100 (MQ=255)                                                                       gTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGga                          >  2:1685743/1‑100 (MQ=255)                                                                       gTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGga                          <  1:1922042/100‑1 (MQ=255)                                                                       gTAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGga                          >  1:956182/1‑100 (MQ=255)                                                                       cgatcTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGga                          <  1:824949/96‑1 (MQ=255)                                                                        tAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgag                         <  2:1883558/100‑1 (MQ=255)                                                                        tAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgag                         <  1:79355/100‑1 (MQ=255)                                                                        tAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgag                         <  1:345333/100‑1 (MQ=255)                                                                        tAACTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgag                         <  1:1480499/100‑1 (MQ=255)                                                                         aaCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgaga                        >  1:1448454/1‑100 (MQ=255)                                                                         aaCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgaga                        <  2:1871864/100‑1 (MQ=255)                                                                          aCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaa                       >  2:243688/1‑100 (MQ=255)                                                                          aCTGGAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaa                       <  1:1591623/100‑1 (MQ=255)                                                                             ggAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGgaagatcg                    >  2:824950/1‑94 (MQ=255)                                                                             ggAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGaaaaa                    >  1:1421090/1‑100 (MQ=255)                                                                              gAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAg                   <  1:1256556/100‑1 (MQ=255)                                                                              gAAGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAg                   <  2:143868/100‑1 (MQ=255)                                                                               aaGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGc                  <  2:402886/100‑1 (MQ=255)                                                                                agtttttGATGTATATCGACCCGGCGTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGcc                 <  2:1134041/94‑1 (MQ=255)                                                                                aGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGcc                 <  2:1267731/100‑1 (MQ=255)                                                                                aGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGcc                 >  2:1579529/1‑100 (MQ=255)                                                                                aGTTGGTGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGcc                 <  2:402432/100‑1 (MQ=255)                                                                                     gTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaaga            <  2:728236/100‑1 (MQ=255)                                                                                      tGATGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTaagaa           >  1:546786/1‑100 (MQ=255)                                                                                         tGTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTc        >  2:718428/1‑100 (MQ=255)                                                                                         tGTATATCGACCCGGCTTTCTTCTCAGATATAAAAGGGTATTTGATGTTAATTTCCTTTGCTTTAATTCCCGTCCACTGGGAGAAAAAGCCTaaga            >  1:1842048/1‑96 (MQ=255)                                                                                          gTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGtaa       >  1:1635649/1‑98 (MQ=255)                                                                                          gTATATCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCa       <  2:839218/100‑1 (MQ=255)                                                                                              atCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTgg   <  1:809437/100‑1 (MQ=255)                                                                                               tCGACCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGc  >  2:28410/1‑100 (MQ=255)                                                                                                  |                                                                                                 GGGGCGGATTTTTGTTTTTAAGGTTTCGGGTCGAAAATATCGGTTCCGAGATGGTTGTAATCCACCTTCTGTAACTGGAAGTTGGTGATGTATATCGGCCCGGCTTTCTTCTCAGATATAAACGGGTATTTGATGTTAATTTCCTTCGCTTTAATCCCCGTCCACTGGGAGAAAAAGCCTAAGAAGTCATTGGC  >  NC_000913/850997‑851190

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