Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,577,694	G→A	H89Y (CAT→TAT)	mcrC ←	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,577,694	0	G	A	100.0%	236.1 / NA	63	H89Y (CAT→TAT)	mcrC	5‑methylcytosine‑specific restriction enzyme McrBC, subunit McrC
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/33);  total (30/33)

GGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATAT  >  NC_000913/4577546‑4577840                                                                                                                                                     |                                                                                                                                                   ggTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAAt                                                                                                                                                  <  1:133724/151‑1 (MQ=255)    aaTTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTaa                                                                                                                                               <  1:499081/151‑1 (MQ=255)     aTTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTa                                                                                                                                                >  2:529541/1‑149 (MQ=255)            ataAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGa                                                                                                                                         >  1:742458/1‑149 (MQ=255)            ataAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGa                                                                                                                                         >  2:859871/1‑149 (MQ=255)              aaaGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGcc                                                                                                                                     <  2:515866/151‑1 (MQ=255)              aaaGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGcc                                                                                                                                     >  2:427115/1‑151 (MQ=255)                   gAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGt                                                                                                                                 <  1:849438/150‑1 (MQ=255)                      cGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTg                                                                                                                             <  1:198688/151‑1 (MQ=255)                      cGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTg                                                                                                                             <  2:325104/151‑1 (MQ=255)                      cGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTg                                                                                                                             <  1:28408/151‑1 (MQ=255)                          cTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGtttt                                                                                                                         >  1:626528/1‑151 (MQ=255)                                    gATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAActct                                                                                                               <  1:375317/151‑1 (MQ=255)                                     aTAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAActct                                                                                                               >  1:441992/1‑150 (MQ=255)                                       aGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTAt                                                                                                             <  1:568957/150‑1 (MQ=255)                                        gTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTc                                                                                                           >  2:335857/1‑151 (MQ=255)                                           gactttAACTTTTCATGCTTGATTAATATGGCTAATGTGCTTTTTATAATTTGGTTAGCCAGCGTGTGTTCATTAATCATATCTAAAGTACTGACGGTTTTCCCATAATTAAGATGAAAGCCGAGTATAGTTTTAGAAAACTCCATTc                                                                                                           >  2:665456/4‑148 (MQ=255)                                            aaTTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGccc                                                                                                       >  1:434405/1‑151 (MQ=255)                                              tttAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCtt                                                                                                     >  2:745193/1‑151 (MQ=255)                                               ttAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCttt                                                                                                    >  2:767765/1‑151 (MQ=255)                                               ttAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCtt                                                                                                     <  2:199543/150‑1 (MQ=255)                                                  aCTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGAt                                                                                                 >  1:225905/1‑151 (MQ=255)                                                  aCTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGAt                                                                                                 <  2:900916/151‑1 (MQ=255)                                                     tttCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATg                                                                                                <  2:263403/149‑1 (MQ=255)                                                     tttCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGcc                                                                                              >  1:83499/1‑151 (MQ=255)                                                       tCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGcc                                                                                              >  1:610383/1‑149 (MQ=255)                                                            cTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATg                                                                                       >  2:301538/1‑151 (MQ=255)                                                            cTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATg                                                                                       >  2:137328/1‑151 (MQ=255)                                                               aattaatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAt                                                                                     <  1:816895/150‑1 (MQ=255)                                                                  taatATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGAtctc                                                                                  <  1:541561/150‑1 (MQ=255)                                                                               gtgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTg                                                                     <  1:444425/150‑1 (MQ=255)                                                                               gtgCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGt                                                                    >  2:451136/1‑151 (MQ=255)                                                                                 gCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTa                                                                   >  2:46596/1‑150 (MQ=255)                                                                                   tttttATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtc                                                                >  1:404737/1‑151 (MQ=255)                                                                                      ttATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtca                                                               >  1:678044/1‑149 (MQ=255)                                                                                        ataATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagct                                                           <  1:259641/151‑1 (MQ=255)                                                                                          aaTTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctc                                                          <  1:504273/150‑1 (MQ=255)                                                                                             tCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAAtcaagctcaag                                                       <  1:523667/150‑1 (MQ=255)                                                                                                   aGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCCTAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCgg                                                >  1:499065/1‑151 (MQ=255)                                                                                                   aGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCgg                                                >  2:197284/1‑151 (MQ=255)                                                                                                       aGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTg                                            >  2:316928/1‑151 (MQ=255)                                                                                                          gtgtCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGaaa                                         <  1:818503/151‑1 (MQ=255)                                                                                                            gtCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAg                                        <  2:572779/150‑1 (MQ=255)                                                                                                                tCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCACGCTCAAGCCCTCGGCGTGAAAGCTGt                                    >  2:79606/1‑150 (MQ=255)                                                                                                                 cATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTa                                   <  2:603791/150‑1 (MQ=255)                                                                                                                 cATTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTa                                   >  2:181897/1‑150 (MQ=255)                                                                                                                  aTTAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaaa                                 <  1:569850/151‑1 (MQ=255)                                                                                                                    tAAGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTaaaa                                <  2:96352/150‑1 (MQ=255)                                                                                                                      aGCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAAccc                             >  2:615246/1‑151 (MQ=255)                                                                                                                       gCATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAAcccc                            <  2:671991/151‑1 (MQ=255)                                                                                                                        cATATCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCt                           <  2:320701/151‑1 (MQ=255)                                                                                                                         atatCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtt                          >  2:529477/1‑151 (MQ=255)                                                                                                                         atatCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtt                          <  1:403683/151‑1 (MQ=255)                                                                                                                           atCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCttt                         <  2:152533/150‑1 (MQ=255)                                                                                                                           atCAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCttt                         >  2:517797/1‑150 (MQ=255)                                                                                                                             cAAAAGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttat                       <  1:504338/150‑1 (MQ=255)                                                                                                                               aaaGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCtttattt                     >  2:840694/1‑150 (MQ=255)                                                                                                                                 aGTACTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAat                  <  2:286611/151‑1 (MQ=255)                                                                                                                                    aCTGACGGTTTTCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAataca               <  1:880941/151‑1 (MQ=255)                                                                                                                                    aCTGACGGTTTTCCCATAATTAAGATGGAAGCCGCCTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCCGTGTTAGGATTGTTATCAAGCTCAAGCCCTCGGCGTGATAGCTGTAAAACCCCTTTATTTAatac                <  2:507940/150‑1 (MQ=255)                                                                                                                                          ggTTTTCCCATAATTAAGATGGAATCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATAcccc         <  2:868498/151‑1 (MQ=255)                                                                                                                                              ttCCCATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAgg      <  2:174660/150‑1 (MQ=255)                                                                                                                                                  cATAATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTAGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGatat  >  2:711407/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGTAATTTTCTATAAAGTGAACGAGCTTCATCTCTGATAGTTGAATTTAACTTTTCATGCTTAATTAATATGGCTAATGTGCTTTTTATAATTCGGTTAGCCAGCGTGTCTTCATTAAGCATATCAAAAGTACTGACGGTTTTCCCATGATTAAGATGGAAGCCGCGTATTGTTTTAGCAAACTCTATTCGCCCTTTGATGCCAGGAATGATCTCGGTGTTAGGATTGTAATCAAGCTCAAGCCCTCGGCGTGAAAGCTGTAAAACCCCTTTATTTAATACATACCCCAGGATAT  >  NC_000913/4577546‑4577840

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