Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,540,785	C→T	intergenic (‑1284/‑172)	nanC ← / → fimB	N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,540,785	0	C	T	100.0%	277.6 / NA	84	intergenic (‑1284/‑172)	nanC/fimB	N‑acetylnuraminic acid outer membrane channel protein/tyrosine recombinase/inversion of on/off regulator of fimA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (44/40);  total (44/40)

AACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAACGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAG  >  NC_000913/4540660‑4540918                                                                                                                              |                                                                                                                                      aaCATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTa                                                                                                                          >  1:387917/1‑139 (MQ=255)             gATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAAc                                                                                                                        <  2:1981809/129‑1 (MQ=255)             gATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAAc                                                                                                                        >  1:1981809/1‑129 (MQ=255)             gATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAAc                                                                                                                        >  1:1884901/1‑129 (MQ=255)             gATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAAc                                                                                                                        <  2:1884901/129‑1 (MQ=255)                 tataGATAAGATTAGTTGTATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCAt                                                                                                          >  2:1969918/1‑139 (MQ=255)                        aaGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGtta                                                                                                                                   <  2:1520712/107‑1 (MQ=255)                        aaGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGtta                                                                                                                                   >  1:1520712/1‑107 (MQ=255)                         aGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGt                                                                                                    <  1:205881/137‑1 (MQ=255)                         aGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGt                                                                                                    >  2:205881/1‑137 (MQ=255)                           aTTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCAc                                                                                                >  2:1442779/1‑139 (MQ=255)                                    aTTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAg                                                                                                     <  1:796536/125‑1 (MQ=255)                                    aTTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAg                                                                                                     >  2:796536/1‑125 (MQ=255)                                             gaGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAAc                                                                                                                        <  1:1652785/97‑1 (MQ=255)                                             gaGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAAc                                                                                                                        >  2:1652785/1‑97 (MQ=255)                                             gaGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTc                                                                                 <  2:1778937/136‑1 (MQ=255)                                             gaGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTc                                                                                 >  1:1778937/1‑136 (MQ=255)                                                ggTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTa                                                                                                                <  1:303999/102‑1 (MQ=255)                                                ggTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTa                                                                                                                >  2:303999/1‑102 (MQ=255)                                                          aGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCa                                                                                                 >  1:1626954/1‑107 (MQ=255)                                                          aGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCa                                                                                                 <  2:1626954/107‑1 (MQ=255)                                                           gATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTgg                                                                                   <  1:1951050/120‑1 (MQ=255)                                                           gATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTgg                                                                                   >  2:1951050/1‑120 (MQ=255)                                                               cGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCa                                                                                                 <  1:327372/102‑1 (MQ=255)                                                               cGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCa                                                                                                 >  2:327372/1‑102 (MQ=255)                                                                 tATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACa                                                                                                                       <  1:876788/78‑1 (MQ=255)                                                                 tATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACa                                                                                                                       >  2:876788/1‑78 (MQ=255)                                                                 tATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTa                                                                                                                >  1:1204926/1‑85 (MQ=255)                                                                 tATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTa                                                                                                                <  2:1204926/85‑1 (MQ=255)                                                                                     aaaGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACAt                                      >  1:1129967/1‑139 (MQ=255)                                                                                          tACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACgt                                                                        >  2:711885/1‑100 (MQ=255)                                                                                          tACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACgt                                                                        <  1:711885/100‑1 (MQ=255)                                                                                          tACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTta                                             >  2:2001441/1‑127 (MQ=255)                                                                                          tACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTta                                             <  1:2001441/127‑1 (MQ=255)                                                                                          tACATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCa                                 >  1:1047863/1‑139 (MQ=255)                                                                                           aCATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCg                                                   <  1:1151272/120‑1 (MQ=255)                                                                                           aCATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCg                                                   >  2:1151272/1‑120 (MQ=255)                                                                                            cATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAAGAGTTGCTACTCCATTTCAGt                                                                                                    <  2:1398176/70‑1 (MQ=255)                                                                                            cATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGt                                                                                                    >  1:1398176/1‑70 (MQ=255)                                                                                            cATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGtt                                                     <  1:1679723/117‑1 (MQ=255)                                                                                            cATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGtt                                                     >  2:1679723/1‑117 (MQ=255)                                                                                            cATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAgg                               >  1:27137/1‑139 (MQ=255)                                                                                            cATAGCATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAgg                               >  1:172078/1‑139 (MQ=255)                                                                                                 cATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACt                                                                                                              >  1:930493/1‑55 (MQ=255)                                                                                                 cATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACt                                                                                                              <  2:930493/55‑1 (MQ=255)                                                                                                 cATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGtt                                                     <  1:174144/112‑1 (MQ=255)                                                                                                 cATGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGtt                                                     >  2:174144/1‑112 (MQ=255)                                                                                                  aTGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCg                                                   >  1:1142298/1‑113 (MQ=255)                                                                                                  aTGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCg                                                   <  2:1142298/113‑1 (MQ=255)                                                                                                   tGCAAGAGCATGGCGTTTGTATGGCAGTGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCa                                                                                                 >  2:1124461/1‑66 (MQ=255)                                                                                                   tGCAAGAGCATGGCGTTTGTATGGCAGTGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCa                                                                                                 <  1:1124461/66‑1 (MQ=255)                                                                                                   tGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGt                                                                                                    <  1:1811744/63‑1 (MQ=255)                                                                                                   tGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGt                                                                                                    >  2:1811744/1‑63 (MQ=255)                                                                                                   tGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCa                                                                                                 >  1:1886758/1‑66 (MQ=255)                                                                                                   tGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCa                                                                                                 <  2:1886758/66‑1 (MQ=255)                                                                                                   tGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACt                                                                                               >  1:824628/1‑68 (MQ=255)                                                                                                   tGCAAGAGCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACt                                                                                               <  2:824628/68‑1 (MQ=255)                                                                                                          gCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGa                 <  2:1129967/139‑1 (MQ=255)                                                                                                          gCATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGa                 <  2:172078/139‑1 (MQ=255)                                                                                                           cATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTa                                                                                                                <  1:757150/43‑1 (MQ=255)                                                                                                           cATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTa                                                                                                                >  2:757150/1‑43 (MQ=255)                                                                                                           cATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCAc                                                                                                <  2:1227049/59‑1 (MQ=255)                                                                                                           cATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACgt                                                                        <  1:955311/83‑1 (MQ=255)                                                                                                           cATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACgt                                                                        >  2:955311/1‑83 (MQ=255)                                                                                                           cATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGat                >  1:990968/1‑139 (MQ=255)                                                                                                           cATGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGat                >  1:744769/1‑139 (MQ=255)                                                                                                           cATGGCGTTTGTATGGCAATGTTATTACAATTAACAGTTGCTACTCCATTTAAGTTCAc                                                                                                >  1:1227049/1‑59 (MQ=255)                                                                                                             tGGCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGatta              >  2:920448/1‑139 (MQ=255)                                                                                                              ggCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAAc                                        <  1:839432/112‑1 (MQ=255)                                                                                                              ggCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAAc                                        >  2:839432/1‑112 (MQ=255)                                                                                                              ggCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGatt               >  1:1855318/1‑137 (MQ=255)                                                                                                              ggCGTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGatt               <  2:1855318/137‑1 (MQ=255)                                                                                                                 gTTTGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGc          <  2:387917/139‑1 (MQ=255)                                                                                                                    tGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCAc                                                                                                <  1:1602889/50‑1 (MQ=255)                                                                                                                    tGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCAc                                                                                                >  2:1602889/1‑50 (MQ=255)                                                                                                                    tGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGt                                                                                  >  1:112501/1‑64 (MQ=255)                                                                                                                    tGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGt                                                                                  <  2:112501/64‑1 (MQ=255)                                                                                                                    tGTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTaa       >  1:1044931/1‑139 (MQ=255)                                                                                                                     gTATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAAc      >  2:681835/1‑139 (MQ=255)                                                                                                                      tATGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAAcc     <  2:1047863/139‑1 (MQ=255)                                                                                                                       aTGGCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAAccc    <  2:744769/139‑1 (MQ=255)                                                                                                                         ggCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAg  <  1:1442779/139‑1 (MQ=255)                                                                                                                         ggCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAg  <  2:27137/139‑1 (MQ=255)                                                                                                                         ggCAATGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAg  <  2:990968/139‑1 (MQ=255)                                                                                                                              |                                                                                                                                      AACATATCACATGATATATAGATAAGATTAGTTGCATTAATGATGAGGGTTATTATTAGATTCGTATCCGATTGATAAATATATAAAGGTACATAGCATGCAAGAGCATGGCGTTTGTATGGCAACGTTATTATAATTAACAGTTGCTACTCCATTTAAGTTCACTCAGAAGAACTGGTCCACTTACGTTAGTTATTAAGCAAACGTTCGCTTTTATAAACATAATCAGGATAAAAATGTTGGATTATTGCTAACCCAG  >  NC_000913/4540660‑4540918

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