Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,494,973	G→A	*92* (TAG→TAA)	ypdI →	putative lipoprotein involved in colanic acid biosynthesis

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,494,973	0	G	A	100.0%	270.3 / NA	82	*92* (TAG→TAA)	ypdI	putative lipoprotein involved in colanic acid biosynthesis
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/45);  total (37/45)

CACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAA  >  NC_000913/2494839‑2495105                                                                                                                                       |                                                                                                                                     cacTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTc                                                                                                                                  <  2:1765611/139‑1 (MQ=255)      tACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTca                                                                                                                             >  1:63589/1‑139 (MQ=255)         aaaaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                                                          <  1:666199/139‑1 (MQ=255)         aaaaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                                                          <  2:1437771/139‑1 (MQ=255)           aaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATgg                                                                                                                        <  2:1341402/139‑1 (MQ=255)           aaCTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATgg                                                                                                                        >  2:470100/1‑139 (MQ=255)              ttCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                     >  1:388020/1‑139 (MQ=255)              ttCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                     >  2:2055717/1‑139 (MQ=255)                cAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                                                          >  1:1718888/1‑132 (MQ=255)                cAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACAt                                                                                                                          <  2:1718888/132‑1 (MQ=255)                  gCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGa                                                                                                                 <  1:627678/139‑1 (MQ=255)                       cacaGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAAt                                                                                                            >  2:598737/1‑139 (MQ=255)                        acaGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATc                                                                                                           <  1:2055717/139‑1 (MQ=255)                             cGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTc                                                                                                      <  2:1203485/139‑1 (MQ=255)                               aTTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAA‑CTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctga                                                                                                   <  1:1808855/139‑1 (MQ=255)                                 tttGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgat                                                                                                  >  2:274216/1‑139 (MQ=255)                                 tttGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTc                                                                                                      <  2:1432713/135‑1 (MQ=255)                                 tttGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTc                                                                                                      >  1:1432713/1‑135 (MQ=255)                                    gTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                               <  1:274216/139‑1 (MQ=255)                                     tttAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                              <  2:1831582/139‑1 (MQ=255)                                     tttAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                              <  2:957127/139‑1 (MQ=255)                                             gagagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAac                                                                                      <  2:22786/139‑1 (MQ=255)                                               gagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAacac                                                                                    <  2:1049436/139‑1 (MQ=255)                                               gagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAacac                                                                                    >  2:944347/1‑139 (MQ=255)                                               gagaAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAacac                                                                                    >  1:702077/1‑139 (MQ=255)                                                    gaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAac                                                                                      <  2:320398/132‑1 (MQ=255)                                                    gaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAac                                                                                      >  1:320398/1‑132 (MQ=255)                                                    gaTGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGcc                                                                               >  2:1546418/1‑139 (MQ=255)                                                        cgcCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                     <  1:1799641/97‑1 (MQ=255)                                                        cgcCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAg                                                                                                                     >  2:1799641/1‑97 (MQ=255)                                                        cgcCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCaaa                                                                           >  2:1839164/1‑139 (MQ=255)                                                         gcCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAAt                                                                          >  1:933389/1‑139 (MQ=255)                                                          ccaccaTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATg                                                                         >  2:1573885/1‑139 (MQ=255)                                                            accaTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                       <  2:388020/139‑1 (MQ=255)                                                             ccaTTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGc                                                                      <  2:933389/139‑1 (MQ=255)                                                                      cTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATAAGCCATCCTGAc                                                             <  2:1582386/139‑1 (MQ=255)                                                                        cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCttt                                                                                                                               >  2:1219322/1‑71 (MQ=255)                                                                        cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCttt                                                                                                                               <  1:1219322/71‑1 (MQ=255)                                                                        cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGaa                                                                                                                >  2:349796/1‑86 (MQ=255)                                                                        cacaGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGaa                                                                                                                <  1:349796/86‑1 (MQ=255)                                                                                aaaTGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGc                                                   <  1:470100/139‑1 (MQ=255)                                                                                 aaTGCCGCCCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGcc                                                  <  2:1718690/139‑1 (MQ=255)                                                                                     ccgccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGa                                                                                                                 >  2:465557/1‑72 (MQ=255)                                                                                     ccgccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGa                                                                                                                 <  1:465557/72‑1 (MQ=255)                                                                                     ccgccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGc                                                                                <  2:544403/105‑1 (MQ=255)                                                                                     ccgccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGc                                                                                >  1:544403/1‑105 (MQ=255)                                                                                       gccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                              <  1:463447/89‑1 (MQ=255)                                                                                       gccCCCTTTTCCGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgctg                                                                                              >  2:463447/1‑89 (MQ=255)                                                                                                 ccGACAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCAt                                  <  1:1573885/139‑1 (MQ=255)                                                                                                     cAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATaa                                >  1:1773728/1‑137 (MQ=255)                                                                                                     cAAAGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATaa                                <  2:1773728/137‑1 (MQ=255)                                                                                                       aaGGCGAATCATCTTATTCATCCTGTGCGATAACGTCttt                                                                                                                               >  2:401267/1‑40 (MQ=37)                                                                                                       aaGGCGAATCATCTTATTCATCCTGTGCGATAACGTCttt                                                                                                                               <  1:401267/40‑1 (MQ=37)                                                                                                        aGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTc                                                                                                      >  2:1405653/1‑64 (MQ=255)                                                                                                        aGGCGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTc                                                                                                      <  1:1405653/64‑1 (MQ=255)                                                                                                           cGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATgg                                                                                                                        >  2:563468/1‑43 (MQ=38)                                                                                                           cGAATCATCTTATTCATCCTGTGCGATAACGTCTTTCACATgg                                                                                                                        <  1:563468/43‑1 (MQ=37)                                                                                                                            ccTGTGCGATAACGTCTTTCACATGGCTGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATcc                                                                 <  1:1913755/81‑1 (MQ=255)                                                                                                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTc                                                                                                      >  1:824465/1‑44 (MQ=255)                                                                                                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTc                                                                                                      <  2:824465/44‑1 (MQ=255)                                                                                                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATcc                                                                 >  2:1913755/1‑81 (MQ=255)                                                                                                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCg                                                    >  2:115612/1‑94 (MQ=255)                                                                                                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCg                                                    <  1:115612/94‑1 (MQ=255)                                                                                                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAatcatc                     <  2:1813330/125‑1 (MQ=255)                                                                                                                            ccTGTGCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAatcatc                     >  1:1813330/1‑125 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                               >  2:1998003/1‑47 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                               >  1:149748/1‑47 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                               <  1:1998003/47‑1 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTctgatgct                                                                                               <  2:149748/47‑1 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                       <  2:323662/71‑1 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAg                                                                       >  1:323662/1‑71 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTg                                                <  2:1297801/94‑1 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTg                                                >  1:1297801/1‑94 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCa                                             <  1:638089/97‑1 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCa                                             >  2:638089/1‑97 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa   >  1:1660842/1‑139 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa   <  2:1660842/139‑1 (MQ=255)                                                                                                                                tgCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCa   <  2:702077/139‑1 (MQ=255)                                                                                                                                 gCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATAccc                                        <  1:1262060/101‑1 (MQ=255)                                                                                                                                 gCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCaa  >  2:600098/1‑139 (MQ=255)                                                                                                                                 gCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCATTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCaa  <  1:944347/139‑1 (MQ=255)                                                                                                                                 gCGATAACGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCACCCTGACCCATTACCGCCTGCCATAccc                                        >  2:1262060/1‑101 (MQ=255)                                                                                                                                       |                                                                                                                                     CACTCTACAAAACTTCAGCCATCACAGGCGATTTTGTTTAAGATGAGAGAAGATGCGCCACCATTAAACCTCACAGAAGAAATGCCGCCCCCTTTTCCGACAAAGGCGAATTATCTTATTCATCCTGTGCGATAGCGTCTTTCACATGGCAGGTGAACAATCATGTCTGATGCTGATGTGAACACCTGCCCAAATGAGCCATCCTGACCCACTACCGCCTGCCATACCCCACCATAAAGCGAATCATCAAGCCACAAGCCGGGGCAA  >  NC_000913/2494839‑2495105

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑