Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	934,839	T→C	V798A (GTT→GCT)	ftsK →	cell division protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	934,839	0	T	C	97.9%	283.2 / ‑3.4	97	V798A (GTT→GCT)	ftsK	cell division protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (49/46);  total (51/46)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.07e-01

ACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATG  >  NC_000913/934750‑934928                                                                                          |                                                                                          acaacaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTgc                                                                                        >  1:1160316/1‑93 (MQ=255)    acaGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGcgcc                                                                                     >  1:7411/1‑93 (MQ=255)     caGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGcgccg                                                                                    >  2:909071/1‑93 (MQ=255)        cgggTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcag                                                                                 <  2:3258150/91‑1 (MQ=255)          ggTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagca                                                                               >  2:1041941/1‑93 (MQ=255)          ggTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGGTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagca                                                                               >  2:197609/1‑93 (MQ=255)           gTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcag                                                                              >  2:2199418/1‑93 (MQ=255)           gTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcag                                                                              >  1:2997912/1‑93 (MQ=255)            ttGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcagc                                                                             >  2:2762353/1‑93 (MQ=255)             tGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAAAAACCGGCTGCGCCgcagcagca                                                                            <  1:3403529/93‑1 (MQ=255)             tGCGCCGCAGCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCgcagcagca                                                                            <  2:4178598/93‑1 (MQ=255)                gccgcAGCCACAAGATCAGCCGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTa                                                                         <  2:3524502/93‑1 (MQ=255)                cccgcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTa                                                                         <  2:2430927/92‑1 (MQ=255)                 ccgcAGCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTAt                                                                        <  2:602430/93‑1 (MQ=255)                  cgcAGCAACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATc                                                                       <  1:3193317/93‑1 (MQ=255)                   gcAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATca                                                                      >  2:1906977/1‑93 (MQ=255)                   cccaccACAATTCCAGCAGCCGCAACAACCGCTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATca                                                                      <  1:259110/89‑1 (MQ=17)                     acccACAATATCAGCAGCCGCAACAGCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagc                                                                    <  2:4294317/91‑1 (MQ=255)                     aGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagc                                                                    >  1:4683146/1‑93 (MQ=255)                     aGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagc                                                                    >  2:3898577/1‑93 (MQ=255)                     aGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCAGCTGCGCCGCAGCAGCAGTATcagc                                                                    >  1:871481/1‑93 (MQ=255)                     aGCCACAATATCAACAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagc                                                                    >  1:871482/1‑93 (MQ=255)                      gCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagca                                                                   >  2:2639772/1‑93 (MQ=255)                      gCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagca                                                                   >  1:3180467/1‑93 (MQ=255)                      gCCACAATACCACCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATcagca                                                                   <  2:67262/93‑1 (MQ=255)                           aaTATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCcac                                                              >  2:198862/1‑93 (MQ=255)                           aaTATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCcac                                                              <  1:1376974/93‑1 (MQ=255)                             tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaa                                                            <  1:4384525/93‑1 (MQ=255)                             tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaa                                                            >  2:1774727/1‑93 (MQ=255)                             tatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaa                                                            >  2:3007434/1‑93 (MQ=255)                             gatCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaa                                                            <  2:242635/92‑1 (MQ=255)                             aatCACCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaa                                                            <  1:2053653/92‑1 (MQ=255)                              atCAGCAGCCGCACCAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCacaac                                                           <  1:1020340/93‑1 (MQ=255)                                caacagcCGCAACAACGCGTTGCCCCACAGCAGCAATATCAGCAGCCGCAACAACCCGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAg                                                         <  1:2671589/90‑1 (MQ=14)                                 agcagcCGCACCACCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACATCCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGc                                                        <  2:2973490/93‑1 (MQ=17)                                  gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGcc                                                       <  2:1351983/93‑1 (MQ=255)                                  gcagcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGcc                                                       >  2:3617974/1‑93 (MQ=255)                                   cagcCGCAACACCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCa                                                      <  2:1500966/93‑1 (MQ=255)                                     gcCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGt                                                    <  1:1351979/93‑1 (MQ=255)                                     cccGCAACAACCGATTGCCCACCAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGt                                                    <  1:1351980/92‑1 (MQ=14)                                      ccGCAACACCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                   <  2:2351877/93‑1 (MQ=255)                                      ccGCAACAACCTGTTGCGCCACACAAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                   <  2:2626429/93‑1 (MQ=17)                                      ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                   <  1:4638655/93‑1 (MQ=255)                                      ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                   <  1:2974730/93‑1 (MQ=255)                                      ccGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGtt                                                   <  2:795706/93‑1 (MQ=255)                                         caacaaCCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcg                                                >  2:2837030/1‑93 (MQ=255)                                          aacaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcgc                                               >  2:4035283/1‑93 (MQ=255)                                          aacaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTgcgc                                               >  2:873448/1‑93 (MQ=255)                                           acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACACCAGCCAGTTGCGcc                                              >  1:1133494/1‑93 (MQ=255)                                           acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGcc                                              >  1:2275833/1‑93 (MQ=255)                                           acaacCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGcc                                              >  2:1715761/1‑93 (MQ=255)                                              acCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCcaca                                           >  2:4297669/1‑93 (MQ=255)                                               ccGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACaa                                          >  1:1702679/1‑93 (MQ=255)                                                cGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAc                                         <  1:996554/93‑1 (MQ=255)                                                  gtcgcgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCg                                       <  2:323870/90‑1 (MQ=255)                                                  gatGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGACAGTTGCGCCACAACCg                                       >  2:3992381/3‑93 (MQ=255)                                                  gTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCg                                       <  1:1733255/93‑1 (MQ=255)                                                    tGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                     <  1:807684/93‑1 (MQ=255)                                                    tGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCa                                     <  2:4428358/93‑1 (MQ=255)                                                     tcgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAg                                    <  2:1682170/92‑1 (MQ=255)                                                      cgcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAgg                                   >  2:2201637/1‑93 (MQ=255)                                                       gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCCGGa                                  >  1:3195678/1‑93 (MQ=255)                                                       gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGa                                  >  1:1823170/1‑93 (MQ=255)                                                       gcCACAGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAAACGCAGGa                                  >  2:4071052/1‑93 (MQ=255)                                                         cacaGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATa                                >  1:1485108/1‑93 (MQ=255)                                                          ccagcagCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAc                               <  1:3096469/92‑1 (MQ=255)                                                          acaGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACACCAGCCAGCTCTGCCACAACCGCAGGATAc                               >  1:4234441/1‑93 (MQ=18)                                                          acaGCAGCAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAc                               >  2:2354890/1‑93 (MQ=255)                                                            agcagcAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                             >  2:3284886/1‑93 (MQ=255)                                                            agcagcAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                             <  2:1755827/93‑1 (MQ=255)                                                            agcagcAATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATAccc                             <  2:1547446/93‑1 (MQ=255)                                                             gcagcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCt                            <  1:3494866/93‑1 (MQ=255)                                                               agcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGc                          >  1:152886/1‑93 (MQ=255)                                                               agcaATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGc                          >  2:683606/1‑93 (MQ=255)                                                                 caATATCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCtt                        <  2:4382959/93‑1 (MQ=255)                                                                    tatCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCAt                     >  2:3590001/1‑93 (MQ=255)                                                                      tCAGCAGCCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATcc                   >  2:3631074/1‑93 (MQ=255)                                                                         gcagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCt                >  2:1237417/1‑93 (MQ=255)                                                                         gcagcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCt                >  2:3295086/1‑93 (MQ=255)                                                                           agcCGCAACAACCGGCTGCGTCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTTCTTCATCCCCTGt              >  2:1288138/1‑93 (MQ=25)                                                                           agcCGCAACAACCGGCTGCGCCTCAGCAGCAATATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt              >  1:845940/1‑93 (MQ=25)                                                                           agcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt              <  1:2985290/93‑1 (MQ=255)                                                                           agcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGt              >  2:2359142/1‑93 (MQ=255)                                                                            gcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGtt             <  1:2337387/93‑1 (MQ=255)                                                                            gcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGtt             <  1:1683164/93‑1 (MQ=255)                                                                            gcCGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGtt             >  1:1862755/1‑93 (MQ=255)                                                                             ccGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGCTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTg            >  1:1535871/1‑93 (MQ=255)                                                                              cGCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGa           <  1:3003069/93‑1 (MQ=255)                                                                               gCAACAACCGGCTGCGCCGCAGCGGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt          >  1:4183402/1‑93 (MQ=255)                                                                               gCAACAACCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt          <  1:1171590/93‑1 (MQ=255)                                                                               ccaacaaCCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGAt          <  1:1551823/92‑1 (MQ=255)                                                                                caacaaCCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATg         <  1:2871551/93‑1 (MQ=255)                                                                                   caacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGt      <  1:725376/93‑1 (MQ=255)                                                                                   caacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGt      <  1:4714178/93‑1 (MQ=255)                                                                                    aacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa     >  2:454533/1‑93 (MQ=255)                                                                                    aacCGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTa     >  2:3658258/1‑93 (MQ=255)                                                                                    aacCGGCTGCGCCGCAGCAGCAGCATCAGCAGCCCCAACGGCCAGTTGCGCCACAACCGCAGGATAACCTGCTTCCTCCGCTGTTGATGCGTa     >  1:2675566/1‑93 (MQ=255)                                                                                       cGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATg  <  2:1689312/93‑1 (MQ=255)                                                                                       cGGCTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATg  >  1:1169626/1‑93 (MQ=255)                                                                                          |                                                                                          ACAACAGCCGGTTGCGCCGCAGCCACAATATCAGCAGCCGCAACAACCGGTTGCGCCACAGCAGCAATATCAGCAGCCGCAACAACCGGTTGCGCCGCAGCAGCAGTATCAGCAGCCACAACAGCCAGTTGCGCCACAACCGCAGGATACCCTGCTTCATCCGCTGTTGATGCGTAATG  >  NC_000913/934750‑934928

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