Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,875,282	A→G	R46R (CGT→CGC)	cysN ←	sulfate adenylyltransferase, subunit 1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,875,282	0	A	G	100.0%	281.5 / NA	77	R46R (CGT→CGC)	cysN	sulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/41);  total (36/41)

GTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTC  >  NC_000913/2875149‑2875417                                                                                                                                      |                                                                                                                                        gTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGAcc                                                                                                                                    >  2:1196259/1‑139 (MQ=255)    cgcgCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTGACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAAt                                                                                                                                 <  1:1034984/139‑1 (MQ=255)       gcTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCag                                                                                                                              <  2:1503539/139‑1 (MQ=255)                 cAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTg                                                                                                                    >  2:1456432/1‑139 (MQ=255)                        tCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCAt                                                                                                             >  2:1323075/1‑139 (MQ=255)                         ccaccaGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTGTTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                            <  2:1154044/139‑1 (MQ=255)                         ccaccaGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                            <  1:209920/139‑1 (MQ=255)                         ccaccaGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                            <  1:1456432/139‑1 (MQ=255)                            ccaGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGAc                                                                                                         <  1:1196259/139‑1 (MQ=255)                            ccaGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGAc                                                                                                         <  2:476318/139‑1 (MQ=255)                              agcagAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGc                                                                                                       <  2:683069/139‑1 (MQ=255)                                  gagCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTAcc                                                                                                   <  2:747838/139‑1 (MQ=255)                                   ggCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACAGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGAGTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGACGTCATCGACGCTACca                                                                                                  >  2:990461/2‑139 (MQ=255)                                                   tCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                            <  1:738342/113‑1 (MQ=255)                                                   tCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                            >  2:738342/1‑113 (MQ=255)                                                   tCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgca                                                                                  >  1:318436/1‑139 (MQ=255)                                                      cccTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgcagca                                                                               <  1:1476031/139‑1 (MQ=255)                                                         tGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGc                                                                            >  1:277518/1‑139 (MQ=255)                                                                cGTGACGCTTACTGTCGTTATGCAGCGCTGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCATCAGGCTTTTAtg                                                                     >  1:960494/1‑139 (MQ=255)                                                                     cGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATgttgtt                                                                >  2:1031130/1‑139 (MQ=255)                                                                       cTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgcag                                                                                 <  1:810325/120‑1 (MQ=255)                                                                       cTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgcag                                                                                 >  2:810325/1‑120 (MQ=255)                                                                                gTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgcagca                                                                               <  1:1006972/113‑1 (MQ=255)                                                                                gTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAatcatc                                                     <  1:969050/139‑1 (MQ=255)                                                                                gTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAatcatc                                                     <  2:318436/139‑1 (MQ=255)                                                                                gTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAatcatc                                                     >  2:969050/1‑139 (MQ=255)                                                                                gTTATGCAGCGATGAGAACTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgcagca                                                                               >  2:1006972/1‑113 (MQ=255)                                                                                   aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAg                                                  >  1:1414606/1‑139 (MQ=255)                                                                                   aTGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTG‑AGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAgg                                                 <  1:191914/139‑1 (MQ=255)                                                                                        gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgcagca                                                                               <  2:509156/105‑1 (MQ=255)                                                                                        gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACgcagca                                                                               >  1:509156/1‑105 (MQ=255)                                                                                        gCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTc                                             >  1:1474677/1‑139 (MQ=255)                                                                                         cGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCg                                            >  2:1501352/1‑139 (MQ=255)                                                                                           aTGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACca                                                                                                  >  2:47545/1‑83 (MQ=255)                                                                                           aTGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACca                                                                                                  <  1:47545/83‑1 (MQ=255)                                                                                           aTGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGAc                                          >  1:719391/1‑139 (MQ=255)                                                                                            tGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGAcg                                         <  2:1474677/139‑1 (MQ=255)                                                                                                agCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGAcgccgc                                     >  2:705887/1‑139 (MQ=255)                                                                                                   tGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGAcc                                                                                                                                    <  2:934550/41‑1 (MQ=255)                                                                                                   tGATCTTCGTACATTTGGCGGGTATCGTGCAGCAGGCGAcc                                                                                                                                    >  1:934550/1‑41 (MQ=38)                                                                                                    gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACg                                                                                                        <  1:847034/68‑1 (MQ=255)                                                                                                    gATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACg                                                                                                        >  2:847034/1‑68 (MQ=255)                                                                                                     aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCagag                                                                                                                            <  2:58141/47‑1 (MQ=255)                                                                                                     aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCagag                                                                                                                            >  1:58141/1‑47 (MQ=255)                                                                                                     aTCTTCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCa                                >  1:95496/1‑139 (MQ=255)                                                                                                        ttCGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTg                             <  1:1031130/139‑1 (MQ=255)                                                                                                          cGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                            >  1:532941/1‑58 (MQ=255)                                                                                                          cGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                            <  2:532941/58‑1 (MQ=255)                                                                                                          cGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGcagcag                                                                              >  1:636133/1‑88 (MQ=255)                                                                                                          cGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGcagcag                                                                              <  2:636133/88‑1 (MQ=255)                                                                                                          cGTAGATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGc                           >  2:615212/1‑139 (MQ=255)                                                                                                              gATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGAtt                       <  2:1414606/139‑1 (MQ=255)                                                                                                              gATTTGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCAGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGAtt                       <  2:960494/139‑1 (MQ=255)                                                                                                                tttGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATttg                     <  2:1333487/139‑1 (MQ=255)                                                                                                                 ttGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGc                                                                                                       >  1:1494222/1‑56 (MQ=255)                                                                                                                 ttGGCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGc                                                                                                       <  2:1494222/56‑1 (MQ=255)                                                                                                                   ggCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGAc                                                                                                         <  1:113187/52‑1 (MQ=255)                                                                                                                   ggCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGAc                                                                                                         >  2:113187/1‑52 (MQ=255)                                                                                                                   ggCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATttgttg                  <  1:705887/139‑1 (MQ=255)                                                                                                                   ggCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATttgttg                  <  2:719391/139‑1 (MQ=255)                                                                                                                   ggCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATttgttg                  <  1:1323075/139‑1 (MQ=255)                                                                                                                    gCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATttg                     >  2:1563708/1‑135 (MQ=255)                                                                                                                    gCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATttg                     >  1:1167890/1‑135 (MQ=255)                                                                                                                    gCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATttg                     <  2:1167890/135‑1 (MQ=255)                                                                                                                    gCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATttg                     <  1:1563708/135‑1 (MQ=255)                                                                                                                    gCGGGTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTtgttgt                 >  1:1481878/1‑139 (MQ=255)                                                                                                                      gggTATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGc               <  2:95496/139‑1 (MQ=255)                                                                                                                         tATCGTGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAg            >  1:354350/1‑139 (MQ=255)                                                                                                                              tGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGc               >  1:926993/1‑131 (MQ=255)                                                                                                                              tGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGc               <  2:926993/131‑1 (MQ=255)                                                                                                                              tGCAGCAGGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCgg       >  1:1614224/1‑139 (MQ=255)                                                                                                                                  gcagGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGttgttg                                                               <  2:610697/79‑1 (MQ=255)                                                                                                                                  gcagGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGttgttg                                                               >  1:610697/1‑79 (MQ=255)                                                                                                                                  gcagGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGtt   <  2:1481878/139‑1 (MQ=255)                                                                                                                                   cagGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAgg                                                                             <  1:806755/64‑1 (MQ=255)                                                                                                                                   cagGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAgg                                                                             >  2:806755/1‑64 (MQ=255)                                                                                                                                   cagGCGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTc  <  2:354350/139‑1 (MQ=255)                                                                                                                                      |                                                                                                                                        GTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTC  >  NC_000913/2875149‑2875417

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